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ABSTRACT 


A theoretical model for analysing the wing with 
submerged lifting fan is developed to determine the theo- 
retical aerodynamic characteristics of vertical/short 
take-off and landing (V/STOL) type aircraft. The model 
replaces the wing by a flat plate which is replaced by a 
net work of horse-shoe vortices distributed on its surface 
and trailing behind. For the ducted fan a series of singu- 
larities representing the fan disc, duct, centre-body and 
the wake are used. Flow tangency condition on the wing is 
used to determine the wing singularities and an iterative 
procedure is used for the ducted fan analysis. Both the 
effect of the wing on the fan and the fan on the wing are 
considered. The effect of the fan on the wing is obtained 
in a simplified manner by considering only the duct bound 
vorticity. The effects of all other ducted fan. singula- 
rities are assumed to be small in comparison. The theore- 
tical results obtained are compared with experiment. 

A computer program in Fortran IV for use on IBM 
7044 is developed for the computer simulation of the prob- 
lem. It uses an existing E'ASA program for the analysis of 
the ducted fan in cross-flow. 
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SYMBOLS 


For the Wing and the Combined Fan-Wing Model: 

p 

A aspect ratio,' 4b /S 

b semi- span of the wing 

c local chord 

aerodynamic mean chord, 


.4 

C F 

C i 

C, 


°m 

Cm 

Cl, 

■UA 

C L: 
G. 


Pi 


u 


~P 

°R 

E(k) 


/ 


9 

C dy 


propeller force coefficient, I/q^S 
local lift coefficient, Eq_. (2i9) 
conventional local lift coefficient, 
lift coefficient, Eq. (2*15) 
conventional lift coefficient 
local moment coefficient 
conventional local moment coefficient 
moment coefficient, Eq. (2*16) 
conventional moment coefficient, 


.JALft 
lco :s 


M( local) 

~ 2 
Oco * 


M ' 


Sc/ 


lower surface pressure coefficient, 
p l~ p o 


Qco 


upper surface pressure coefficient, 

P ~P 
_u ° 

coefficient of the pressure difference across 


the wing surface , C 
wing root chord 


n -CL 

Pi p 


u 


complete elliptic integral of the second kind, 
%/2 


]f 1~k 2 sin 2 9 


d9 


nondimen sional down- wash function, Eq. (2.6). 
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X 
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JTY 


Y' 


nondimensional down-wash function, Sq. (2.7) 

complete elliptic integral of the first kind, 

n/2 , 

f ae 

0 fl-k 2 sin 2 G 

lift 

moment 

matrix order or number of vortices 
static pressure 
free stream dynamic pressure 
radius vector 

vortex semi-span, non-dimensional with respect 

to Cj^ 

fan thrust 

induced velocity 

Forward velocity 

downwash 'velocity 

chordwise coordinate, non-dimensional with respect 
to or axial distance measured from duct c/2 
axial distance measured from duct c/2 
fraction of chord over which Cp is assumed 
constant, non-dimensions,! with respect to Cj^ 
chordwise location of the control point with 
respect to vortex centre, x -x v 

spanwise location of the control point with respect 
to vortex centre , y^-y^ 

spanwise location of the control point with res- 
pect to the image vortex centre, Yp+7 V 
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r 
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angle of incidence of the wing s radians 
angle of incidence of the fan with respect to 
the fan axis of .rotation;, radians 
circulation 
air density 


Matrices: 


[ ] 

[ 3 

< 1 
L j 

>C < 

l p'. 

W 


a square or rectangular matrix 
a diagonal matrix 
a column matrix 

matrix of pressure coefficients 

matrix of downwash influence coefficients 

matrix of downwash velocity at the control points 

matrix of distances (in root chords) of the 

vortex centre from the line about which the 

moment ? are to be determined. 

modulus 


Subscripts : 


f fan 

i matrix element identification by row 

0 matrix element identification by column or jet 

1 lower surface or local 

le leading edge 

o ' ambient 

p control point or propeller 

radial 


r 



V 


te trailing edge 
v vortex centre 
w wing 

vector 


Ror the fan <= 
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A 
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C 
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P 
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D 



h 
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J 

^3 

n 

q. 

R 

% 

r 

S 


p 

area of duct exit plan* %'D /\ 

p O 

propeller disc area, x(R„-R,Rq) 

p VA) 

propeller chord length. Rig. 9. 
chord length of the duct , / •* . 

G-lauert series coefficients for Eq. ( 3 * 6 ) 

pressure coefficient 

lift coefficient for fan blade section 
lift curve slope for fan blade section 
diameter of duct in the exit plane 
propeller diameter 
fan blade thickness 

incidence angle of fan blade section measured 
from the line of zero lift. Rig. 9. 
advance ratio, 

P 

length of centre body 
fan rotational speed, rev, /sec. 
free stream dynamic pressure 
radius of duct exit plane, D/2 
radius of fan tip 

mean radius of equal area element of fan disc, 
wing area 
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x 

z 


OCj: 


P 


duct thickness 
induced axial velocity 

axial velocity induced bj?- centre body singularity 
distribution 

axial velocity induced by duct thickness 
distribution 

axial velocity induced by the vortex cylinder 

trailing from the duct trailing edge- 

axial velocity induced by duct-bound vortici'ty, Yjj 

total axial velocity induced by all internal 

vortex cylinders trailing from the fan 

free stream velocity 

induced radial velocity 

radial velocity induced by centre body singularity 
distribution 

radial velocity induced by the vortex cylinder 

trailing from the duct trailing edge 

total radial velocity induced by all the internal 

vortex cylinders trailing from the fan 

radial velocity induced by duct-bound vorticity, 

*D 

axial distance from duct leading edge 

number of equal area annuli making up fan disc 

area 

free stream angle of attack with respect to the 

fan axis of rotation, 

fan blade section pitch angle, Fig. 9. 



circulation bound to fan blade, Eq.(3»l) 
convergence criterion 

strength of outer trailing vortex cylinder 
axially symmetric component of duct bound 
■’orticity , i)q. (3-6) 

strength of w-th inner trailing vortex cylinder 

duct bound vorticity component due to angle 
of attack, Eq. (3.3) 

-1 

transformed axial distance, x = (l~c cos 0) 

free stream density 

azimuthal angle 

fan rotational 3 peed, rad/sec. 



CHAPTER I 


IITROiXJCTION 

Of all the known aircraft configurations the fixed 
wing aircraft is the "best logistically , tactically and econo- 
mically. Its utility would he greatly increased if it could 
also take eff and land vertically, hover, and move to the side 
and rear at low speeds. Phis can be achieved if the forward 
flight aerodynamics of the aircraft is not compromised and 
an efficient method of converting the propulsive power to 
vertical lift is designed. Amongst the many VTOL configu- 
rations, the fan in wing system is very promising because it 
enhances the utility of the aircraft with the least compro- 
mise to its present characteristics. This system consists of 
a ducted fan or fans submerged in the wing or the fuselage 
with an arrangement for the distribution of the' propulsive 
power to drive the. ducted fans. The ducts have an arrangement 
of vanes which smoothly close the duct inlet and exit during 
forward flight. 

The flight mode of a VT01 aircraft can be divided into 
three regimes. They are lift-off, transition ahd cruise. 
Transition comprises the complete flight regime lying between 
the static hovering aircraft at one end and the cruising 
aircraft at the other. The performance of the fan-wing system 
in forward flight is not very different from the conventional 
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aircraft. However, the complexity comes in evaluating the 
aerodynamics during lift-off and transition. This report 
concerns itself only with the transition regime. From the 
aerodynamics point of view some of the major , problems faced 
by the fan-wing system during transition are: 

1) Ducted fan in cross-flow, 

2) Interference effects of fan on wing and wing on fan, 

3) Effect of vanes and inlet and outlet shapes of the 
duct on lifting efficiency of the ducted fan. 

4) Interference of free stream with fan jet and its 
influence on the wing undersurface, 

5) The effect of ground proximity on the system. 

A considerable amount of experimental work has been 
carried out on wing/body submerged lift-fans for V/STOL 
applications. Some of the conclusions, which are relevent 
to the present study, are presented. Wardlaw and McEachern £ 3 ] 
and Duvivier and Me Cal bum £ 4 } observed that under static con- 
ditions the combined lift of the fan-wing system is lower 
than the thrust of the fan alone. Hickey and Ellis [ 5 ] 
report considerable lift for the wing fan system over the 
fan thrust. Hackett [6J observed a slight improvement in the 
static lift due to the attachment of the wing to a fan- 
nacelle system. References J 7-12] do not report any loss of 
lift for the wing fan combination under static conditions. 
However, References [5»7,10,1l] report a lift loss for the 
wing-fan system at low flight speeds for a given fan rpm. 
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Spreeman [lO] reports no change in fan-wing total lift with 
forward speed, Trebble [13] mentions that T/yatt [l4,15] 
observed loss of lift at all forward speeds investigated. 

Goldsmith and Hickey [16] observed no loss of lift on full 
scale models. All other experimental work surveyed on the 
fan-wing system indicates that there is a rapid recovery of 
the total lift and a subsequent increase even beyond the static 
lift at zero wing incidence. 

The behaviour of the fan-wing system at incidence, the 
effect of inlet and exit vanes on the total lift, drag, 
pitching moment and fan power, effect of ground proximity have 
all been studied in some of the references given above. However, 
all experimental work suffers from a degree of unreliability 
because of complicated wall interference effects on models with 
a near vertical jet efflux. Almost all the results quoted 
above are without any correction for the tunnel wall interference. 

Theoretical studies suggest, in general, that the loss 
of lift at low forward speeds is due to the fan jet induced 
suction on the underside of the wing. The subsequent recovery 
and increase of lift at higher speeds is due to the fan induced 
upwash in the plane of the wing and the associated additional 
circulation of the wing. Norland [l9] discusses the case of 
a two dimensional wing with a submerged lifting fan. He rep- 
laces the portion of the wing ahead of the hole by an ellipse 
and the portion behind by a flat plate. The gap between the 
ellipse and the flat plate is equal to the fan diameter. He 
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suggests that the lift loss at low flight speeds is due to 
flow seperation on the under side of the wing behind the fan 
and the subsequent recovery and increase in lift is due to 
the reattachment of the flow. He determines the circulation 
of the ellipse and the flat plate by using the condition of 
the jet eff Lax velocity over the fan disc and the Kutta con- 
dition at the leading edge or the trailing edge depending on 
the jet efflux velocity. The main purpose of the paper was 
to investigate the lift- discrepancy-decrease of lift at low 
V/V. and a subsequent increase. 

v 

Martin [20] extends the two dimensional approach of 
Norland to a three dimensional case using the lifting line 
theory. He does not consider the effect of the wing on the 
fan. The wing is approximated by a discrete number of rec- 
tangular .sections each having a constant strength bound vortex 
and a corresponding trailing vortex system. The effect of 
the fan is considered as an additional induced circulation 
in the rectangular section containing the fan. The strength 
of the vortex system is determined by applying the condition 
of flow tangency on the wing along with the proper jet efflux 
velocity on the fan disc. 

McEa chern and Currie [2l] replace the two dimensional 
wing with a submerged lifting fan by a vortex lattice. The 
strength of the vortices in the lattice is determined by 
applying the proper surface boundary conditions at control'.*;, 
points chosen on the wing surface and the fan disc. 
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Monical [_22j uses the same vortex lattice model for 
a finite wing but considers the effect of the jet efflux in 
a different manner. He replaces the jet plume by a curved 
tube of square cross-section made up of a network of vortices, 
their strengths being determined by taking proper boundary 
conditions at a number of control points on the surface of 
the tube. Both McEachern [2lj and Monical [22] replace the 
wing-fan system by a vortex network. This network though 
adequate for representing the wing does not represent the fan 
and the fan jet adequately. 

Murti [24] uses a modified lifting line for the wing 
and tries two different models, a vortex ring and a vortex 
tube, for the fan. He considers the effect of the fan on the 
wing which results in a modified lifting line, but does not 
consider the effect of the wing on the fan. The vortex ring 
model predicted some fan induced wing lift but the trend was 
very different from what was observed experimentally [ 7 ] . 

The straight, semi-infinite vortex tube model did not predict 
any induced lift on the wing. The vortex tube model for 
the deflected jet gave some induced lift but it was much smaller 
than what was observed. The results of the above studies agree 
at best qualitatively with the experimental results. 

In this report, a detailed model for the fan-wing 
system has been proposed and studied. The wing with a sub- 
merged lifting fan has been separated into a wing with a duct 
and a ducted fan. The wing with a duct has been replaced by 
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a network of vortices and a series of singularities have "been 
used for representing the fan. A computer program in Fortran 
IV for use on IBM 7044 has £een developed for the computer 
simulation of the problem. The program uses an existing 
lliSA computer program jjs] for analysing the ducted fan in 
cross-f low. 

The total lift produced by a fan-wing system can be 
separated into its components as the lift produced by the 
fan, the lift produced by the wing, the fan induced wing lift 
and the lift induced due to the ducting of the propeller. The 
theoretical model calculates all the components taking into 
account the interference effects. 

The interference of the wing on the ducted fan can be 
considered to be due to two effects. The first is due to the 
wing induced upwash in the duct. This has been found to be 
negligible as compared to the fan induced inflow at moderate 
and. high values of Cj,. The second is due to the turning of 
the flow by the wing. The flow on the wing surface is tan- 
gential to it and thus perpendicular to the fan axis of rotation 

For the interference effect of the fan on the wing the 
induced upwash on the wing due to the fan singularity system 
can be determined. This upwash in turn determines the wing 
singularity system. However, a simplified model for the 
effect of the fan on the wing has been suggested and tried 
here. It assumes that the effect on the wing due to the other 
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fan singularities is negligible as compared to the duct 
bound vorticity effect. Earlier work [24] which takes into 
account the deflection of the ’fan slipstream due to the 
forward velocity suggests that the effect of the deflection 
is negligible. Hence, in the present model, the curvature 
of the fan slipstream is neglected. The theoretically pre- 
dicted .’fan induced wing lift has been found to be in agreement 
with the experiment [5] . 

In Chapter II, a vortex lattice model for the wing 
with a duct and the wing is discussed. The wing is replaced 
by a flat plate which in turn is replaced by a vortex lattice. 
Elow tangency condition is used to determine the strengths 
of the vortices in the lattice. Eor the wing with a duct, 
the same model is used except that no vortices are assumed to 
lie in the duct. 

In Chapter III, the fan model is discussed. The fan 
has been replaced by a series of singularities representing 
the fan duct, fan disc, fan wake and the fan centre body. 

These singularities are determined by an iterative procedure 
and they are used to get the fan inflow and forces and moments 
on the fan. 

Chapter IV discusses the proposed model for the fan- 
wing system and the various interference effects. 

Chapter V gives the details of the example considered 
for comparing the theory with the experiment. The model is 
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a semi-span tapered wing with a span of 60 inches, root 
chord of 40 inches and a taper ratio of half. It has a single 
wing submerged lifting fan. 

Chapter VI gives the results and the conclusions along 
with suggestions for further work and improvement on the 
suggested model. 

The computer program is given in Appendix C. 



CHAPTER II 


WI KG MOREL 

The wing is assumed to be thin and is replaced by a 
flat plate at an angle of attack. A vortex lattice represen- 
tation is used for representing the flat plate. The lattice 
is made up of a series of horse shoe vortices, as shown in 
Eigure 4. The arrangement is so chosen as to concentrate 
the vortices in regions of steep pressure gradients i.e. near 
the leading edge and near the hole. They are arranged In 
chordwise strips for ease of computation. For the case of 
the wing with a hole, no vortices are assumed to lie in the 
hole, Figure 5* 

Having determined the geometry of the vortex lattice, 
it is necessary to determine the strengths of all the indivi- 
dual vortices. This is done by suitably choosing as many con- 
trol points as the number of vortices and enforcing the flow 
tangency condition at each control point. This provides an 
unknown vortex strength and a known boundary condition for 
each of the horse-shoe vortices in the lattice. 

1. MATHEMATICAL FORMULATION: 

From Biot Savart law the induced velocity. dv due to 
an element ds of a vortex of strength f* at a point distant 
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R from the vortex element, is given by, 

ds x R 


dv 


_n 

4ix 


...( 2 . 1 ) 


R 


Where d_s is the coordinate along the vortex such that looking 
along it the circulation is clockwise. For a vortex of 
length 1^1^, Figure 2, the magnitude of the induced velocity 
is, 

1 V j = (cos e 1 + cos e 2 ) , ...(2.2) 

the direction being given by ds x ia. 


A right handed cartesian coordinate systems with 
origin at the leading edge of the wing root chord, C R , is 
used to represent the wing, The positive x axis lies along 
the wing chord pointing towards the wing trailing edge and the 
y axis is positive to the starboard side, Fig. 1. The horse- 
shoe vortex abed of strength f* induces a down-wash at the 
control point P lying in the plane of the vortex. The down- 
wash 'w* is taken positive in the negative Z direction. Mow, 

w ’ w ab +w bc +w cd’ ...C2-3) 


where the subscript denotes the portion of the horse-shoe 
vortex abed. Let, 


X = y p - ...(2.4) 

x f = y p + y v , 
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with x and y being non-dimensional quantities with respect 
to the root chord C^. From Biot Savart law, 


P 


w ab “ TP ' c”’T/+sT ^ 


X 


rrrrzr + i) 


{ X 2 +(Y+s) 2 


Wi 


be 




• (_Li + 

4u C X 1 

x *- •'P -\T . _ \ if uo , /u 4- 


w 


cd 


-'- Sy “), 

f X^+CY+s) ^ if X^+(Y~s) 

p _ / „x_ 

4tc C^IT-sT V —o — 


+ 1) 


]f X 2 +(Y-s ) 2 


Hence , 


w 


Jd 

4-tc C 


“ f(X,Y,s), 

R 


..(2.5) 


where , 

f (X,Y ,s) 


Y+s 

XfX 2 +( Y+s) 2 



fX 2 +(Y+s) 2 


ll=sl 


XVX 2 +( Y-s) 2 


Ti-iT 


(i + 


X 


fX 2 +( Y-s) 2 


. . . (2.6) 


Similarly the down wash at P due to the image vortex 
is obtained by replacing Y aby Y 1 in Bq. (2.5). Denoting 
f(X, Y',s) by f', the total downwash at P due to vortex abed 
and its image is, 

w = P(X,Y,Y‘,s) ...(2.7) 

R . : 

where , 

B = f + f» 

Ek e funotion-B — dep e n d m only on lir5 ~~ggome try i7e~. — the 
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'The function I' depends only on the geometry i.e. the 
vortex span , location and the location of the control point, 
generalising the above argument for a lattice of 2n vortices 
and n control points, we have. 


{wl - |JF] {-^ J , ...(2.8) 

where, the element F . . of the matrix |F| , represents the 
influence of the j-th vortex on the i-th control point. 

For a number of vortices, at a given spanwise station, the 
sectional lift is, 

i = pv zr , 

so that the section lift coefficient is, 




1 

2 


L 



...( 2 . 9 ) 


Also , 


8-j^ — 2 ( 0^ Ax) , ...(2., 10) 

where, C is the pressure coefficient assumed to be constant 

Jr 

over a small fraction of the chord Ax and the summation is 
over the chord of the assumed spanwise station. It can be 
shown that for each vortex, circulation P can be related 
to a pressure coefficient C at the same chordwise position 

Jr 

as the vortex centre, cf. Appendix A. Thus, 



_L 

Ax 


2 P 


YC 


B. 


.. .( 2 . 11 ) 
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Now from Eq. (2.8), we have, 


8te 

“V 1 


{ w \ = 




or 




. . . ( 2 . 12 ) 


Multiplying Doth sides of the Eq. (2.12) hy [i/ax] and 
using Eq. (2*11) s we obtain, 


^ = rfU ixj - 1 if} * 


rpJ 


... (2.13) 


Also the local moment coefficient can be obtained as, 

f*’ c p j = [x'J {Cp2 ...(z.14) 

t 

where , x^ is the distance of the vortex centre from the 
leading edge of the wing or the distance from the line about 
which the moment coefficient has to be determined. Now to 
obtain the overall lift and moment coefficients, Eqs.(2.l3) 
and (2.14) can be numerically integrated first in the ©taord- 
wise and then in the spanwise directions. Thus, 

b/CR x te 

Ct = / / { C | dxdy ...(2.15) 

-b/CR x le l PJ 

and, 

b/CR x te 

C Cm = f f { X' Cj dxdy ...(2.16) 

-b/C R x le 1 P i 

To get the conventional C^ and C M- C i should be multiplied 
by XC^/S and C^ should be multiplied by C^ASe^),, 
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where , 




While choosing the vortex network for the wing the 
following must be borne in mind. The control points should 
not lie on the bound or trailing portion of any vortex. 

The choice of Ax determines the vortex and control point 
locations. It should be so selected as to concentrate the 
bound vortices in the portions of steep pressure gradients 
e.g. leading edge and near the duct, for selecting ax, 
the wing was divided into spanwise sections whose widths 
were chosen so as to concentrate more vortices near the 
duct. Each spanwise section was then divided into a number 
of rectangular sections whose length was small near the 
leading edge of the wing and large near the trailing edge. 

In each of these rectangular sections a bound vortex was 
placed at a quarter of the length of the section and the 
control point was placed at three quarters, from the leading 
edge of the section. The number of rectangular sections 
determines the number of vortices. The only limitation on 
the number of vortices is the ease of computation and the 
time available. 

A computer program in Eortran 17 has been made for 
the vortex lattice method, cf. Appendix C. The method used 
for numerical integration is Simpson’s method with unequal 
interval, cf. Appendix B. 
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Since the theoretical pressure distribution for a 
flat plate is infinite at the leading edge a curve of the 
form ax ' + bx ^ f was fitted to the first 

three calculated points of the spanwise station, and inte- 
grated from the leading edge to the first point. For the 
remaining points^ starting from the first point to the 
trailing edge? Simpson^ method was used for integration. 



CHAPTER III 


PAN MO PEL 

A detailed singularity distribution is used for repre- 
senting the fan because of the availability of a 17ASA Computer 
Program [18] . A source and a sink is used for representing 
the fan centre body. The propeller disc is replaced by a 
series of equal area annuli each having a constant strength 
bound vortex and a cylindrical shed vortex. The duct is 
replaced by a duct bound vorticity distributed from the 
leading edge to the trailing edge. The effect of the duct 
thickness Is represented by a series of sources and sinks 
along the duct chord, Pig. 8. 

Por a given angle of attack a^, advance ratio J and 
fan and duct configuration an iterative procedure is used 

* 

to arrive at the fan input velocity profile, thrust, normal 
force and pitching moment coefficients of the fan. The 
method used treats the ducted fan at an angle of attack as 
ducted fan with axial flow and ducted fan with crossflow. 

The two solutions are superposed to obtain the final solution. 

1. AXIAL PLOW: 

Por the ducted fan in axial flow the local blade per- 
formance is obtained by using blade element theory. The bound 
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circulation , 

Pw 

RY 


p , in the w-th annulus is, 


w 

_1 p 

2 1 R 


h u r u) rs 

+ 1 ] 


1/2 


V 


...(3.1) 


w 


Where, is the blade sectional lift coefficient, b is 

the blade chord, R is the propeller radius, u is the induced 
axial •velocity, r w is the mean radius of the w-th annulus, 

V is the free stream -velocity and (X) is the angular velocity. 
Strength of the w-th internal trailing vortex cylinder is, 

1/2 




= [0 


z 

2 


li 


}jh 2 , 

Y ~ i- v 1 " r “ y > + , 

m=w+1 ixJ 


( 


fw 

RY 


Jw+1 
RY 




z 

(1 + 2 
m=w+1 


7m \ 
V ' 


ydsere 


J ' = t>ih 


jL 

tub 


. ..(3.2) 

.. .(3.3) 


T m is the strength of the m-th shed vortex cylinder, N is 
the number of fan blades and z is the number of equal area 
annuli making up the fan disc area. The strength, 7^, of 
the vortex cylinder shed from the duct trailing edge is, 


*z 


= D + 




(~) 1 
V RV ; J 


1/2 


- 1 


(3.4) 


The Eqs. (3-l)-(3.4) are functions only of the fan 
blade geometry, advance ratio J and the fan inflow profile. 
The fan blade geometry and advance ratio are known. Thus 
knowing the inflow profile the above circulations can be 
determined. Row the inflow velocity to the fan, is 

given by the axial components of the velocity induced by all 
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the fan singularities and the free stream axial velocity, 


% 

Y 


* D u 
T” + ~ 


U 1L 

'V 




z-1 

£ 

w=m 



...( 3 . 5 ) 


where, u denotes the axial component of the induced velocity 
and the subscripts ¥ , q D and CB denote the duct bound 

vorticity, z , duct thickness and centre body respectively. 


The free stream, trailing vortex cylinders and the 
centre body induce components of flow through the duct cam- 
ber line. To cancel this flow and make the net flow tangential 
to the camber line a bound vorticity, ^ , is placed on the 
duct, it is assumed to be in the form of a Glauert series. 

~T = °o Got "i" + 2 sin ne ...(3.6) 

. n=l 

The C n coefficients are found by the flow tangency condition 
on the duct reference cylinder. The duct reference cylinder 

is defined as that cylinder having the same radius as the duct 
trailing edge. Thus, 


v 


"JD 


4- 7 . + V 




+ V 


CB 


= (V + u + u + u + u rp ) 
*D * Y w CB 


dr 

e 

dx 


W h ere JS. ls the effective slope of the camber line and v 
denotes the radial component of induced velocity. 

2. CROSS-BLOW: 

The ducted fan in a cross flow is replaced by a thin 
cylindrical duct in a cross flow Y sin a f . A non-axisymme- 
tric vorticity distribution, y a , is placed on the duct to 
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cancel the cross flow through the duct. Thus y > the 

a 

duct bound vorticity has the form, 

~Y = sin cos 0 |_C Q cot + 2 C n sin n8 J ...(3.8) 

n=1 

In addition , there is a distribution of free trailing 
vortex filaments because of the variation of bound vortex 
strength around the circumference of the duct. The effect 
of the angle of attack on the fan loading is not included in 
the computer program because of relatively small contribution , J 18] 

3. HERAT I'VE CALCUIAIIOH PROCEDURE: 

'An initial inflow profile of 2V is assumed and the 
blade element calculations are made to obtain the fan and the 
fan- wake characteristics. The flow tangency condition Is 
then applied to determine the duct bound vorticity. Having 
found all the singularity distributions the fan inflow profile 
is calculated. This is compared with the previous inflow 
profile. If the two are not the same, then a new inflow, 
mean of the two, is assumed and the procedure repeated. A 
series of iterations are gone through till the flow profile 
converges to within the desired accuracy. The forces and 
moments on the duct and the fan are then calculated using the 
final values of the axisymmetric and nonaxisyrometric singu- 
larity distributions. . 



CHAPTER IV 


COMBINING PAN- WING MODEL 

The combined fan-wing system is replaced by a vortex 
lattice for the wing portion and a set of singularities for 
the ducted fan portion and its wake. Ihe wing singularity 
system is the same as the singularity system for a wing with a 
duct , cf. Chapter II. The ducted fan singularity system is 
the same as given in Chapter III. 

The total lift produced by the fan-wing system is made- 
up of the lift produced by the wing, the fan induced wing lift 
and the lift of the fan in the presence of the wing. The last 
mentioned lift component consists of the lift of the fan and 
the lift induced due to ducting of the fan. 

The method used for solving the problem, considers each 
of the lift components separately and sums them up to get the 
total lift produced by the fan-wing system. The strength of 
the wing singularities and hence the forces are determined by 
using the flow tangency condition on the surface of the flat 
plate used for representing the wing. For the forces on the 
fan in the presence of the wing, the ducted fan is considered 
together with the interference of the wing on the fan. The 
wing interference effect on the fan is considered due to two 
effects. The first is the wing induced upwash in the duct. 
This has been found to be negligible as compared to the fan 



induced inflow at moderate and high values of C^. The second 
is due to the turning of the flow by the wing. The flow on 
the wing surface is tangential to it and thus perpendicular 
to the fan axis of rotation. Thus for the ducted fan in the 
presence of the wing, it is sufficient to consider the fan 
at approximately ninety degrees incidence with respect to the 
fan axis of rotation. For this incidence and given advance 
ratio the fan program determines the strengths of the various 
singularities of the ducted fan and the resulting inflow, forces 
and moments. 

For the fan induced wing lift, the induced upwash on 
the wing due to the ducted fan singularities is first deter- 
mined and the corresponding wing singularity system found. 

This is the interference effect of the fan on the wing. It is 
considered in more detail in subsection 1. 

Now the strength of the singularity system representing 
the fan-wing system is completely determined. The correspond- 
ing velocities, pressures and forces can be computed. 

1. EFFECT OF THE .DUCTED FIN ON THE WIND: 

The singularities used for representing the ducted 
fan and its wake are: source and sink distribution along the 
duct chord representing the duct thickness, axi symmetric 
duct bound vorticity due to the axial flow component of the 
free stream, non-axi symmetric duct bound vorticity due to the 
cross-flow component of the free stream, trailing vortices 
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due to variation of the duct hound vorticity along the 
circumference ? the propeller hound vorticity representing 
the propeller,- the shed vortex cylinders representing the 
propeller wake, and a point source and a point sink represent- 
ing the propeller centre body. All these singularities 
influence the flow on the wing. 'Their influences are derived 
separately and then super-imposed to get the effect of the 
ducted fan on the wing. 

EFFECT OF DUCT THICKNESS: 


The duct is assumed to he thin and as such the effect 
of its thickness on the wing is neglected as a first approxi- 
mation. 


INDUCED 'VELOCITY DUE TO THE AXI SYMMETRIC DUCT BOUND VORTICITY: 

The induced velocity at a point P(x, y, z) or (x, r, 0) , 
Fig. 13, due to a vortex ring of strength P , is given hy v 

w X 

the axial component and v r the radial component, [2S] . 


v x (x,r) 


1 


2-rcr ' 


fx^+(r+l) 1 


^F(k)- [l + “T^: \ 72 3- 


x 2 +(r-l) 2 


S(k) 


\ * 


... (4.1) 


and, 


v r (x,r) 


r 


2xr' 


rfx 2 +(r+l) 2 ^ 


\ F (k) 


[1 + — pt ^— — p ] 
x 2 +(r-l) 2 


E(k)| 


.. .(4.2) 
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where s r' is radius of the vortex ring, 


P(k) and E(k) an 
first and second 


the complete elliptic 
kind respectively and 


4 r ___ 

- (r+l) 11 


x = x/r' , 
integrals 


r = r/r 1 , 
of the 


...( 4 - 3 ) 


Tiovs t.'f duct hound vorticity is represented as a cylin- 
drical vortex sheet with vorticity varying continuously from 
the leading edge to the trailing edge of the duct, let ^ 

be the circulation per unit length of the duct, from [js] , 

5 

= h 'z ( C o c0_t ~c~ + 2 G r. sin n9), ...(4.4) 

n=1 

where, 9 = cos (•—£—) 

I'hus the strength of the elemental vortex ring, Pig. 14 , 


i ’ = (x ’ ) dx ’ 

Q 5 

= tf (C Cot ~~ + £ C n sin ne) dx J ...(4.5) 

20 2 n=1 

putting, X = x/r’ , X' = x’/r’ , r = r/r’ 


= r ! (C 0 Cot -§- + 2 G n sin n 9 ) < 3 X’ 

n=1 


'z ' J ° 2 


Prom Eq. (4.1) 


dv x (X,r) = - - 


r’Y^ dX’ 
2mr ’ 


f(X-X') 2 +(r+l) 2 


{ P(k') 
L 


[1 


2Lz~1 1 . 


2 ./.. u\2 


(X-X’) 2 +( r- 1 ) 


iOc*)] 


where , 



(X-X') 2 +(r+l) 



* . .(4.6) 




c/£r 1 


r~ 

q 5 

Cot + £ C n sin n9j 

^ n=1 

or 

v (0 , r) = - / 0 

x --c/2r' ,c 

v 2 (Co 


1 

>rx |2 +(r+l) 2 

(F(k')- 

fu — 2 f r -Ul- ] 

- (X')t(r-I) 2 


1 ( k ' ) ) dX’ 

$ 

. . . (4.7) 

and, 

c/2r ! 

v (0 , -- d * Y 

-c, dr 

(c 0 Cot • 

q 5 . 

- r + b C r , sin n0) 
n=1 



t:: 0?(k 5 
) 2 

X' t “+(r-l) d 




ryx 1 +(r+1 


B ( k ' ; ) dX’ 

0 

. . . (4-8) 


It should he noted that v is positive in the positive 

-X. 

x direction , which is the z direction in the coordinate sy sin- 
used for representing the wing. Also note that downwash is 
taken positive in the wing coordinate system. 


IN DUG Ah VB10CI1I hUS 10 THE NON AXISYNiNSlitlC DUCi BOUND VORTIC: 

1'he duct bound vorticity, per unit length of the duct, 
due to the angle of attack of the duct is given by |j8] , 


Y a = Y sin a f cos 0 1 f c o Cot - "f* + 


+ £ C n sin n8 

n=l 


. . .(4.9) 


__ -t O t 

where 0 = cos (--*-) s a _^ is the angle of attack of the fan, 

and 0 and 0* are shown in fig. 13. 
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for a single vortex ring of strength P [26] , the 
induced axial velocity is, 


v x (x 5 r) = 


2tc 

f 


p 


4icr ’ 


, r„cos. . J_. . 


X 2 +r 2 +1-2r cos ( 0—0 7 ) 




d0 1 


. . .(4. io: 


refer Fig. 13* 


Now for an elemental vortex ring, 


1’ = Y r ' dX ' 

c ' a 

and the induced axial velocity v^ is, 

C/Zr ’ V 

0, r) = / J " ""d 

^ -c/2r ! oo 


.X..,c2,s.iAvll)™l d0 'dX 

((X-X 1 ) 2 + r 2 +1-2r cos ( 0 - 0 ’)) 5/ 2 

... (4. 11) 


refer Fis. 14 


Or the influence of the duct bound vorticity due to 


angle of attack, on the wing is 


\x (0 > r) 


Y sin a 


c/ 2r ' 


e 


4* i/ 2r . (C ° ^ " 


+ 2 Cjj sin n9 

n=1 


2 % 

J cob 0 ' 
0 


( (X-X ' ) 2 -fr 2 +1-2r cos (0-0')) 


d0 . 


.. .( 4 . 12 ) 
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Similarly, the radial velocity component due to a 
single vortex of strength f , [26], refer fig. 13, 


2 % 

= / 
o 


r 


X cos 0 


4it:r ’ 


(X 2 +r 2 +1-2r cos (0-0’)) 


~/2 


d0 * 


( 4 -. 13 ) 


i'hus the induced radial velocity component of the duct hound 
vorticity due to angle of attack, v^, on the wing is, 


v or(0’ 0 ’ r) 


X...sia-£s c y 2r ' T <( 

4tc ‘ J - U o 


■ c/ 2r ! 


iHO. cot 


2 % 


+ 2 Bill 110) f COB 0 1 (X 5 

n=tl o 

+ r 2 + 1-2r cos (0-0'))" 5 ' /2 d0* dX* 

. ..(4.14) 


Ihe integrals involved were integrated numerically to 
get the influence of the duct on the wing. 

INDUCED VELOCITIES DUE TO THE TRAILING- VORTICES AND PROPELLER 
BOUND VORTICITY: 

As a first approximation the effect of the trailing 
vortex filaments, due to the variation of duct hound vorticity 
around the duct circumference, and the propeller bound vorti- 
city on the wing has heen neglected. 

INDUCED VELOCITIES DUE TO THE PROPELLER SHED VORTEX CYLINDERS: 

The shed vortex cylinders are assumed to he straight and 
extending to infinity. As such their effect is confined to the 
area within the cylinder and have no influence on the wing 
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flow field. However at higher forward velocities the shed 
vortex cylinders tend to bend downstream and they induce a 
flow field on the lower surface of the wing. Murti j 24- 1 consi- 
ders the effect of the deflected fan jet and finds its influence 
to be very small. Hence in the present model the curvature 
of the fan jet is not considered. 


INDUCE I) VELOCITIES DUE 10 THE CENTRE BODY: 


The model assumed for the centre body, in the fan 
program, is a Ran'kine body. G-iven the dimensions of a body of 
revolution ( ^ and h) , Eig. 11 , it can be approximated by a 
source and a sink whose locations are given by, 


C? 0 2 - a 2 ) 2 


h 2 fa 2 + h 2 


Induced velocity at a point ( § , r) is given hy : 


V 


(T 2 „2v2 

“ a / 5 


_ (f + a} 

4 5 0 a ' l [(5«) 2 +r 2 ] 3/2 


(1 -a) __ 1 

LCS-a) 2 ^ 3/5 J 


v(S, r ) 

V 


<Sc 


a 2 ) 


4J 0 a 


J. 




r 

o 3/2 


... (4-. 15) 


...(4.16) 

. . . (4. 16) 


( 5 -a ) 2 + r 2 p75 J 
where u and v are the axial and radial components. However, 
for the example considered there is no centre body and hence 
the induced velocity due to the centre body is not taken into 
a cco~unt. 
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From the above discussion, it will be seen that in- 
effect only the duct bound vorticity is considered for getting 
the effect of the fan on the wing. The effect of other singu- 
larities is assumed to be negligible as compared to the duct 
bound vorticity. The deflection of the fan slip stream with 
forward velocity is also neglected and so the fan slip stream 
is semi infinite and nearly perpendicular to the free stream 
for the assumed fan-wing model. 



CHAPTER V 


EXAMPLE 

As ci n example we consider a ssmi-span trapizoidal wing 
model, Fig. 3, with a submerged fan, which has been treated 
experimentally elsewhere [ 5 ] . The details of the model are 
given below. 

WING-: 

The wing is a semi-span trapizoidal wins of aspect 

ratio four, semi-span area of 12.5 square feet, taper ratio 

of 0.5 and a 0° sweep of the half-chord line. The airfoil 

section of the wing is MCA 16-015 (employed to give nearly 

equal fore and aft duct length). The aerodynamic axis of the 

fan is on the mean aerodynamic chord at 0.5 c 

A 

SAN : 

The fan has four blades. The blade form characteris- 
tics are given in Pig. 7. The airfoil section of the propeller 
bladexxis NACA 16-309 on a 66 series mean line. The 
clearance between the propeller tip and the duct is 0.006 
times propeller radius. 

for the purpose of mathematical simulation, the semi- 
span wing is replaced by a vortex lattice of 64 horse-shoe 
vortices as shown in Pig. 4 . For simulating the wing with a 
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vortex lattice of 48 horse-shoe vortices as shown in Fig. 5 • 
Figs. 4 and 5 also show the controls points location on the 
wing, fable 2 gives the details of the vortices and the 
control points for the wing and 'I'able 3 gives the details 
for the wing with a duct, cf. Appendix C. Fig. 8 gives a 
schematic view of the singularities used for representing 
the fan. fhe details are presented in input data, cf. Appendix 
C and ‘Table 1 . 



CHAPTER VI 


RESULTS AliD CONCLUSIONS 

The experimental results given by Hickey and Ellis [5] 
have been used here for comparison with the theory. It should 
be noted that wind tunnel corrections have not been incorpora- 
ted in the experimental results given in [5] . Some correc- 
tions can be expected because at high propeller loadings as 
much as ten percent of the wind tunnel airflow passes through 
the propeller duct. There is a blockage effect due to this 
wake. , the propeller drive motor and the supporting struts 
which are also exposed to the wind tunnel flow. Wind tunnel 
boundary corrections were also not considered. 

The vortex lattice method used for the simulation of 
the wing predicted higher values for the wing lift coefficient 
and the lift curve slope as compared to the experiment, 

Pig. 17 * Prom the experimental curves, it is seen that a 
wing with a duct gives some lift even at zero incidence 
unlike the wing without the duct. As the theoretical model 
replaces the wing with a flat plate no such effect can be 
expected. The theory also does not take into account any duct 
induced separation on the wing upper surface. The wing with 
a duct corresponds to the case of the wing with a fan, fan 
not working or working with very low values of C^. This 
flight region of the fan-wing system is not of interest. Our 
interest lies in the transition regime. 
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The predictions of the fan program for the propeller 
thrust agree very well with the experiment at low rpm s hut 
at higher rpm it predicts higher values for the propeller 
thrust, fig. 18. It is not justifiable to comment on the 
theoretical predictions at high rpm as then the experimental 
results are amenable to large wind tunnel interference and 
blockage corrections. The variation of propeller thrust with 
for a given rpm is small for both theory and experiment. 
However, for the range of q^ considered the variations in the 
theoretical values are negligible. 

for a given experimental value of Cj, and q^, J can be 

found, for this J the theoretical Cj> is found to be different 

from the experimental C-p. This is because the theory predicts 

higher thrust than the experiment at the same rpm. Thus for 

comparison of the results it is a choice between J and C^. 

The theoretical and experimental comparisons are made at the 

same value ’of Cj,. This essentially means that the theoretical 

and experimental propeller thrusts are the same. Thus in this 

comparison, the fact that the propeller gives higher thrust 

\ 

at high rpm has been circumvented, fig. 20 compares the 
theoretical and experimental incremental forces due to the 
fan in the wing. Jit zero incidence the lift produced by the 
wing is made up of the fan induced wing lift and the lift due 
to the ducting of the propeller, from the figure we see that 
except at very low and very high values of C^, the predictions 
of the theoretical model are very encouraging. This high-lights 
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"the fact that for the fan— wing system the duct has a signi- 
ficant effect on the fan induced lift. 


figure 19 compares the total! lift produced by the 
fan-wing system at nearly same value s of C^. The model 
predicts comparable results except at the higher value of C-^. 


The experimental resuLts do not give the fan induced 
upwash distribution and the pressurre distribution over the 
wing. fig. 21 gives the ch or d wise fan induced upwash dis- 
tribution on the wing at three spaxirise stations, fig. 22 
gives the chordwise overall pressurre distribution at the 
same three spanwise stations on the wing. As expected the 
effect of the fan has been found to be considerable near the 
fan but small away from it. Tig- 23 gives the variation of 

C M le with C L* 

It Is suggested that i_n future work, the effects of the 
fan singularities, which have not ~be en considered, should be 
taken into account for obtaining tie effect of the fan on the 
wing, further detailed in flight experimental investigations 
should be carried out on a model ozr "the prototype and the 
theoretical predictions compared a-g&inst the results obtained. 
Experiments devised with the minimum of wind tunnel inter- 
ference or determinable tunnel interference will also be very 
useful for comparison. 
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APPENDIX A 


SLAP ION BETWEEN VORTEX STRENGTH AND COEFFICIENT OE PRESSURE: 


It can be shown that for each vortex, the circulation 
can be related to the pressure coefficient at the same 
chordwise position as the vortex centre. . Take a vanishingly 
thin element of area of length l>x surrounding a typical 
vortex, of. Pig. 12. Constant over A.x is equivalent to 
the assumptions that the induced velocity u over 2A x can be 
replaced by a mean velocity u where u' = f udx. If w u and 
w d are the downwash at the upstream and downstream ends, we 
have , 



f (u 


dx 

Cr 


dz \ 

w cl > 


2u _ + (w d - w u ) 




2u 


^x 

C R* 


. . . (A— 1 ) 


as u and w are of the same order and Az can be made negligi- 
bly small with respect to A.x. Applying Bernoulli f s Eq. 
to the lower and upper surfaces gives, 


r?v 2 + p o = -r ?(T+ 5)2 + Pu = -L § (y-u ) 2 + p x 


= i- ? (4V© 


or 


Pl - Pu 


. . . (A-2) 
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Defining 


P u p o 


'P 


u 


\ ?v 2 


C-r 


G P “ C P1 


C p , 
F u 


' p l 

we obtain for C 


p l - p c 


p : 


C p = 4 (u/V) 


Substituting for u in terms of P we obtain. 


C 7 


_J_ ( 2_C ) 
Ax v vc R ' 
0r“ 


and 


. ..(A-3) 


(A— 4) 


Thus subject to the .,bove assum.ptions every term in Eq.(2.9) 
can be equated to corresponding term in Eq. (2.10). 



APPENDIX B 


SIMPSON'S METHOD WITH UNEQUAL INTERVALS: 

The formula of Simpson's method with unequal intervals for ' 

x 2 

numerical evaluation of an integral, f f(x) dx, is derived 
below. 0 

Take a two dimensional rectangular cartesian coordinate 
system. Let x Q , x^ and x^ be three arbitrary points, on the 
x-axis such that x^ > x^ > x Q . Let the values of the function 
y = f(x) at these points be y Q , y^ and y^. Shift the origin 
to x^ such that, 

x 1 =0, 

x o = -a 1 H, ... (B-1) 

X ^ ■— CCr, H . 

Let the function y = f(x) be approximated by a quadratic, 

p 

y = Ax + Bx + C 


in the interval x„ x 0 . Therefore, 

o 2 J 

y Q = A( a 1 H) 2 - B( a 1 H) + C, 

y 1 = c, 

yg = A(«2 H)^ + B(«2 H) + 0, 


whe re , 

A = 

B = 

C = 



. . . (B-2) 


(B-5) 
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This is the desired formula. for using the above method the 
function value must be given at an odd number of points greater 
than three. 



APPENDIX C 


PROGRAM DESCRIPTION : 

The program has been written in Fortran IT for IBM 7044 
computer. It has been divided into two parts and if it is run 
as indicated, it does not require any other tapes except the 
normal system tapes. The reason for dividing the program into 
two separate parts is computer memory considerations and conve- 
nience of using the system as it exists in the Computer Centre 
at IIT Kanpur. However, it can be run as a single program if 
a scratch tape is provided. 

The computation proceeds as follows. For a given fan 
in wing combination Part I of the program calculates the values 
of the various fan singularities, the thrust developed and the 
fan inflow velocities. Part II uses some of the results 
obtained in the first part to calculate the wing singularities 
in the presence of the fan and gives the overall pressure and 
forces on the wing. Part II can also be used to calculate the 
wing singularities both for a wing with or without a duct in 
the absence of the fan. 

PART I DESCRIPTION: 

Part I of the program is the NASA program, Ref. £18} 
with only slight modifications in the main program so that 
the strength GAM of the vortex cylinder trailing from the duct 
trailing edge and the coefficients C n of the Glauert series 
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representing the duct bound vorticity, are also printed out. 
these are used as part of the input to Part II. 

This consists of a main program and twenty subroutines. 

The relationship between the main program and the subroutines 
is given in I'ig. 15. The subroutines used and their functions 
are given below. 

MAIN controls the flow of program, calculates the local 
blade element performance, calculates the total fan- in-flow, 
compares the initial and computed inflows for convergence and 
calculates the new inflow profile if convergence is not achieved. 
It is in two parts - the initialisation part and the iteration 
and convergence controlling part. 

INPUT determines the quantities used in the calculation 
that are functions only of the duct configuration, centre body 
geometry fan blade characteristics, number of fan annuli. This 
includes fourier series coefficients for velocities and vorti- 
cities which are functions only of duct chord to diameter 
ratio or duct thickness. It also reads the input data, checks 
it for certain specific errors and prints it out. 

PKL does a table look up of the P k ]_ coefficient for 
various values of the duct chord to diameter ratio. The influence 
.coefficients, P^j_ are used in computing the radial velocity 
induced by a distribution of vortex rings. These are functions 
only of the duct chord to diameter ratio. 
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EL1IPS does a table look up for the complete elliptic 
integrals of the first and second kind given the arguments of 
the elliptic integrals. 

-UaivlBfA does a table look up of Heuman Lambda 'function 
given the two arguments of the function.-. 

HUB calculates the location of the point source and the 
point sink used for representing the centre body model. It also 
computes the nondimen sional axial inflow velocity induced at 
the mean radii of the fan annuli by the source sink distribu- 
tion. 

CAMBER calculates the induced camber coefficients combined 
with the geometric camber of the duct to give the duct effective 
camber . 

PROP divides the fan annulus into the required number 
of annuli, computes the inner, outer and mean radii of the 
annuli and interpolates in the blade characteristics from the 
input table to get their values at the mean annuli radii. 

CLA1P does a table look up of the local blade lift 
coefficients taking into account the fact that after blade stall 
the blade lift coefficient is assumed to be C-^ max * 

ARC SIN computes the principal value of the angle given 
the value of the sine of the angle. 
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jsiAi'RUC performs matrix inversion and prints error 
messages if the matrix is singular. 

5‘OURCS does a Courier analysis of a function given a 
table of the function values. It computes an n term (n < 50) 
cosine series to fit the given function. However , only the 
first six terms of the series are used for the calculations 
in the program. 

SRGfUJGr computes the non-dimensional axial inflow velocity, 
u/V, induced at the mean radii of the fan annuli by the source 
rings representing the duct thickness distribution. 

IfXRl'TG- computes the non-dimensional induced axial velocity 
to the fan due to the axi- symmetric duct bound vorticity Y p- 

AIFRUG- calculates the non-dimensional induced axial 
velocity to the fan due to the non-axi symmetric duct bound 

vorticity Y . 

J ° a 

GAiiCYL gives the non-dimensional axial inflow velocity 
to the fan induced by the vortex cylinder trailing from the 
duct trailing edge. 

BICOEF calculates the 1'ourier series coefficients for 
the radial velocity induced along the duct reference cylinder 
by the vortex cylinder trailing from the duct trailing edge. 

AI'ICOEf calculates the Fourier series coefficients for 
the radial and axial velocities induced along the duct reference 
cylinder due to shed vortex cylinders and the vortex cylinder 
trailing from the duct trailing edge. 
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CNCOEE applies the boundary condition of no flow through 
the duct camber line to get the GlaUert series coefficients 
representing the duct bound vorticity. 

PRESS computes the velocity induced on the duct 
reference cylinder by all the singularity distributions and then 
calculates the duct surface pressure coefficients. 

OUTPUT computes the force and moment coefficients and 
prints out all the output of the program. 

1. INPUT RESCRIPT ION; 

This section gives the description of the form in which 
the input to Part I is to be given. It gives the card number, 
the variable name in order of punching and the format in which 
the various duct and fan parameters are to be punched. 

CARR NO. 1 contains any alphabetic and numeric information that 
is desired for identification purposes. This information is 
printed on the top of each page of output. 

CARR HO. 2 contains the duct geometry information. The variables 
to be punched are duct chord-to-diameter ratio (C/D) s the 
axial location of the fan within the duct (Xp./C) , the maximum 
thickness to chord ratio of the airfoil section of the duct 
(t/C), the ratio of the duct trailing edge radius to the fan 
radius (R/Rp) , the ratio of the centre body radius at the fan 
station to the fan tip radius (^03/Rp) and the convergence 
criterion (£). The format for all the variables in this card 
is 1 10.4. 
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CARD NO. 3 contains the four Courier coefficients (R^) of a 
cosine series fit to the duct camberline. If the duct is 
un cambered the four values are all zero. The format for all 
the four coefficients is S' 10. 4 . 


oaRD NO. 4 contains the centre body geometry information. The 
variables to be punched are the centre body length as a frac- 


tion of the duct chord (Iq^/C), the location of the centre body 
nose in the duct coordinate system (X^/C). If the centre body 
extends forward of the duct leading edge then X^-g/C Is nega- 
tive. Ratio of the maximum radius of the centre body to its 

length (r Y /l np ) and the x location of r (X„ ) in the 

max or max r max 

duct coordinate system are punched next. All the variables 
are punched in S' 10 . 4 format. 


CARJ NO. 5 contains six quantities. Number of fan blades (NBLT) , 
number of equal area annuli (NZ) , into which the fan annulus 
is divided, the number of input stations in the fan blade cha- 
racteristics (NZP) , the number of x stations in which the duct 
surface pressure distribution is to be calculated (IR), an index 
(NPRES) which controls the form of the surface pressure coeffi- 
cient and an index (NPRINT) which controls the quantity of 
output. All these are punched in the format I 5. 

CARD NO. 6 contains the fan radii (r/Rp) at which the blade 
characteristics are to be input. There are NZP entries in 
this table and the format in which the values are punched 
is F 10.4. If NZP > 8 then the remaining values are punched 
on extra cards. The values must be put in the order of increas- 
ing radius ratio.. 
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CaHD NO. 7 contains the values of the fan blade chord (b/R p ) 
at the radial stations corresponding to the values in Card No. 6 . 
The format for all the values is P 10.4. 

CARD NO. 8 contains the fan blade pitch angle ( (3) in degrees 
at the radial stations corresponding to the values given in 
Card No. 6. The format for: all the values is P 10*4. 

CARD NO. 9 contains the fan blade thickness to chord ratio at 
the above mentioned radial stations, the format being P 10. 4 . 

CARD NO. 10 contains the values of x/C at which the duct- 
surface pressure coefficients are to be calculated. This card 
contains IR values. The values should be put in the ascending 
order with the limitation 0 < x/C £ 1.0. The format is P 10. 4 . 

CARD NO. 11 is the last card required for a particular confi- 
guration. The first value is the run identification no. (URDU), 
next is the index (UP HI) which is the number of azimuth station 
around the duct at which the surface pressure coefficients are 
to be calculated. Advance ratio (j) is given next. It is follow- 
ed by the angle of attack ( a^) in degrees and the remaining 
five fields contain the values of the azimuth angles 0 in degrees. 
The formats for the variables are I 5,1 5, P 10.4, P 10.4 and 
5 P 10.4 respectively. 

If runs are to be stacked then the following procedure 
should be used. If the same duct fan configuration is to be 
investigated for different J, a and 0 then a new card with the 
same format as Card No. 11 is punched and placed after Card 
No. 11. This can be repeated for as many runs as desired. 
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If a new configuration is to be investigated then a 
blank card should follow Card No. 11 and if should be followed 
by the deck of the new case consisting of Card Nos. 1 through 11. 
The input to Part I for the fan of Reference [ 5 ] is given in 
Table 1 . 

2. LIMITATIONS ON THE INPUT DATA TO PART Is 

In Card No. 2 the admissible range of C/D is 0 < C/D <2.5 
The upper limit is an approximate upper limit because the 
coefficients built into the program have to be extrapolated for 
C/D greater than 2. As P^ do not vary linearly with C/D an 
arbitrary limit of C/D = 2.5 is fixed beyond which extrapola- 
tion is meaningless. The admissible range for the duct t/C 
is 0.06 < t/C < 0.24. This is because of the correction factors 
used in the calculation of duct pressure distribution. If t/C 
is outside this range then the pressure distribution will be 
slightly in error but all other calculations will be correct. 

The suggested range for 6 is 0.001 < 6 < 0.01. The iteration 
will stop when the profiles agree to withiri E x 100 percent. 

The limitations on the variables gi ven i n Oard No. 5 
are that NB1D should be greater than zero- The range for NS 
is 2 < NZ < 24. The range for NZP is 2 £ N2P 1 24. mus t 
lie in the range 1 < IR < 25. If the duct surface pressure 
is not required then also JR should lie in the indicated range. 

If NPRES = 1 , then the surface pressure coefficient is based 
on the free stream dynamic pressure, if ITPRPS = 2 then it is 
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based on the propeller tip speed. For normal runs IMPRINT = 0. 

If IMPRINT = 1 , an extra page of output is developed which 
contains a table of the components of the inflow to the 'pro- 
peller. If IMPRINT = 10 then a page listing the Fourier cosine 
series coefficients for the u and v velocities induced at the 
duct reference cylinder is printed. If IMPRINT = 11 then both 
the above optional outputs are printed. 

The first value of the r/R punched in card No. 6 should 

P 

be that corresponding to the hub and the last value should be 
approximately 1. 

As there is a discontinuity in the pressure across the 
fan thus it is desirable to have the x/C values, punched in 
Card No. 10, to be those corresponding to the points ahead and 
behind the fan. However, no numerical problem occurs if the 
fan station itself is given as the value of x/C. This is be- 
cause the program uses the total pressure upstream of the 
fan to compute the duct inner surface pressure. 

In the Card No. 11 NRUN should lie between 0 and 999, 
NPHI should lie between 0 and 5- If NPHI = 0 then no pressure 
calculations are performed. The advance ratio . ".J must be 
greater than zero. For a the only requirement is that 
cos a / 0. 0 values should lie between 0 <_ 0 <_ 180° because 


of symmetry. 
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3. PART II DESCRIPTION: 

Part II of the program determines the pressure and 
forces on tne wing for three cases. They are for the complete 
wing, the wing with a duct in the absence of the fan and the 
wing in the presence of the fan. For thelast case the program 
uses the Part I results for the values of the C n coefficients 
of the ulauert series representing the duct bound vorticity, 
the strength of the vortex cylinder trailing from the duct trail- 
ing edge and the fan angle of attack measured from the fan axis 
of rotation. The choice as to which of the three cases the 
program should analyse is governed by the values of the two 
indices NNN and NEXT built into the main program. 

Part II is made-up of a main program and six subroutines. 
The relationship between the main program and the subroutines 
is given in Pig. 16. The subroutines used and their functions 
are given below. 

IAIN program controls the flow of the program and reads 
and prints the input data. It nondimensionalise s the vortex 
lattice data with respect to the root chord. The flow of the 
program is governed by the indeces Iffl and NEXT. If NM is 
positive then -only velocities in the duct due to the wing are 
calculated. If NM is zero then only the forces on the wing 
are calculated. If NNN is negative then both the forces on 
the wing and the upwash in the duct due to the wing are computed. 
If the index NEXT is zero or negative then the effect of the 
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fan on the wing is not calculated. If NEXT is positive then 
the pressure and forces on the wing in the presence of the 
fan are calculated. 

WlN'U computes the pressure and forces on the various 
points on the wing and calls other subroutines to get the over 
all lift and moment coefficients. This subroutine also prints 
out the output. 

js&aTINV performs matrix inversion and solves simultaneous 
equations. 

VAR SI M uses the Simpson's method with unequal intervals 
for numerical integration. Some statements have been built into 
the program to tacl-cle the specific example considered. 


CURFIT does . some curve fitting to avoid the numerical 
problems due to theoretically predicted infinite pressure at 
the flat plate leading edge. 

EFFECT computes the effect of the fan on the wing. This 
uses the fan data and some of the results obtained in Part I. 


EL LIPS does a table look up for the complete elliptic 
integrals of the first and second kind given the arguments of 
the elliptic integrals. This subroutine is the same as the 


subroutine ELLIPS of Part I. 




V KANPUR 


f 


A 


233821 
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4. IHPUT DESCRIPTION: 

This sub-section describes the form in which the input 
to Part II is to be given. It gives the card number, the vari- 
able name in order of punching and the format in which the 
various wing parameters are to be input. 

CaRD NO. 1 contains the number of vortices N in the wing 
vortex lattice, the number of control points HIT, on the wing 
and the number of spanwise stations M. The format for all the 
three variables is I 3. 

CARD No. 2 contains more detailed information regarding the 
vortex lattice. There are as many cards in this as the number 
of vortices. Each card contains the x location of the vortex 
(XV(l)), the distance between the successive control point, 
(DELTA X (I)) , the x location of the control point (X P(I)), 
the y location of the vortex centre (YY (I)), the span of the 
vortex (S(I)) and the y location of the control point (YP(l)). 
The format for all the variables is E 10.6. 

CARD NO. 3 is present only when velocity at various control 
points' in the duct in the absence of the fan, is to be deter- 
mined. It contains the x location of the control point in the 
duct (X P(I)) followed by the y location (Y P(l)). The format 
is P 10.6 for both the variables. There are as many cards as 
the number of control points In. the duct at which the velocity 
is to be determined. 
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CARD NO. 4 gives the number of vortices plus one (N 1 ( I ) ) at 
each spanwise station. This card contains M values. The format 
for all the values is I 5. 

CaRD NO. 5 gives the identification number L 1 (I) of the 
first vortex in each spanwise station. There are M values in 
this card and the format in which they are to be punched is I 5 

CARD NO. 6 is present only when the effect of the fan on the 
wing is to be determined. It contains the six C^ coefficients 
determined in Part I of the program. The format in which each 
of them is to be punched is P 10 . 5 . 

The data for the wing of Reference j 5 1 , without the 
duct, is given in Table 2. Table 3 gives the input data for 
the above wing with the duct. It should be noted that values 
of the vortex and control point locations are one fifth of the 
actual values. This was because they were measured from a 
scaled down diagram. The program converts the value to the 
actual values and then non-dimensionalises them with respect 
to the root chord. 

In addition to the above input some information has to 
be input in the body of the program. In the main program, 
the following have to be given, the value of the wing angle 
of attack (Al) in degrees, the value of the fan angle of attack 
(ALP) in radians and the value GAM of the vortex cylinder trail 
ing from the duct trailing edge. The latter two values are 
to be given only when the effect of the fan on the wing is 
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to be determined. In subroutine EFFECT the following data . 
has to be given. The x location of the fan (XF) , the y loca- 
tion of the fan (YF) , both non-dimensionalised with respect to 
the rest chord, the root chord (CROOT) , the C/D of the duct (CD) 
and the duct trailing edge radius (RTS) . 

The only limitations on the wing input data are the 
memory limitations of the computer and size of the matrix 
MAI IN V can invert. 

5. PROGRAM LISTING: 


The following pages give the program listing. The table 
below will act as a table of contents for the progr am listing. 


Program 
Part I 


Part II 



Program 

Page 

MAIN 

55 

INPUT 

59 

PKL 

65 

ELLIPS 

66 

LAMBDA 

73 

HUB 

81 

CAMBER 

83 

PROP 

64 

CLALF 

85 

ABC SIN 

87 

matrix 

63 

POURCS 

89 

SRCRNG 

91 

VTXRNG 

93 

ALFRl'l G 

97 

GAMCIIi 

101 

BNCOEF 

104 

A1MC0EF 

105 

CNCOEF 

108 

PRESS 

110 

OUTPUT 

114 

MAIN 

120 

WING 

122 

MATINV 

125 

VAR SIIvI 

127 

CURE IT 

126 

EFFECT 

129 

ELLIPS 

131 
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SUBROUTINE MAIN(NNNI) 
er-FTC v i A I N 

c ‘ ' JUCTED PROPELLER ANALYSIS PROGRAM 

D I MENS I ON B ( 6 ) 5 BS ( 6 } » 3 A ( 6 ) * S AS ( 6 ) 9 P ( 6 5 6 ) 
u I M E !\i o 1 0 M 0 ( 5 ) 9 h ( 2 5 51 6 ) s a A ( 2 5 9 6 ^ 9 i-> ( 5 ) » D 3 ( 6 I * S C ( 5 ) •> G S ( 5 ) 9 H ( 5 I 
DIMENSION rB( 25 ) »&R < 25 ) »STA( 25 } >TCdlD( 25 ) »rA( 25 ) »XprE J ( 2 5 ) 9 pH I (5) 
[ U I MENTION UG ( 25 ) » U QB ( 25 > s UGD ( 25 > *uv ( 25 ) sUCB ( 25 } 9 vcb ( 2 > ) 

2 s J G p ( 2 5 5 2 5 ) 9 U G A ( 2 5 ) » G V ( 2 5 ) 9 u R V ( 2 5 ) »Cpp( 5»2 5 ) 9 C P ,M ( 5 9 2 ? } 

DIMENSION ALPHA ( 25 ) >STaLL< 25 ) jJSTl( 25 ) »TALk( 20 ) 

Common vZZZ»CD»R‘U» ^ 1 »r 2 »R 3 »PI »b»BS»SA»SAS»p 

Cq v.mqn/ NEAP 1 / NR Um » NBLU » NZ a mZ 9 mPrE 5 9 Ip » n T I ME 5 nErr 9 MP AG »MpH I 9 i\jPR InT 
C 0 M M CN/NEAR 2 / C 9 A. 9 A S » D 9 Di> 9 SC 9 GS »H 

Oo?’’«ON/ NEA p 3 / oppsXp 9 Z*BLO»pb 9 Bp 9 bTA 9 TCBLU»TC »pCBpp » Ap A » A[_F »Xpp ES 9 
1 _rA 9XC6 »Xp 9 ElCBC 9 rmAX , pHI 9C0RJ 9 C0RC6 

OoMMON/NEAr a/UG 9 Uql) 9 UGD 9 U V sUCB 9 VCb 9 UGp $ UGA 9 G V 9 GAm >GrV »Cpp »CpM 
LO.'sMQ.N/ ;mEAr 5 / AR J 9 Ap j p 9 EPS 9 R Ai) 9 C L » ALPHA 9 S I A [_ [_ 9 jbT L 3 T A L \ 

3; j F 0 RMAT( 2 I 5 » 7 Fl 0 «, 6 ) 

9 , ' FCkMAT ( / / / 1 OX 5 42 HPRQGR AM DIu jmqT ConVErGE ON INFlO^ PROFI[_E///) 

9'-: Format (/// 1 ox 9 -x-errop - hxcecutiom terminated*// / ) 

< initialization of subroutines 

c 

it 

T DUM= 0 o a 

MZZZ =0 

CALL E L L I P S ( D U M 9 DUM 9 DU M ) 

call la m b d a (dum*du m 9 d um ) 

CALL PKL ( DUm 9 p ) 

CALL CLALF(O) 

f '-’ZLZ-l 

? I = 3 0 1415 R 26 

RAD= 180 o /pi 

c 

■2 : - NER.P =0 

CALL INPUT 

IF l NERR) 2 b 9 26 s 999 

h CSalF=C0S ( ALF/RAD ) 

ARJV=ARJ 

| ARJVP =AR JP 

I A p J = A R J -K- C S A L F 

s ARJP«AR-JP*CSALF 

I DO 30 K = 1 9 N 7 

I 3 -' UV ( « ) = 2 „ 0 

I NT I ME - 1 

S CORCB =2 o A 

P 1 DO 3 2 K = 1 9 n Z 

I = RB ( K ) / A R J P / U .V ( ’< ) 

ABPI=ATAN(BP I ) 

ALPHA ( K > = ( A.B.P I-BTA ( K ) ) *.RAD 
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api-BPi*BPi 

5dI=SQRT (B pI+loO) 


J = K 

CALL CLALF(J) 

jp ( mERR** - ! ) 25 » 999 ? 999 

C 0 9 T I N U t 

GR V ( < ) =0 o 5*CL*Br { K ) *UV ( K ) * 8p I 
CONTINUE 

GV(NZ)=CR'V (NZ ) *BLD/PI / ARJP+1 c „ 
AGiV=l o 0 

iF(GVUiZ) ) 133 » 233 »233 
GV(NZ)=-GV(mZ) 

AGV=-1 o 

GV (\iZ)=AGV*SqpT (GV(nD ) — 1 ,, 
AGV=lo 

SGV=1o0 
00350*1 ^ M Z 
K=NZ-J 
L = '<+1 

SgV=SGV+GV ( L ) 

SGV2*SGV*5GV. 

GRVD=GRV ( K ) -Gr V ( L ) 
GrVD=GR'/D-x-3lD/P I /AN JP 
DUM»SGV2+GRVD 
AGValoO 

I F ( DUf'io LT o 0 c 0 5 AGV=— l a 0 

GV(k)*AGV*SORT(AGV*DUM)-SGV 

CONTINUE 

0038 N»l *6 

Si (N )=0, 

S A S ( N ) = 0 c 
DO 40 N*l »6 
D03?M=1 »MZ 

SA ( N ) =SA { N ) +A { M 9 N ) *GV ( M ) 

SaS ( k; ) =SaS { k; ) + AS ( jv| 9 pi ) * G ( M ) 
CONTINUE 


DQ41 M = 1 9 6 

SA ( N ) = SA ( N ) +D ( N ) *CORCB 
S AS (■ M } =SAS { N 5 +DS ( N ) *C0RC3 
SAN=GV ( NZ ) 

■"ALL CNCOEF(GAM»C) 

141 N*1 ,6 

, ( ,v. ) =5A { N 5 — D ( M ) *CORCB 
0 o ( m ) =S AS ( N ) -DS ( n ) *C0RCB 

scnpute corrections for low advance ratios 

0U:M=XPRES( 1 ) 

DUmm=ra { 1 ) 

XPRES ( 1 } =xp 
F A ( 1 ) = 1 0 0 

CA U GA^CYl I CD 9 Xp 9 1 9 rB aUG s 1 » XdrEp 9 UG 0 9 PRO 9 pA ) 
IRC®— 1 

CALL VTXppiG ( CD 9 Xp 9 1 pC 9 Rb > C 9 UGD 9 XppES 9 p ) 


(Cor j> 
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CoRJ=UGD ( i ) *GAM + uGp ( 1 > t ) -k-GA.m+i « 0 

i f ( cor J o lt o i o o ) cor j= i - o 

compute correction for centre body induced velocities (Copcb> 

r A ( i ) = o ■- 0 

Call GAiMCYl (CO»Xp>i »PA$U6 »U» xpres » UGP , r p d , d A ) 

IP Cal 

CALL \/TXRq;G(CD»XP»i p<_ »RA $ C »U60 » XPRhS »P ) 

C 0 R C b = U G D ( ] ) •’f'GAM + OGp ( 1 » 1 ) #-GA*' + i c U 
DO 1 A 2 J = 1 5 M Z 
C o k C b = C o R C B + G V ( . | ) / ? o 0 
if (CORCiioLTo 1 - o ) COR Cl = i e 0 
XPRES ( l ) sDUM 
RA ( l )=DUMM 

CDP=CD*RRP 

CALL VTXpjsjG ( Cop » Xp * m2 js- ? C » UGD » XppES >F ) 

I rt = o 

SUMGsOo 
DO 42 J=1 »MZ 

SUmG=Su m g+GV ( j ) 

DO 45J*1»M2 

u N v *l- + ( ULj( J) *GA« ) + l UCsD C j ) *6 AM ) +UqD ( J ) #CCRJ+UCB ( J ) *C( RCB+SUmG/? a 
w U»*;G a P Uf’i V ( 0 ) 

DEL V = 0 V ( j ) — Li i 4 V 
UELVaDtL V/UNV 
DELV =ABS(DELV) 

I F ( DEL V— EPS )45>45»44 
IRT*IrT+i 

pV(j) = {UV( j)+UmV)/? 0 

J = M Z 

L 'N v »lc>+ (UG( J) *GAM ) + (UGD( J) *GAM ) +UQD ( J ) -*COR J+UC6 ( J ) *CCRCB 
u E L V = u V ( Z ) —U \i V 
DELV=DELV/UniV ' 

UbLVsAeS (OE lV) 
u V l i'll ) = ( UV ( N2 ) +UNV) / 2 - 

IF(DElV-EPS)46>46*47 

I RT= I RT +1 
I F { I RT ) 50 5 g0 9 kO 
•'•'TImE = nT IME +1 
i R ( NT I ME- 5 i } -5 i s 60s 60 
CALL OUTPUT 
- r ( N T I M t— 5 1 ) g pi v 6 9 s 6 S 
E • R=0 

CP To 31 

I F ( m E R R } 7 0 ? 7 G 9 9 q 6 
- A I M T 6 00 9 GAy. 

format c j x,* gam—* » £i *> 07 

PkInT 60 is(C(I) 9 I=is 6 ) 

F OR, MAT ( 1 x 3 * C { I ) * 4 X * 6 L 1 5 o 7 ) 


DUC Olf SC 
DUCOl c d] 
DUCOl ! 3 ' 
DUCOl! Sr 


A, 'V J 


58 


pEaDi2 u 5 NRUf\i » v.pH I sAp j s AlF 9 C pH I ( j ) » j = x » c ) DUCOl c 6 4 

ARJP=ARJ/PI/RRP 

NP AG =0 

IF ( N! P H I cGT o 5 ) NPri r = 5 
r F ( MRUM ) 2 8 5 2 9 9 2 8 

ERROR STOP 
RRlTE (59007} 

V. R I T E ( 6 > c > Q 8 5 
ST OP 
E N 0 
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c 5U0i 

SUBROUTINE: INPUT 

0 1 MEN 5 1 ON B( 6 )^BS( 6 )s^Al 6 ) »SaS(6)9P(6»6) 

D I wEnB ION C ( 6 ) * A ( 2 5 * 6 > > AS ( 2 5 9 6 > 9 D ( 6 ) * DS ( 6 ) ’ BC ( 6 ) s GS ( 6 ) » H < 6 ) 

i> I ME^piON rB( 25 ) »bp ( 25 ) »BTA( 25 ) »TCB|_D{ 25 ) » pA ( 25 ) »XppES ( 25 ) » pHI ( 5 ) 

j 1 ytNi> I ON UG ( 2 5 ) » UqD t 2 5 > »U 6 D ( 25 ) » UV ( 25 ) » UC 8 ( 35 ) » VCB (21) ♦ 

; SJGP ( 25 »25 ) »UGA I 25 > »GV ( 25 ) »GpV ( 25 ) »Cpp ( 5 » 25 ) »Cpm < 5 *25 ) 

I MENTION ALPHA (25 3 » STALL ( 2 5 ) 9 jSTl (25) 9 TALK ( 20) 

PENSION RE ( 4 ) 

"ON MZZZ 9 CD s> R 0 » R ]_ »R2*R.3»Pl »B»BS»SA»SAS9p 
C 0 : •; m ON/ N EAR1/NRUH*NBL0»NZ»M2,NPRES>IR,NTIHE»NERR*NPA 5»NPHI »NPRINT 
C r / NEAR2/C » A » AS » D» DS »SC »G$ »H 

. ' m :->o N / N £ A p 3 / p R P » X p »Z»3L0 » Rb *BR »BTA »TCBLD»TC»rCBrp» Apa* alF »XrrES , 
pA»XCB>XR »ElCBC »RMAX>PHI sCORJ» CORGB 
Common^ NEAR4/UG sUqD »UcjD»UV»UCB » VCB »UGp »UGA»GV »GAv]»GrV >Cpp ,CpM 
CO-M-ON/ NEAR 5 / aRJ»aRJP»EFS»RAD»CL»ALPHA»p*TALL»JGiLsTaL.C 


FORMAT ( 15H1. RUN nUMBEr»I5»49X»4HpAGE»I3//) 

Format ( 12H Input /5X *6shduct giiomet pVe o 0 c/o xp/c 

i T/C RTE/RP RCB/RP) 

FORMAT ( 10X 5 2GhCAMB£r COEFFICIENTS, ? 4Fl0o6// ) 

Format (5X»22 h prcp^ll^r geometry- « r 9 13 »ix»6HBi_ades/ / ) 

FORMAT! 2 5X»34HR/RP B/op BETA TH/CHD) 

Format ( / 5X s 47HDEF i n i ti on of symdqls used im tabular o'lTpuf o«»//) 
formak iox,66hr/rp radial propeller- statuon in fraction cf pr 

10PELLER RADIUS) 

FORMAT l 10X»57HB/RP PROPELLER CHORD IN FRACTION Or PRC^DLLER R 

1ADIUS) 


FORMAT! l A X»36HotTA PrIdEllEr PITCH IN DEGREES) 

format (iox, 5 uhth/*-hd prop^ll^r ^lade thic^wcss— to-chqrd ratio) 

FrrmaT ( / 5X , 22 HCEnjTEr do uY GdqiviETr Y ... « ‘> 2X51-1 lCo/C , 5X5HXC6/C ,3X3 Hr v AX /l 
1C5 9 2 X9HX { rm a X ) / C/2 5 X » 4Fl 0 , 5 ) 

FORMAT ( / 5X»*C0NVERGENCE CR I TER I ON o u EPS I LON=* ,Fl0 o 5 ) 

FORMAT! 1OX1HV9X2UHFPEE ST REAM VELOCITY) 

FORMAT! 10X,31HU T/Tal INFLO« VELOCITY) 

r 0 rm A T ( 1 0 X * » 4 3 H G A M / V STrEngTH OF INTERNAL VnpTEX CYlInCEr N) 

r G:\f-. AT i 1 OX , 34H A LPHA ANGLE OR ATTAC^sUtGREE' 3 ) 

FORMAT ( lOX ,6 °HDELTa P/q RISE IN TqTAL PRESSURE ACROSS PROPELLER NO 
1RMa|_ IZED/ 25X9 3 1H0N FREE STREAM DYNAMIC PRESSURE) 

Fu9 ! 'aT ( lOX 9 53HCTp ID ) THpUsT COEFFICIENT qN PpqpEllSr In THE DUC 


aqrvAT { lOX 94UHCTD ( p ) THRUST CoEFFECiEnT on THE DUCT) 
rC"’ M AT ( 1OX34 H.CTDP TqTAL thrust COEFFICIENT) 

FORMaT ( lOX 965HCTD ( p ) 1 THrUST CoEFFECiEniT qM DUCT InSlUDI\j 6 ppESS 
IUpd THRUST ON/25X,2?HTHc DUCT AFT OF The' PROPELLER ) 

Formate i0x»* ctddi t^tal thrust coeffecieniT including prEssl 
ire thrust*) 

FORMAT! i0X,*Cn'Dp TOTAL NoFMAl FORCE COEFFICIENT*) 

*"ORmAT ( lOX 9 * CmDd TqTAl plTCHlNG mqnEimT CotFFUCl t\|T* ) 

FORMAT! 10X,*J ADVANCE PATIO *> 

■ ORMAT ( l n X 9*J 1 RATIO OF V TO P POP ELLER TIP SPEED *) 

format 1 20A4) 
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li9 FORMAT (10X 5 20.A4//) 

r 

FORMAT ( o 3 » F8 j 3 » F lu o 4 s- f g „ 3 5 Fp. “ 4 s 2 t 2 ^ j ipE 1 2 o 5 ) ) 

[l'< FORMAT ( I 1 0 9 6 F I 0 6 ) 

2-, FORMAT ( 5X»6 ( 1PE13.6 ) ) 

243 FORMAT ( 20X 9 5 F 10 0 6 // ) 

251 FORMAT ( 20 X 9 2 F 1 0 o 6 9 F 1 0 c 3 9 F ]. 0 9 5 ) 

5 2 FORMAT ( 2 I 5 ,7Fl0 o 4 ) 

521 FORMAT ( 8F10.4) 

522 FORMA I (1^15) 

75 0 FORMATt///iOX,*NU,M- 3 tR o r BLADES In ERROR » N^lD »* * 1 5 ) 

751 FORMAT ( /// 1 OX 9 ^-C-ON v cRuEi'iC L CRITERION MUST BE GREATER than OoO*) 

752 FORMAT (///I OX ,* A ; -i G L E OF ATTACK MUST BE LESS THAM 9„D IGREES* ) 
7*3 FqrnAT ( / // 1 OX 9 -t-a u v'Af> ^i_ RmTIq mUST bii GoEATEo THA,\] 0 « 0-* ) 

c 

20 READ ( 5 » 148 ) TALK 

r tAD (5*521) CD »XP» T'.sR 0 p 9 p ' 0 R p » bps 

READ (5 » 52 ? ) R 0 9 P 1 9 R 2 ? R 3 

READ ( 5 , 521 ) El CB C 9 X C S 9 d ha X 9 X P 

pEAD ( 5 9 522 ) NAlD » ,nZ 9 mZp » I p » n p RE S 9 NPR I NT 

READ ( 5 9 521 ) ( R B ( J ) 9 J = 1 9 N Z p ) 

READ (5 9 521} (Bp ( J) , j=l,NZp ) 

READ (59521) (BTA(J) 9 J = 1 ,NZp ) 

READ( 5 ? 521 ) (TCSLD(J) »J»l»NZpj 
READ (5 9 521 ) ( XpRES ( M ) sN«l-* IP ) 

READ ( 5 9 52 0 } NRUM »NPHI »ARJ*AlF» (PHI ( J) » j = i 95 ) 

f 

C CHECK INPUT DATA 

r 

IF ( IN PH I »GT c 5 ) N P H I = 5 
IF(NBLD) 700,7*0,701 
PRINT750 »MBLD 
NERR=1 

IF(EPS)702,7O2,703 

" -2 PRINT 751 
N E R R = 1 

IF(ALF-90. 0)704, 705.705 

D R I M T 7 5 2 

NERR-l 

I Ft ARJ) 706,706,707 
-5 PRINT 753 
N r ’’ R = 1 

IF (NERR) 7 ^ 8 , 708,999 
CONTINUE 

PRINT101 ,NRUN . NPAG 
RR I NT 1 49 , T ALK " 

INITIALISATION OF SUBROUTINES 
MZZZ=0 

CA|_|_ HUS ( CD 9 XR , XCb » ElCoC , PHAX , I p 9 XppES 9 pB , RRP 9 UCb , VC 1 ) 

IFiNERR.GEcl) go To 999 

CALL AlFrng ( CD 9 50 9 SC 9 UQA , XPRES ) 
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MZZZ=1 
DO 4-0 K-l»6 
D { K ) = V C B ( K) 

DS(K) =UCB ( K ) 

CALL PROP ( M Z P » Ro > oR 9 b Ta»TcBl*^»PRP 9 R C B R p ? M Z 9 R A ) 

compute the constant portion of the inflow profile 

rDP“-OD*RR 5 

Call 0 A M C Y L ( C 0 P v X P a .m L 9 p ; d t v 6 9 „ 5 X P p E S 9 U G p 9 p p p 9 p A ) 

CALL SrCr w G (CD» X p » TC 9 m Z 9 -■ B » U p D ) 

Call HUo (CDsXp sXCbj E l b C 9 rmAX , MZ 9 XppEo , rB 9 rrp 9 UCB » VCB) 
lompUTc. THc Four - co cq^FFilIEnTo DUE T / THE TrAIlIwG VorT ICI TY 
CRBN=CD*2 0 

CALL BNCOEF (CR3N»lJ) 

NOUT = 6 
■'4=50 

NCYL=N2-1 

CALL. AnCoEF (NCYl »N »rA»XP ♦ CD 9 pop 9 NRUf-,; ,0 »BS*A*AS) 

DO 30 K = 1 9 N Z 
BR(K)=8R(K)/RPP 
RB(K)=RB(K)/RpP 
B T A ( K ) = 9 ^ 0 0— B T A ( K ) 

BTA(k)=3TA (K) /RAD 

return 

END 



63 


TC SU02 

SUBROUTINE PKL ( CD » p ) 

D I MENS I ON p ( 6 * 6 ) , A ( 6 ? 6 * 
COMMON MZZZ 


IF ( m Z Z Z ) 1,2,1 

DO 1 0 K - 1 9 6 
DO 10 L = 1 , 6 
DO 10 f<= 1,7 
A ! K , L , M ) — o . 0 

~ ( 1 s 1 3 2 ) = o 0 2 6 3 3 
( i } 3 s 2 ) = o 0 1 3 4 3 


( I 9 

5 

, 2 1 =»— 

cOOoOl 

A ( 2 , 

1 

,2 ) = 

0 0 P 3 6 6 

A ( 2 , 

2 

9 2 ) = 

0 029 87 

A ( 2 , 

4 

, 2 ) =- 

0 00304 

A (3, 

1 

» 2 ) » 

•j 0 0 5 ^ 8 

K { -2 9 

3 

, 2 ) = 

o 004 n 3 

A ( 3 9 

5 

» 2 ) = n 

0 00099 

A (4, 

1 

9 2 ) C ^ 

,000*5 

A ( 4 9 

2 

s 2 ) = „ 

O 

O 

O 

K— 4 

O 

A ( A <, 

4 

9 2 ) = 

0 00 148 

A (4, 

6 

5 P } i 

0 000-4-9 

A ( 5 

5 

1 , 2 )=. 

- 0 c 0 0 


A ( 5 9 3 s 2 ) = ^ c 0 0 0 4 9 
A { 5 9 5 s 2 ) = o 000 79 
A ( 6 » 4 » 2 ) =— 0 0 0 n 2 9 
A (6,6,2)= s 00049 

A (1 ,1,3)= 3 072 31 
A ( 1 , 3 9 3 ) = 0 0 3 6 4 4 
,A ( 1 9 5 ,3 ) =**„ Oq^04 
A { 2 , 1 , r ) = o 1 4 5 5 1 
A ( 2 s 2 , .5 ) = 3 0 o 5 2 6 
A (2,4,3)=,,, 01252 
A (2,6,3)= 0 000 06 
A ( 3 , 1 , 3 ) = 0 o 2 4 9 1 
(3,5,3)= „ 0 1 6 5 5 
4(3, 5,3) =- c 00411 

A { a- 9 1 , 3 ) =-. o o 0 n 1 6 

A ( 4, 2 9 3 ) = „ 0 004 1 7 
A (4,4,3)= e 0 0 6 /■. 9 

A(4,6,3) = Ao 002/\0 

A ( 5 , 1 , 3 } = „ o 000 1 3 
A ( 5 , 3 , 3 5 =_ o oo 2 06 
A ( 5 , 5 , 3 ) = c 003 13 
A(6,2,3)= o 0 0 o 0 2 
A ( 6 , 4 , 3 ) =— 0 0 o 1 2 0 
A(6,6»3)= o 00l97 

All, 1,4)= 1 1 Q 2 3 

A( 1,3, 4) =,o5 94 5 
A ( 1 , 5 ,4 ) = o 0G019 
A( 2,1,4)= „ 2 3 8 4 9 


7 ) , 0 ( 7 ) 



A { 2 , 2 , 4 ) => o 1 4 6 6 S 
A. (2,4,4)= A o02765 
A ( 2 , 6 9 4 ) = , 0 0 0 £ 1 

A(3, 1 »4) » o 05488 
A ( 3 » 3 » 4 ) = o 0 3 7 2 
A(3,5*4)®-.00965 
A(4,l,4)« 0 000 71 
A ( 4 , 2 » 4 ) « n 0 00 92 2 
A ( 4 , 4 » 4 } ® ,01418 
A ( 4 , 6 9 4 ) = n 0 0046 3 
A ( 3 , 1 ,4 ) = A c,00q42 
a (3, 3, 4} <*~ e 00482 
.(5,5,4)= e 00732 
A ! 6 , 1 s 4 ) =■» o 00 0 0 5 
A ( & -5 2 j 4 ) a o 00005 
A ( o 9 4 9 4 ) =— 0 0 0 2 7 8 
A ( 6 , 5 » 4 ) a o 0 04 5 Q 

A(l,l»5)a.l60l6 


Ail, 

3 , 

5 ) = 

,07Q07 

Ail 

, -> 

, 5 ) 

= 0 0 0 1 0 7 

A ( 2 5 

1, 

5)» 

c32Q32 

A ( 2 9 

2, 

5) = 

c , 2 0 6 a 9 

A (2, 

4, 

5)=- 

“0 04629 

A (2, 

6 , 

5) = 

.00035 

A ( 3 9 1 9 

5)« 

0 0 9 1 8 6 

A ( 3 9 

3 9 

5 Jar 

,06338 

A( 3, 

5 , 

5 ] =• 

-0 u 0 1 76 6 


A ( 4 » 1 » 5 ) » 0 0 0 4 0 2 
A ( 4 j 2 > 5 ) =— 0 0 1 5 4 4 
A (4,4,5)= „ 02 592 
4(4, 6 , 5 ) =^ 0 00 855 
A ! 5 , 1 , 5 ) = A , 0 0 o 7 2 
A!5,3»5 ) =-,0o883 
A ( 5 , 5 , 5 ) = .,01338 
A £ 6 9 1 ,5 } = A , 0 0 0 2 5 

A (6, 2 , 5 ) a s o 0 0 0 8 
A { fc 94 9 5 ) =— 0 OO 5 12 

A ( 6 , 6 9 5 ) = ,00814 


A ( 1 9 1 9 6 ) = ,22278 
A (1,3,6)= 0 1 0 6 3 4 
A{ 1,5,6) = ,00536 

A t 0 1 4 < _ , ■ r~ , 


A / O 

H l C. tj 

1 

,6) = 

0 4 4 5 5 6 

A ( 2 , 

2 

>6 ) = 

A 0 Q " r 

0 w 0 w 

H ( £ , 

4 

*6 )=A 

,08513 

^ ( 2 9 

6 

>6 ) = A 

,00059 

A ( 3 , 

1 

6 ) = 

,17105 

A ( 3 , 

3 

6) = 

,12346 

A i 3 9 

5 

6 )=„ 

.03345 

A(4, 

1 

6 5 = 

0 0 0 2 u 

A(4, 

2 

6 )=o 

,02784 

M ■( 4 9 4 i 

6 ) a 

,05692 

A { 4 a 6 < 

6 ) =- 

o 0 1 9 3 2 

4(5, 

li 

6 ! = o 

00079 

A ( 5 , 

3 * 

6 }= n 

,01858 

A ( 5 

9 - 

> 96 ) = 

,0293 

A f 6 , 

7 


,00123 

4(6, 

2j 

6 )=a 

.00011 
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» 

A ( 1 3 1 ? 7 } = <3 2662 4 
A ( 1 , 3 , 7 ) = 0 l2l46 
A { 1 5 5 f> 7 ) = O 0l2l4 
A. (2,1,7) = o 53248 
A (2,2 ? 7 }* 0 33891 
A { 2 ,4 »7 ) “** o 1 1 8 9 s 
A l 2 so »7) * *^U'->429 
- { - , 1 ,7 )« o 245 3 3 
A {3,3,7)= ,18468 
: ( 3 j 5 s 7 )•»*•& 06 2 1 0 
. i - , 1 , 7 / » „ 0*4-662 
■ (4,2,7 )- fto 03871 
1 ( 4 , 497 }= o 09563 
a (4,6, 7}= Aio 03416 
A (5,1,7)= 0 00738 
~ { 3 , 3 , 7 ) = r 02992 
A S3, 5,7)= o 05 270 
s ( c , 1 s 7 ) » A o 00 194 
A ( 6 , 2 , 7 ) = „ o OUQ 88 
A ( 6 , 4 , 7 ) =— .. 0 1 9 1 8 
A ( 6 , 6 » 7 ) = e 0 2 9 9 8 

C( l)=0o0 

c C 2 )=0 c 25 

C ( 3 ) =0 „ 5 n 
C ( 4 ) =0 o 75 
C ( 5 ) = 1 o 0 n 
C ( 6 ) = 1 o 5 Q 
C(7)=2o0n 
GOTO 71 


DO 20 J = 1 , 7 
IF ( C ( J ) —C D ) 20,19,21 
CO* TI NO L 
GO TO 21 
M = J 

GOTO 100 
. ’ = J— 1 

M = J 

CELT =C(N)-C(M) 

DIFF =CD-C(M) 

0 E L T A = D i F F / D E L T 
DO 30K= 1 , 6 
DO 30 L = 1 , 6 

P-t K 9 L ) = A ( k , l » M ) + ( kt£ L T A* ( A ( K * L » N ) — A ( K » L » M ) 5 } 
continue 

go TO 71 
DO 40 K= 1 5 6 
DO 40L = 1 , 6 
9 ( < , L ) = A ( K , L , M. ) 

CONTINUE 

return 

END 


A 1 1 

p 1 

9 7 5 

as 

0 2 6 6 2 4 

A< 1 

9 3 

9 7 ) 

£Z 

0 12146 

A ! 1 

*5 

97} 

a 

0 0 1 2 1 4 

A (2 

,1 

,7) 

s 

,53248 

•-(2 

9 2 

,7) 

*0 

38691 

A ( 2 

9 

9 7 } 

=— 

,11.898 

A ( 2 

v C 

9?-) 

=» - 

- no /. 7 0 

A ( 3 

*1 

5 JL 

97} 


,2^5 3 3 

A ( 3 

9 3 

,7} 

- 

- 1 8 *+ 0 n 

A (3 

9 5 

97) ( 

06 214 

A { 4 

9 1 

9 7 ) 

a 

n 0466 2 

A ( 4 

9 2 

,7) 


,03871 

A (4 

9 4 

9 7 } 


0 09563 

A (4 

96 

9 7 ) 

“ r\ 

o 0 3 4 1 6 

A ( 5 

$ 1 

97) 

= 

,007 38 

A ( b 

5 3 

97} 


,02992 

A ( 5 

9 5 

97) 

= 

0 05 270 

A ( £ 

9 1 

,7) 

®n 

A n -) O /, 

0 i -C 

A ( 6 

92 

97) 

= a 

0 000 88 

A ( 6 

94 

97) 

=5— 

,01918 

A { 6 

96 

97) 

= 

c 02998 

CU 

) = 

0 0 0 



CC2 ) 

= 0 0 

25 


C < 3 } 

=o„ 

5o 


C ( 4 ) = 

0 c 75 


C ( 5 

) = 

1.0 

n 


C ( 5 5 = 

1 0 5q 


C( 

7) 

= 2 . 

On 


GOTO 

71 



r\ r. 

L J W 

2 

0 J 

=1 

,7 


IF (C(J)-CD) 20,1 9» 21 
C0‘‘T I NuL 
GO TO 21 

*■■*= j 

10 TO 100 

'lT =C(N)-C(M) 

O::-:- =CD“C (M) 

. T~T A=DI FF/DElT 
■ * 30k= 1,6 
• ' 30 L=l,6 

3 ( K 5 l ) =A { k 9 1_ , ^ ) + { D£.(_T A* C a < <. , L ,N ) -A ( K » L , M ) ) ) 

continue 

GO TO 71 
00 40 K= 1 , 6 
Cu 40|_ = 1 5 6 
F { .< »L_ ) =A ( K ,L »M ) 

CONTINUE 

RETURN 

END 


mmwmmmmmmmm 
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| c f T 0 3 

SUBROUTINE ELLlPS (A<Sq»Tk*TE) 

sudo ell i ps — .table look up og Elliptic integeralS 

DIMENSION CKK( 100) ,CK(100) >CE(10 a ) 

I COM M 0 N M 2 Z Z 
[ if (MZZZ) 73 > 10,3 
p continue 

CkK« ARGUMENT qF ELLIPTIC INTEGRALS 
i C<K( 1 ) =0 o 00 

CKK(2)=0 oa 1 

C'<K(3)=0 o 02 


CkK( 4)=0 o 03 

CKK(5)=0 c ,0A 


CKK { 6 ) = 0 e n 5 
CKK { T ) =0 0 n 6 
C K K ( 8 ) = 0 „ n 7 
CKK ( 9 ) =0 0 A 8 


CKK 1 10 ) = A 0 09 
CKK{ll)= n olO 
CKK ( 12 ) =0 o 1 1 
CKK { 13 ) =0 - 12 
CKK(14)=0,13 
CKK.( 15) =0 o 14 
CKK ( 16 ) =0 o 1 5 


CKK{ 17 ) =0 o 16 


CKK(18)=0ol7 
CKK(19)= ao 18 
CKK { 20 ) =0 o 1 9 
CKK(21)= AO 20 
CKK ( 22 ) =0 o 2 1 
CKK { 23 ) =0 o 22 
CKK { 24 ) =o „ 23 
CKK< 25 ) =0 o 24 
CKK(26)=0 o 25 
CKK{27)=0 o26 
C K \ ( 2 8 3 = 0 o 2 7 
O ''(29J=0c28 
CKK { 30 3 =0 o 29 
CKK{31)=0 o 30 
CKK( 32 ) =0 o 31 
CKK ( 33 ) = A 0 32 
CKK( 34) = AO 33 
CKK( 35 ) =0 o 34 
CKK ( 36 3 „ 35 

CKK( 37 ) =0 o 36 
CKK ( 38 } =0 o 37 
CKK( 39) = 0o 38 
CKK 140 ) =0 0 39 
CKK(41 ) =o o 40 
CKK{ 42 3 =0 o 41 
CKK{ 43 } =o o 42 


CKKt44) -0 o 43 
CKK C 4-5 ) =0 o A-A- 
CK'K ( 46 ) = 0 o 45 
CKK ( 47 ) =0 o 46 
C<K ( 48 ) =0 o 47 
CKK { 49 ) =0 o 48 
CKK(50) =^ 0 49 
C<K(51)=^o50 
CKK.( 52 ) =0 o 5 1 
CKK{53)= A o52 
CKK ( 54 ) =0 o 53 
CKK( 55 ) =0 o 54 
CKK( 56 ) =0 o 55 
CKK ( 57 ) =0 o 56 
CKK(58)=0o57 
CKK { 59 ) =0 o 58 
CKK { 60 ) = A o 59 
CKK ( 61 } =0 o 6 0 
CKK ( 62 ) =0 o 6 1 
CKK(63)= ao 62 
CKK ( 64 ) = A r 63 
CKK( 65 ) =0 o 64 
CKK(66)=0c65 
CkK< 67 ) =0 o 66 
CKK ( 68 ) = A c 6 7 
CKK ( 69 ) =0 0 68 
CKK ( 70 ) =0 o 69 
CKK ( 71 ) =0 o 70 
CKK ( 72 ) =0 o 71 
CKK ( 73 ) =0 o 72 
CKK ( 74 5 = A o 73 
CKK ( 75 ) = A o 74 
CKK ( 76 ) =0 o 75 
CKK{ 77 ) =0 o 76 
CKK ( 78 ) =0 o 77 
CKK { 79 ) =0 o 78 
CKK{ 80 ) =0 „ 79 
CKK ( 81 ) =0 „ 80 
CKK ( 82 } = 0 * 8 1 
CKK{ 83 ) =0o82 
CKK ( 84 ) =0 o 83 
CKK ( 85 ) =0 o 84 
CKK{ 86 } = A0 8§ 
CKK < 87 ) =0 o 86 
CKK ( 88 ) =0 c S7 
CKK{ 89 ) =0 c 88 
CKK{ 90 ) =0o 89 
CKK (91 5 =0 0 90 
CKK{ 92 ) = A 0 91 
CKK ( 93 ) =0 0 92 
CKK{ 94 ) =0 „ 93 
CKK( 95 ) =0 0 94 
CKK ( 96 ) =0 o 95 
CKK ( 97 ) =0 o 96 
CKK{ 98 ) =q 0 97 
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CKK. ( 99 ) =n ° 98 
CKK!l00)=0 o 99 

ck=complete elliptic integral of first kind 

C k'(l)-I°57079& 

Ck ( 2 )=1°574746 
CK ( 3 ) = 1 o 5787A-Q 
rK(A)-l 5 5b278o 
(X {5)-l 0 S86S6S 
C :< ! 6 ) B 1 c- 5 9 1 0 0 3 
CK(7)-1o595188 
CK ( 8 } “1 ° 5994^3 
CM9) = lo6037l0 
CK(10) = 1o608049 
CKUD=1o612441 
CK (12 ) = 1 o 6 16 8 89 
C< (13) = 1o621393 
CK(141=1o625955 
CK ( 15 ) = 1 o 6305 76 
CK(16)=lo635257 
Ck(17)=1 o 64000^ 

CK ( 1 8 ) = l o 644 80 6 
CK { 19 ) = 1 o 649673 
CK ( 20 ) = 1 o 6546 1 7 
CK(21)=1 c 659624 
CK ( 22 ) = 1 0 6647 q 1 
Ck(23)=1 c 669S50 
CK(24)=1 o 675073 
CK(25)=1 o 680373 
CK.(26) = 1 o68575a 
CK ( 27 } = 1 o 691 20 8 
CK(28)=1 0 696749 
CK(29)=1 o 702374 
CK ( 30 ) = 1 0 7080 87 
CK(31)=1 0 7136S9 
CK(32)=1 0 7197 s 5 
CK(33)=1 0 725776 
C;< (345 = 1^731865 
CK ( 35 3 =1 c 7 3 80 5 5 
CK ( 36 ) = 1 „ 7443 5 1 
CK ( 37 ) = 1 o 75 07 54 
CK(38) = 1' 0 757269 
CK(39}=1 0 763898 
CK ( 40 ) =1 o 770647 
CK Ul) = l *777519 
CK ( 42 ) = 1 0 "7 gio 1 9 
CK(43)=l 0 79165o 
CK. (44) = 1*798918 
CK(45)=1 o 806328 
CK(46)=1 0 813884 
CK(47)=1 0 821593 
CK(48)=1*82946„ 

CK ( 49 ) = 1 0 y 3749 1 
CK ( 50 )=1 *345694 
CK(51 5=1.854075 
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CK(52| = 1o 862641 
CK<53)=1o8714(U 

C<(54)=1o880361 
CK(55 l = lo 8S9&33 
CKt56)*l» 898925 
CK C 57 ) = 1 * 90 3647 
CK<58) = 1 0 918MO 
CKt 59 )=1.92 85 26 
CK (60 ) = 1, 938908 
CK t 61 1 = 1 e 9495 68 
C<(62 ) 0 1 o 96 05 21 
CK ( 63 ) = 1 0 971783 
CK f 64] = 1„ 98 33 71 
CK ( 65 ]*1 0 93 53 03 
CK £ 66 ) = 2 n 0 75 98 
CK £ 67 j * 2 0 02 0279 
CK f 68 I = 2 0 03 33 69 
CK ( 691=2 o 046894 
CK170 ) = 2„ 060&82 
CK.(71 ! = 2 o 075363 
CK(?2)=2„090373 
CK ( 73 1 B 2 o 105948 
CK.£74) = 2 0 122132 
CK t 75 )»2 0 138970 
CK(761=2 0 156516 
CK(77)=2 0 174527 
CK ( 78 ) = 2 0 193971 
CK(79)»2. 214022 
Ck( 80 ) »2« 235068 
CK C 81 | = 2. 257205 
CK l 82 ] = c 0 28 05 49 
CK t 83 ) = 2 o 305232 
CK ( 84 ) “ 2 0 33 1409 
CK( 85 ) = 2 0 359264 
CK t 86 ) 3 2 0 3 8 90 16 
CK £ 87 ] * 2 „ 42 09 33 
CK ( 88 ) - 2 0 45 53 38 
CM89) = 2 0 492635 
CK ( 90 ) » 2 0 53 33 35 
CK { 91 ) = 2 0 5780 92 
C !/ £ 92 ) = 2„ 627773 
C:< [ S3 1 = 2 0 68 35 5 1 
C^ ( 94 ) = 2 0 7470 73 
CK C 95 ) = 2 0 820752 
CK ( 96 ) = 2 0 90 83 37 
CK £ 97 1 = 3 „ ^ 1 6 1 12 
CK(S8!«3 0 155875 
CK £ 9 9 ) « 3 0 354141 
CK(100}=3„695637 

ce=Compl£te Elliptic inTEgerAl o f second <4«d 
CE( 1 ) =1 „ 5 70796 

CE t 2 } = 1 0 5 66 862 

CE(3)al 0 56 2 9l3 



c f(4)=1°558948 
C f ( 5 ) = 1 , 5 5 4 9 6 9 
CF(6)=1c550973 
rc{ 7 ) = 1 0 546963 
C E £ 8 ) 3 1 o 542936 
ci { 9 ) . 53 8 S 93- 

CE(l0)=i o 534833 
r, .EUl) = lo530758 
12) = 1 0 5 26665 
13 1 = 1 = 5 2255 50 
=1 = 5 18428 
( 1 5 ) = 1 . 514284 
IF ( 16)“1, 510122 
17 ) = 1 0 505942 
. ' : ( 1 8 ) = 1 o 5 0 1 7 4 3 
i •; 1 9 ) = i o 4 9 7 5 2 6 
.. • ! 2 C } = 1 o 4 9 3 2 9 0 
CE £ 2 i ) = 1 c 4890 35 
:E!221=1o4847610 
CE(23 )=1 o 480466 
CE(24)-1 0 476152 
C E f 2 5 ) = lo4718l8 
CE(26)=1c^67462 
CE ( 27 ) = 1 a 4630 85 
CE(28)“1 0 458688 
CE(29)=1 0 454269 
CE ( 30 ) = 1 o 449 8 2 7 
CE(31)=1 0 445363 
CE ( 32 ] = 1 o 440876 
CE( 33 J = 1 0 436366 
CE(34)=1 0 431832 
CE(35)=l a 477274 
CE ( 36 1 = 1 o 4-2269 1 
CE(37)=l o 4l80s3 
C£{38)=1 o 413450 
CE(39)=1 o 408791 
CE(40 )=l o 404l05 
CE(41)=1 0 399392 
CE ( 42 } = 1 0 3 94652 
CE(43)=1 c 389883 
CE( 44 } = 1 0 3 85086 
CE(45)=1 o 380259 
CE £ 46 ) = 1 c 3 7540 2 

C E £ 4 7 ) = l o 3 7 0 5 ]_ 5 

CE(43)=1 0 365596 
Ct(49)=l o 360645 
':E(50} = l o 355661 
C :! 51 ) = 1 0 350644 
CE £ 52 )=> 1 0 3455g2 
CE ! 53 )*1 0 3405 q5 
CE £ 54 ) = 1 o 3 35 3e2 
CE£ 55 ) ® 1 0 33022 3 
CE { 56 ) = 1 0 325024 
CE( 57 ) = 1 0 3 197 8 8 
CE( 58 ) * 1 „ 3 145 1 1 



CE(59)=l c 3C9l92 
CEf60)=lo303832 
C-(61)=1c293428 
GE(62)=1. 292979 
CE ( 63 } = 1 o 2 57484 
CE(6 4)“1 0 281942 
CE(63)=1o27635a 
CE (66 ) = 1 * 27070 7 
CE(67)=1 c 2650i3 
CE (63 ) =1 o 259263 
CE(69) = 1 o 2534-5 8 
CE(7C} = ls24/293 
CE(71) = Ic24167.1 
CEi72)=I»235634 
CE(73)“1 0 229632 
CE(74)=l a 2235l2 
CE(75)=l 0 21732l 
CE176)=1o211056 
CE ( 77 ) = 1 o 2047 14 
CE(78)=1, 198290 
CE(79)=lol9l78l 
CE ( 80 5 = 1 o 1 85 1 8 3 
CE f 8 1 ) = 1 o 1 78490 
CE(82)=1 0 171697 
CE(83)=1 0 164798 
CE(84)=1 0 1577S7 
CE(85)=1 U 150656 
CE(86)=1 0 143396 
CE(87}=1 0 135998 
CE ( 88 ) = 1 0 128451 
CE(89)=1 0 1SS742 
CE(90)=1 o 112856 
CE l 91 ) =1 o 104775 
CE{92)=1o096478 
CE(93) =l on 87938 
CE(94)»l o 079l2l 
CE( 95 ) =1 „ 069986 
CE ( 96 ) = 1 0 060474 
CE(97)-1.050502 
CE ( 9 6 ) = 1 „ n 39 947 
CE(99)=l ari 28595 
r r i 100 ) -1 0 0l5994 

'0 TO 30 

{ AKSQ— e 99)20>20j21 
1 ARA=0.2 5* (loO-AKSQ ) 

! -ST =1 o OOE-O 7 
IE (PARA- TEST) 701 ,702 »7o2 
PARA =TEST 
2lP=AL0G(4o/PARA ) 

Tk =ZLP^0 o 5* ( 1 o + PARA ) “PARA 
TE-1 o 0+ ( 2LP*PaRA )— PAR A 
GO TO 30 

ja=ioo„o*aksq 

J A = 1 + j A 
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IF (C«K O AJ—AKSQ )* 21-2 » 23 5 22 
TK=C K( JA ) 

TE=C EC JA ) 

GO TO 3 0 

C0N= CAKSq-CKK( JA I W (CKK t JA+1 )-CKK'( JA) ) 
TK=CK (JA)+CON«- COKI JA+ 1 ) -CK ( JA ) > 

• TE=C E( JA ) ft CON*(«CC CjA+1 )~CE( JA) ) 

GO TO 30 

3 IF (AKSQ- on 1) 221 »r.2l 9 72 q 
1 D£RA=o 23*ASSQ 
GO TO 70 0 
0 AKSQ»1 o- AK ?.?. 

GO TO20 
CONT IMUE 
RETURN 
END 
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-FTC 04 

SUBROUTINE LAMBDA (XP»ZP»YP) 

taslT lookup of heuman lambda function 

DIMENSION Y ( 1 9 .» 1 9 ) » X ( 1 9 ) » Z ( 1 9 ) 
COMMON N.ZZZ 
IF ( MZZZ ) 20 > 10 5 ^0 
CONTINUE 

X=ARC SIN K( DEGREES) 

X(l)=0 o 00 rt 000 
X{~)=5o AAA 000 
X ( 3 ) =19 oaa 00uq 
X ( 4 ) =1 5 o 0 A 000 0 
X ( 3 } -20 o 0 A 000 0 
X ( 6 ) -25o AA OOOn 
X { 7 ) =30 o OC OOOq 
X { a ) =35^000000 
X { 9 ) =40 o 000000 
X( 10 ) =45 0 000O0 A 
X(li)=50 o 0000n A 
X(12)=55 o 00000a 
X( l3)=60 OA 000b A 
X( 14)=65 0 A 0000 A 
X{ 15) = 70 oa 0000/v 
X ( 16 ) = 75 a r 0000 A 
X { 17 ) = 80 0 n 000Q A 
X{ 18 ) = 85 o 0000n0 
X(19)=90 o 00C00 a 

Z=BETa (DEGREES) 

Z{ l)=0 o 000000 
Z ( 2 ) =5 o 000000 
K 3) =10 o 0 a 000q 
Z( 4) =15c 0 a 000q 

Z ( 5 ) = 20 0 a qOQo a 
Z(6)=25o AA 0000 
It 7) =39 o 0 a 0000 
Z ( 8) =35 o OOOOOn 
Z(9) =40 <= 00 OOOq 
K 10 ) =45 a A OOOn A 
Z{ 11 ) = 50 OA 0000 A 
Z( 12 ) = 55 o 0000q a 
It 13 )=60 o 00000 
Z( 14 )=65 c QOOOnO 
Z(15 )=7Q o 000U00 
Zt 16 } = 75 0 A OOOq A 

Z{17)=8Q oa 0000a 
Z(18)=85 (>a o00o a 
Z( 19 )=90 o 0000r'0 


y=heuman lambda function 


y ( 1 , 1 ) -0 onOOCOn 
Y(l»2)* 037156 

Y(1,3)=0o 173&48 
Y(1,4)=0 o 258819 
Y ( 1 5 5 ) = o 342020 
Y(l»6) 3 -=422&18 

Y ( i 9 7 ) = o 5 
Y(1,8)=o573576 
Y(1,9)=o642788 
Y(l,10) -o707ln7 
Y{1 9 11)=,76&044 
Y(l 9 12)=. 819152 
Y ( 1 * 13 ) = o 866025 
Y(1s 14)= o 906303 

Y ( 1 s 1 5 ) = 0 9 3 9 6 Q 3 
Y(1,16)=o965926 
Y ( 1 s 17 ) = o 934808 
Y(l,18)= 0 996l95 

Y ( 1 » 19 ) = 1 o 

Y ( 2 » 1 ) = o 0 

Y{ 2» 2 ) = o 086990 
Y(2,3)=. 173318 

Y { 2 9 4 ) = o 2 5 8327 
Y ( 2 s 5 ) = o 34 13 70 
Y(2,6 J»o42l8l5 
Y ( 2 » 7 ) = o 49 905 0 

Y ( 2 , 8 ) = o 5 7 2487 
Y(2,9)=o 641567 
Y(2,10)= o 705765 

Y ! 2 s 1 1 ) = o 7645g2 
Y!2,12)= 0 8175 
Y(2, 13)=. 864388 
Y(2>14)= o 904599 
Y(2, 15 )= ,937930 

Y ( 2 , 16 ) - 0 964135 
Y(2, 17) = , ,983037 
Y(2»18 ) = . 99462-4 
Y ! 2 s 1 9 ) = 1 0 

Y ( 3 , 1 ) = 0 0 

( 3 g 2 ) = ,086495 
Y(3,3)= 0 172332 
V{3,4)=. 256858 

Y ( 5 s 5 ) =0 c 3 39430 
Y(3 9 6)=0 o 4]_94^9 

Y < 3 , 7 ) =0 0 496 2 1 9 
v (3 5 8)=0, 569244 

Y i 3 , 9 ) =0 0 637940 

Y(3,10)=0 o 70l786 

Y<3,ll )=Ao760298 

Y(3 »12}=„ o 813034 

Y ( 3 , 13 ) = 0 , 85 9602 
Y(3,14)= no 899660 

Y<3,i5)=o,932934 
Y(3,16)a0 o9 5 9 244 
Yf3,l7)=o, 978597 



Y ( 3 

9 1 8 ) 

a 0 , 99 1 5 1 1 

Y ( 3 

,19) 

=1 o0000 A 0 

Y(4 

9 1 ) = 

o A o o o n A 

Y(4 

9 2 ) = 

0 o a 85677 

Y (4 

9 3 ) = 

0, 1707Q4 

Y ( 4 

,4) = 

0 0 2 54434 

Y ( 4 j 3 ) = 

0,336231 

Y ( 4 

9 6 ) “ 

0, 415475 

Y(4 

9 7 ) = 

0,491565 

Y ( 4 

9 8 ) = 

0,563926 

Y ( 4 

9 9 ) = 

C 0 6320.10 

Y { 4 

,10) 

= 0 0 695 3 a'7 

Y ( 4 

.■11 ) 

= A s 75 3 346 

Y ( 4 

,12) 

= 0 08 057^3 

Y ( 4 

,13) 

=a=S52010 

Y ( 4 

,14) 

= ac,891969 

Y { 4 

,15) 

=0*925384 

Y ( 4 

,16) 

-0 0 952226 

Y ( 4 

,17) 

“a°972787 

Y ( 4 

,18 ) 

=0o98e015 

Y ! 4 

,19) 

= 1 0 OOOQaO 

Y (5 

,1) = 

0,00000a 

Y ( 5 

9 2 ) = 

0,084549 

Y { 5 

, 3 ] “ 

0,168458 

Y ( 5 

9 4 ) = 

0,251092 

Y i 5 

9 5 ) = 

0,331827 

Y(5 

9 6 ) = 

0 0 4 1 00 5 4 

Y { 5 

9 7 ) = 

0,485184 

Y(5 

.81- 

0,556657 

Y ( 5 

» 9 ) = 

0,623939 

Y(5 

9 1 0 ) 

= 0 0 686540 

Y ( 5 

,11 ) 

=0,744012 

Y ( 5 

,12 ) 

=0,795963 

Y ( 5 

,13 ) 

= 0 0 842073 

Y ( 5 

,14) 

•0. 882 119 

Y ( 5 

,15) 

=0*916018 

Y ( 5 

9 16 ) 

=0,943 Q 18 

Y ( 5 

,17) 

=0,966343 

Y ( 5 

,18 ) 

= a 0 984410 

Y{5 

,19) 

= 1 0 OOOOaC 


Y(& 9 1)=0, A 0000 A 
Y{6,2)=0, 083124 
Y(6,3)=0 o 165625 
Y(6,4)=0 o 246S82 

Y (6,5 )=0 o 3 262,8 8 
Y(6,6)=0,403252 
Y{6,7)=0 0 4772o3 
Y { 6 » 8 ) = 0 0 5 47 6 0 0 

Y < 6 , 9 > = O o613936 
Y(6»10) =0 o 675748 
Y{6,11)= Ao 732623 
Y(6,12)= rtO 784220 
Y ( 6 , 13 ) =0 „ 8302 62 
Y(6,14) =0,870676 

Y (6» 15 ) =0c 905441 



76 


Y(6,16)=0„934867 
Y(6,17)=0o9596a7 
Y{6»18)=0o930779 
Y{6,19)=1„00C0a0 
Y ( 7 , 1 ) =0 „ nOOOo a 

Y(7,2)=0 on 8l425 

Y ( 7 , 3 ) -0 o 1622^-7 
Y(7»4)“0 o 241870 
Y (7»5)“0 e 319707 
Y(7»6)=0 o 395l9l 

Y ( 7 * 7 ) “0 o 467777 

Y ! 7 j 3 } = 0 o 5 3 6 9 5 3 
v (7,9)=0 o 602244 
V(7 s l0)=0 o 663225 

Y ( 7 9 1 1 > = 0 o 7 1 9 5 3 3 
Y(7,12)=0. 770883 
Y(7,13)=ao817q93 

Y ( 7 s 1 4 ) = 0 o 8 5 8 1 1 7 
Y ( 7 , 15 ) =a o 8940 95 
Y(7,16)=0o925409 
Y(7»17)-0o952751 
Y(7,18)»0o977159 
Y(7*19)“l o 0000,,0 

Y ( 8 s 1 ) =0 0 C'OOOoO 
Y { 8 » 2 ) =0 0 079476 
Y{8»3)=0 o 158377 
Y<3»4)=0 o 236l34 
Y(8,5)=0 o 3l2l92 
Y(8,6)=0 o 3860l3 
Y(8,7)=0 o 4570p. 6 
Y(3,8)=0 o 524935 
Y{8»9)=0 o 589127 
Y(3sl0)«0 .,649233 
Y(8,11)=0o 705094 
Y(8,12)=0 o 756337 
Y(8,13)=0 o8029a3 
Y ( 8 » 1 4 ) = 0 o 8 44 8 2 0 
Y(3,15)=0o882297 

Y ( S » 1 6 ) = 0 o 9 1 5 7 5 7 
Y ( 8 » 1 7 } = A c> 945 8 7 3 

Y(8,18)= ao 973573 

Y ( 8 9 1 9 ) = 1 o 0 0 0 a 0 
Y (9, 1 ) =0 „ OOOOoO 
Y(9,2)=0 o 0773r>7 

Y { 9 9 3 ) =0 0 1 54U 73 

Y{9,4)“0 o 229767 

Y(9,55=0 o 303869 

Y(9»6)=0 o 3758p a 

Y ( 9 , 7 ) -0 0 4453 30 

Y(9»81=0 o 5H786 

Y(9,9)=0 o 574862 

Y{9,10)=0 o 634231 

Y{9,11}=c : =689642 

Y'(9, 12 ) =0 C 740932 

Y(9,13)= ao 788d51 
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Y < 9 , 145=0. 831085 
Y(9, 155=0. 870277 

Y { 9» 16 ) =0.906056 
Y(9, 175=0. 939042 

Y ! 9 , 18 5 =0 . 970039 
Y ( 9 9 19 5=1.0000^0 
Y ( 10 » 1 5 = 0 . 0 00 n A 0 

Y { 10 5 2 ) - 0 o 0 7 4 9 5 3 

Y ( 10 ,3)=„. 149408 
Y(10,4)= 1 A o222873 

Y ( 1 0 9 5 5=0.294 & 84 
Y( 10 96 >=0,364976 
7(10,7 5=0. 432729 
7(10,8 ) *0.497760 
Y(10»9)“0o559735 
Y( 10 ♦13)*0 o 6 18381 
Y{10»ll)= A ,67350l 
Y( 10*12 72498 5 
Y(l0,l3)= no 77283Q 
Y{10»14J* fto 8l7l55 

Y ( 1 0 , 1 5 ) = A „ 8 5 3 2 1 7 
Y(10»16) =0.896419 
Y!10,17)= A . 932311 
Y(10*18)»0r 966575 
Y(10»19)*1o000a00 
Y(ll*l)=0o0000 A 0 
Y(ll, 25=0.072455 
Y(ll, 35*0.144464 
Y(ll,4)=0o215587 
Y(ll, 55=0.235399 
Y(ll ,6 ) = 0 . 3 5 3 5 A 0 
Y(ll,7)= A . 419519 
Y( 11 , 85=0.483126 

Y ( 11 *9 5 =0 o 544o38 
Y(ll9lU)«0 o 6U203S 
YlllslDoO, 656976 
Ylll, 12)=0 c 708785 
7(11,135=0.757496 
7(11,145*0.803241 

Y (11,155=0.84626? 
7(11,165= A .886942 
Y(ll, 17)» rt , 925731 

Y ( 1 1 , 1 8 ) = A , 9 6 3 2 0 4 
Y( H ,19 ) = 1 o 000„0C 

Y ( 12 , 1 5 = 0 o 000 000 

Y ( 1 2 , 2 ) = A ,, 0 6 9 8 6 .1 
7(12,3) =0.139334 
7 ( 12 ,45=0.208034 
7(12, 55=0.275597 
7(12,65=0.341676 
7(12,75=0.405956 
V(12, 85=0.468167 
7(12,95=0.528076 
Y(l2»10)= no 585512 
7(12,11 }= Ao 64o369 
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YP?»12)=^o692612 

Y{12,13)»^o742291 

Y(12 9l4)=0 = 739- 1 3 / 
Y(12»15)**0r.834576 
Y(12»16)*»0 o 877717 
Y(12»17)*0 o 9l9353 


Y C 12 s18)»ao959944 

Y ( 12 s 1 9 3 « 1 o O'Oq a 0 0 

■.Ml3,l)=no0000«0 
Y( 13 *2) =0-067226 


/( 13,3)=0,.134i26 

Y i 13 » 4 ) = a o 2 0.0 ? S 0 

Y ( ].3 »5') =0 = 265664 
v{13»6)=0. 329751 
' ( 13 »7 ) =0 = 392328 
Y ( 13 »3)=0 o 45 3 192 
Y ( 12 » 9 ) =0 = 512 167 


Y(13,10)= A „ 569122 
Y!13 5 11)=a» 623985 
Y{13»12)= ao 676745 
Y(13»13 )=n<,727455 
Y(13 9 14)= Ao 776237 
Y(13»15)« rt =823283 
Y<13, 161=0=863346 

Y ( 1 3 » 1 7 ) »0 c 9 1 3 2 4 0 
Y{13»18)»0 e 956626 
Y(13»19)*1 o000a00 

Y ( 14s 1 ) =0 = G0O0 a 0 
Y(14)2)=0 o 064614 
Y(14,3)=0 o128968 
Y(14 9 4)=0=192s a 9 
Y ( 14 * 5 ) =0 o 255 897 
Y{14»6)=0o318Ca 9 
Y(14,7)= ao 378946 
Y(14,8)=0 o 438541 
Y(14,9)=0 s 496661 
Y(l4»10)a0o553214 

Y ( 1 4 s 1 1 ) »0 o 6 0 8 1 5 3 

Y i 1 4 j 1 2 ) =0 o 6 6 1 48 0 

Y ( 14 s 1 3 1=0 = 713246 
Y( 14,145=0 = 763552 

Y ( 1 4 > 1 5 ) « 0 o 8 1 2 5 5 2 

Y ( 14 9 1 6 ) = A =. 860 44 3 

Y ( 14 » 1 7 ) =■ A . 907464 

Y P4 9 18 )*a = 95388 5 

Y ( 1 4 j 1 9 ) = 1 0 0 0 0 a 0 0 
Y!15,1)= ao 0000a0 
Y{15 9 2)=0o062l00 
Y(15.3)»0 c124009 

Y ( 15 94) =0 = 185540 
Y ( 15 s 5 } =0 o 2465 17 
Y( 15 96 ) =0 = 306778 
Y(15,7)«0o366180 
Y{15,8)=0=4246a4 
Y(l5»9)=0 o 48lQ59 



Y( 15 . 1 O)»*no 530 l 83 
Y( 15 jll 593247 

Y ( 15 »12 )*» A ,647159 
Y( 15 j 13 )* 0 o 699961 

Y(15,14)«0 3 75i731 

Y{15»15)a0o802581 

Y ( 1 5 s 1 6 ) a 0 o 8 5 2 6 5 4 
Y(15 ? 17)=0o902ll9 

Y( 15 > 18 )**o < 35 H 66 

Y ( 1 5 j 1 9 ) = 1 o 00 o „ 0 0 
Y ( 16 j 1 )=Qo 0000 a 0 

Y( 16 * 2 )= 0, 059779 
Y( 16 9 3)= n0 119433 
Y( 16 . 4 )= 0 o 178339 

Y ( 1 6 s 5 ) = 0 o 2 3 7 8 S 3 
Y ( 16 9 6 )= 0 o 296459 

Y ( 1 5 j 7 ) = 0 o 3 5 4 4 1 _ 
Y{ 16 * 8 ) = 0,41 18 57 
Y ( 16 9 9 ) =0 o 468 5 46 
Y( 16 9 10 )= 0, 524596 
Y{l 6 sll )= />0 579721 
Y(l 6 ,l 2)= O0 634200 
Y( 16 , 13] =„ ,6 87972 

Y ( 16 9 1 4 ) = A e 7 4 1 0 8 9 
Y(l 6 » 15 )a ft0 793624 
Y{ 16 , 16)= ft , 845669 
Y( 16 » 17 )« 0 o 897332 
Y( 16 > 18 )= 0 o 943733 

Y { 1 6 » 1 9 )=l o 00 Q n 00 
'i ( 17 * 1 ) = 0 . 0000*0 
Y( 17 , 2 )= 0 c 057773 
Y(l 7 , 3 )= 0 dl 5479 
Y( 17 * 4 ) =0 , 173054 
Y( 17 , 5 ) = ao 230436 

Y ( 1? ,6 ) =0o 287571 

Y ( 17 ,7 ) =■„, 344410 
Y( 17 3 8 5 = ao 400915 

Y ( 17 9 9 ) =0 o 4570 55 
Y(l7»10)a0 o 5l2813 
Y(17»11)»0 o 568181 
'■ (17,12) = 0 o 623166 
Y(17,13)»0, 677782 
Y(17»14)»0 o 732059 
Y(17,15 )=o o 786036 
Y(17 9 16)=^ 0 839759 
Y(17,17)o A o 893286 
Y (17 j IS ) 0 94667 7 

Y(l7»19)sl o 00c*00 

Y f 18 9 1 ) a 0OO0OO 
Y(18,2)=0o056256 

18 s3) =0 D 112490 
Yds, 4 )= 0 o 168682 
Y{18,5'!=0 0 224814 
Y( 18 96 ) = 0 , 28086.7 
Y(l8,7)= 0o 336S26 
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Y(i8>8)- An 392679 

Y(18>9)= a o4484I7 

Y(18,l f 31 ss AO 304034 

Y( 18 9 H ) °C 0 5 59529 

Y ( 18 » 12 )*0 c 6 14903 
Y(18>13 )“ A o670162 

Y ( 1 8 » 1 4 ) - 0 o 7-2 5 3 1 5 
Y(18»15 J»0. ,780373 

Y ( 1 3 » 1 6 ) *0 , 8 3 5 3 5 2 

Y ( i 8 » 1 7 ) = 0 0 8 9 0 2 7 0 

Y ( 1 8 » 1 S ) — A o 9 4 5 1 4 b 
Y ! 1 8 5 1 9 ) = 1 0 00 0 A 0 0 

Y ( 1 9 3 1 ) = 0 o 0 0 0 n ^ ' 

Y ( 1 9 s 2 ) = A o 0 5 5 5 5 6 
Y(l9»3)=0clllli: 

Y( 19 »4)=0 o 16666" 

Y(19»5) = A « 222,222 

Y ( 1 9 » 6 ) = 0 c 2 7 7 7 7 8 
Y(19>7) =0 o 333?33 

Y ( 19 >8 ) =0 > 388839 
Y( 19 »9)-0» 444444 

Y ( 1 9 9 1 0 ) = 0 o 5 0 0 0 0 0 
Y{ I 99 11 )« Ae 55 5 5 56 
Y(19»12 )« Ao 6H111 
Y{19,13 )« Ao 666667 
Y(19914)« ao 722222 
Y(l9»15)=0 o 777778 
Y|l9»16)a0o833333 
Y( 19»17)»0 o 888889 

Y ! 19 9 18 )*0 o 944444 
Y!19919) = 1 0 OOo a (jO 

60 TO 90 
DX=5 o 0 
DZ =5 o 0 
FNX =XP/DX 
7=FNX 

f! =|v+l 

f;:Z=zp/dz 

7=FnZ 

tm+i 

Cl- ( ZP-Z ( M } )/DZ 
C2- ( XP— X(N) } /DX 
i = (Y(N»M+i 5 -Y { n ,m ) ) *C1 
' = ( Y ( N+ 1 9 M ) -Y ( N 9 M ) ) #C2 

■ i = Cl*C2* ( Y (N+l »M+1 > -Y ( N+l »M ) +Y ( N »M.) -Y ( N *M„1 ) ) 

Y 'P = Y ( N » M ) A A+B+D 

RETURN 

END 
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■H pj giTC 05 . . 

SucjRQOT X fvj L 


HUB 


( C U * Xr ^X^ bsELCBCsH^IpjX 9 p 6 ? p p p ? U 9 V ) 


calculation of velocities induced by the clhTrE body 
#* centre body is represented as rankine body 


0 1 hEnsS Ion x ( 25 ) » RU ( 25 ) * U ( 2 5 ) 9 v ( ? 5 ) , Fu ( 5 „ ) 9 FV ( 5 0 ) 9 D ( 5 ) 9 DS ( 5 ) 


u FORMAT .(/// 10 X,*IMPUT CENTrEboDY DIMENSION IM ERROR*/// ) 

1 format c ///iox ,-K-suDRQbi ine hub unAolE Tq Compute cenTrEboDy gom-Tr 

1Y*///) 


IF 1 MZ 22 ) 3 * 1 9 2 
XlZs (X r-XCB) /ElCBC 

COMPUTE LOCATION Of POINT SOURCE 

IF(XIZ-H) 10 »ln »11 
NERR *1 


ERROR MESSAGE 
PRINT 700 

return 

1 H 2 =H*H 

X I Z 2 =X I Z*X I Z 
A=Xi 2 -HR 
DO 20 J = 1 9 20 
AA=XIZ*H 2 *SQRT ( A*A+H 2 ) 

AA= X I Z 2 «^ 6 qpT ( AA ) 

AA=SQRT(AA) 

DEL = A B S ( A A— A ) / A A 
I F ( D E L — oO rt l )30 9 3 0 s 1 9 
A= (AA+A ) / 2 o 
CONTINUE 

Error message 

PRINT 701 
f'ERR = l 
RETURN 
A = A A 

EMV=( ( X I Z 2 — A*A } **2 ) / l 4 “*X I Z*A ) 

RETURN 

I F ( IR)l Aft »lOo »200 

Compute velocities Induced at the duct refrEnce cylinder 
C= 0 o 5 / cd/elcbc 

NX «»-IR 

DC 40 J= 1 ,MX 
Xl=X { j )-XCB~Xi Z^E lCBC 
XI “XI/ELCBC 
XPA »XI+A 

xma«xi-a 

RXPA=Sqrt ( XPA*XPA+W*W ) 

RXMA b SQRT (XMA*XMA+W*W) 

RXPA3*RXpa#*3 



’vX^A3=RXMA*-*3 

j ( j ) = EM V* { X P A /R X P A 3-X M A / RX M A 3 ) 

V ( J ) =»EMV*W# ( 1 o /RXPA3-1 = ^ RXMA3 ) 

CONTINUE 

RETURN 

COMPUTE IMFLO’U VELOCITIES 

Xj=XR* J ' XCB“XI Z*ELCBC 
X I = X I /ELCBC 
DO 60 J= 1 » I R 

■;, S 5B ( J ) * 5 /CD/El CBC/ RRP 

XDAaXI+A 

XMA “XI -A 

i?XPA*SQRT ( XpA#XpA+'.J*W ) 

R X M A = S Q R T ( X m A *• X M A + ' ■ J * W ) 

RXPA3=RXPA#*3 

rXMA3=RXMA**3 

u t J ) =EMV* ( X P A / R X P A 3 “ X M A / R X M A 3 ) 

V ( J ) »EMV*W* ( 1 o / R XPA 3- 1 s / R XMA 3 ) 
pSIttW*W/2 c-EMV* t XPA/RXpA-XMA/pXMA ) 

CONTINUE 

RETURN 

compute dim) and d-star(n) fqurier coefficientd 

W=0o5/CD/ElCBC 
N=IR 
CN *N 

PI “3 0 1A15926 
DTH«PI/(CN-1» 5 
TH “-DTH 
DO 80 J = 1 9 N 

th =th +dth 
CSTH =COS(TH) 

-<G=0 e 5*(l ( »_t s TH) 

XI=XG-XCB“X I Z*ELCBC 
XI “XI /ELCBC 
X A = X I -f A 
X: A “XI -a 

'XP'A “SORT ( XPA*XPA+W*W ) 

"XmA=SqrT ( XwA*XmA+W*W ) 

‘ V .P; 3*RXP a**3 
XXMA3»RXMA**3 

~U(j)= EmV * ( XpA/rXpA3~XmA/RXMA3 ) 

F'v ( j ) = EMV#W* ( l 0 /rXpa3—1 o /rXmA3 ) 

CONTINUE 
NOUT =6 

CALL F0UrCS(FU»FU*N»NOUT,0') 

CALL FoUrC ( (FV *FV»N »NOUT 9 0 ) 

DO 81 J=l 9 N 0 UT 
0(J}-FV( j) 

DS(J)*FU{ J) 
u ! J} =DS ( j } 

V(J)“D( J) 

RETURN 

END 
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jDpic 06 _ 

1 SUBROUTINE CAm ; dEr (CD^TCsrp) 

suoroetinE to compute induced camber coefficients 

DIMENSION CDRAT (5) » DCD (4,5) »Rp(4) 

O.aTA CDRAT/0*0»0»25»0o50»0 o 75»lo rt / 

OATADCD ( l s 1 ) * DCD (1,2 ) »DCD(i,3 ) , DCD ( 1 , 4 ) » DCD ( 1 , 5 ) /0 * „ , v 
i 0 .. ^,0 3 1 6 9 0 o 0 0 5 0 5 , 0 o o 0 6 9 6 / 

K-rtTA t-iCt/ { 2 » 1 ) , UCD ( 2 »2 ) , ^CDf 2 , 3 ) , DCD ( 2 , 4 ) ,UCU( 2 » 5 ) / n ». 
I ) ; 22236 s 0.334 14, ()„ 43 940/ _ 

.rtlA UCu ( 3 * 1 ) 9 ECO ( 3 9 2 •) ^ ^CD ( 393 ) , DCD ( 394 .) 9 DCD ( 3 , 5 ) / «>( 
;.2..07223 3 0ol 1022 5 Oo 14923/ 

j,hTA 0 C D ( 4 , 1 ) , D C D ( 4 , 2 ) , <>CD ( 4 , 3 ) 9 D (_ l/ { 4 , 4 ) , D C L > ( 4 , 5 ) /On 
1-0 « 00539 9-0 o O0943 9-0 o 0x37 5/ 

DO 2 J=1 » 5 
N=J 

I F (CDRAT ( J ) —CD ) 2 , 3 , 4 
CONTINUE 
4 M=N-1 

OEL=CDRAT (N) -CDRAT ( M ) 

DIFbCD-CDRAT(M) 

delta=dif/del 

DO 10 K=l,4 

RP(K) = DCD ('S»M) +DEL.TA* ( DCD ( ION ) -DCD ( < ,M) ) 

.0 continue 
C-0 TO 20 
3 . M=N 

| DO 6K®1 9 4 
| R D I K ) = DCD ( K »M ) 

DO 21 K= 1 , 4 
21 RP I K ) ®RP ( K ) #Tc 
| RETURN 
END 


.001499 

»0o 1 0897 » 
90oC35589 
9— O0OO2O3 9 



RXMA3 = RX-MA**3 

,J < j )- = EMV* { XPA/RXPA3-XMA/RX.MA3 ) 

V ( J) =EMV*W* ( I 0 /RXPA 3 -I 0 XRXMA 3 ) 

40 CONTINUE 

RETURN 
C 

C COMPUTE INFLOW VELOCITIES 

C 

2 a 0 Xi=XR^ XCB~XIZ*ELCBC 

X I =X I /ELCBC 
DO 60 J=1,IR 
W = RB ( J } * 0 5 /CD/ElC BC / RRP 

63 XPAoXI+A 
XMA “XI —A 

64 rXp,A=SqrT (XpA*XP'A+W*W) 

R X M A = S Q R T { X M A *• X M A + ;1 * W ) 

RXPA3=RXPA**3 

R.XmA3 = RXMA**3 

u ( J ) =E mV* ( X P A / R X p A 3 — X M A / P X M A 3 ) 

V £ J ) aEMV*W*( 1 o / R X P A 3*“ 1 s / R X M A 3 ) 
pS I sW#W/ 2 e —EMV* t X P A / RX p A-X M A / P X M A ) 

60 CONTINUE 

RETURN 
C 

C COMPUTE DIM) AND D-STAr(N) FOURIER CoEFFlCTENTD 

C 

3 W=0o5/CD/ELCBC 

N-IR 
CN =N 

PI *3, 1415926 
, DTH»PI/(CN-lo> 

TH a-DTH 
DO 8 0 J = 1 9 N 
TH «=TH + DT H 
CSTH =COS£TH) 
xg*o„5*( i>-t; s th) 

X I =XG-XCB“X I Z*ELCBC 

XI aXI/ELCBC 
X P A a- X I + A 

X. A -XI -A 

RXP.A =SQRT (XPA*XPA+W*W) 
pXMA-SoRT ( XmA*XmA+W*W ) 

R X P A 3 ® RX P A#*3 
R X M A 3 = R X M A * * 3 

FU(j}= EMV *( X P A / R X p A 3 ~ X m A / P X M A 3 ) 

FV ( j ). = EMV*W*.( 1 o / rXpA3—1 » 7rX,mA3 ) 

80 CONTINUE 
N 0 U T =6 

CALL F0U-R.CS ( FU jFU »N sNOUT » 0 ) 

CALL FqUrC ( (FV»FV,N»NO'JT ,0) 

DO 81 J=1 sNoUT 
D £ J ) °FV ( J ) 

DS(J)-FU{JJ 
U( J)=DS( J) 

81 V ( J ) =D £ J ) 

RETURN 

END 
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SIBFTC 06 

SUbROUTlfsjE CA.mdEr (CD <> Tc 9 Rp J 

C subroutine to Compute induced ca m be r colfficiemt 

dimension CDRaT(5) ,DCD(4,5) ,RP(4 ) 


/0*0»0o25»0 o 50,0o 75 9 1 o ^/ 

uh TAiXO ( x 5 i ) »DCO ( 1 » 2 ) »0CD( i , 3 ) ? DCD( i , 4 ) ,DCD( 1 ,«)/r 

1V»«0316 »0 o 00505. 0 cJO696/ 1 

,uCU( 2»3> 9 DOj ( 2 » 4 ) »UC0(2.5) 
i. J c 22 236 s 0 0 334l 4 , 0 0 43 940/ * 5 

- * 1 ! ,ULU( 3 , 2 ) s u CD { 393 } 9 DCD ( 394 ) 9 DCD( 39 R) 

-'-'.,07223 9 0 0 1 1 u 2 2 9 0 o 1 49 2 3 / 

^JhTa DCD( 4 5l ) 5 UCD( 4 , 2 } ,oCD (4 9 ^) 9 Uuu( 494 ) ,uci>U»A) 

.1 — o - O' 05 39 9—0 oU 0 9439 — 0 u 0i375/ 


- A 9O0OO149 5 
/ rt “0 .0 «■ 1 0897 9 
/^oOsOoC 3558 s 
/O 0 J 9 — 0c00203 9 


DO 2 J=l*5 
N = J 

I F (CDRAT ( J )~CD ) 2 .3 9 4 
2 CONTINUE 

4 M = N— 1 

DEL=CDRAT ( N )— CDRAT ( m ) 

DlFaCD-CDRAT(M) 

delta=dif/del 

DO 10 K=l »4 

RP ( K ) =D0D ( S » M) +DELTA* ( DCD ( K sN ) -DCD ( if »m 1 1 
ID CONTINUE 

GO TO 20 
3 . M=N 

DO 6K= 1 94 

6 R P ( K ) - D C D ( K » M ) 

2a DO 21 K = 1 9 4 

21 RP(K)=RP(K)*Tc 
RETURN 
END 
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SIBFTC 07 

SUBROUTINE PROP(NZp »RC5»0R.»bTA»TCb L D»RRp» R CBRp,NZ > RA) 

c subroutine Tq compute prorEllEr geometry parameters 

c 

DIMENSION rB( 2 5) »Br ( 25 ) s'TCBlD(? 5 ) »-X'« 25 ) *Y( 2 5) ? Z ( 2*^ 5 >24(25) >dTA(25) 

1 »RA(2'5) 

L ■ * 

P=RCBRP#RCRRP 

D=l--P 

X ( 1 ) =RCBRP 

ZN=NZ 

M'Z = N Z + 1 

DO 20 J=2,MZ 

K = J-1 

AK.=K 

E=AK/ZN*DaP 
2 o X(J)=S0RT(E) 

DO 21 J=1»MZ 
21 RA(J)»X(J a 1) 

X ( 1 ) = ( RCBRP+RA ( 1 ) )/2o 
DO 22 J = 2 s M Z 

22 X(J)=(RA{J)+RA(J— 1) J/2» 

■ K® 1 

RB (n2P+1 ) =RB ( mZp ) 

BR(NZP+1)=BR(NZP ) 

BTA ( NZP+1 ) =BTa (NZP } 

T CB|_D ( N.Zp n i ) = TCBlD(nZP ) 

DO 39 J = 1 » N Z 

3^ I F ( RS ( K ) ^ X C J ) ) 3 1 , 32,33 

32 Y(J)=BR(K) 

Z( J)=BTA(K) . 

ZZ( J)=TCBLD(K) 

GO TO 39 

31 K*K+1 

GO TO 30 

33 DEL =X( j)-rB(k-I) 

DEL=DEL/ t RB ( K) — RB ( K— 1 ) ) 

Y(J)=DEL*(BR( K ) — B R ( K-l ) ) 

Y(J)»Y( j) a br(<-1) 

Z( J)=DEL*(BTA'<K)-BTA(K-l) ) 

Z(J)=Z(J) a BTA(K-1) 

ZZ( j) = DEl*(TCBlD(k)-T.CBLD(k-1) ) 

Z Z £ J)=ZZ( J)+TCBLDC<-1) 

39 CONTINUE 
DO 40 j = 1 > M Z 
RB ( J ) =X ( J ) 

BR(J)=Y( J) 

BTA ( J ) =Z ( j ) 

TCBLD( J)»ZZ( J ) 

40 CONTINUE 
RETURN 
END 
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BIBFTC SU08 

SUBROUTINE CLAlF (j) 

c 

C COMPUTATION OF BLADE SECTION LIFT COEFFICIENT 


r 


r 


c 

lOi 


1 

r 

r 

V. 

c 


c 


c 

2 

19 

20 
21 


DIMENSION TCB( 10) ,CLMX(1C ) 


DIMENSION 
D I MEnS I ON 
DIMENSION 


B ( 6 ) sBS (6 ) sSa ( 6 ) sSAS ( 6 ) »P ( 6 » 6 ') 

rB( 25) »Br (25 ) »BTA( 2 F ) sTCBLDC ? 5 ) a pA ( 2 5 ) » Xpp.ES ( 2 3 ) » dHI ( 5 ) 

ALPHA(25T »STALL(25) »JSTl( 25) ,TALk(2 A ) 


^Common mezz > cd » R u 9 pi 9 p 2 s rs 9 p 1 » b sbs a sa , sas , p 

common /n£ari/nrun»nBld»n2»mz»npres»ir, n time,n£rr»npag»np^i »nprint 

common /neap 3/ ppp 9 xp s z 9 Bld s rb S sp »bta »-tcbld aTC * rcbrp 9 apa»alF *xpr£s» 

1 . RA » XCB » XR » E’LCBC » RMAX » PH I » COR J » CORCB 

COMMON /N Ear 5/aR J » AR JP » EPS » RAD » CL » ALPHA » STaLL a JSTL 9 TALK. 


FORMAT!//////////* THE BLADE THI CkNESS^TO-CHORD RATIO IS OUTSID 
IE RANGE .0,0 To „34#) 

I F ( MZZZ ) 2 a 1 9 2 
CONTINUE 


TAclE OF VALUES qF ClMAX VErSUS BlADE THICknESS_Tq Ch )rD rATIq 

TCB( l)=0 o 0 
TCB(2)=0 o 06 
TCB(3)=0 o 08 
TC8(4)=0 o 10 
TCB(5}=0 3 12 
TCB (6 ) =O f , 15 
TCB(7)=0 o 18 
TCB { 8 ) =0 0 2 1 
TCB(9)=0„24 
TCB(10)=0=34 

CLMX ( 1 ) = p, o 9 
CLMX { 2 ) = A o 9 
CLMX ( 3 ) =1 <3 2 
C LMX ( 4 ) = 1 o 45 
CLMX ( 5 ) =1 0 6 
CLMX ( 6 ) = 1 o 5 
CLMX ! 7 ) = 1 . 35 
CLMX (8) *1.3 
CLMX (9 ) =1 0 25 
CLMX (10)«l o 1 
RETURN 

CONTINUE 

IF (TCB(IO)-TCBLD(J) )2l»2C,20 

IF(TCB(1)-TCB L D( j) ) 22’22a21 

PRINT 10l 

NERR*1 

RETURN 
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N-l 

IF<TCB(N) n TCBLD< J))24s25,26 
N = N + 1 
60 TO 23 
CLMAX«CLMX (N ) 

GO TO 27 

BEL = ( I<-BLD(J)-TCB{ N -1 ) ) / (TCB (N ) „TC6 (N-l ) ) 
0 L M A X = C L M X ( N— 1 ) +DEL* ( C(_MX ( N ) „CLMX ( N — i ) 1 

continue 

C L = 2 „ * P I * A L P H A ( J ) / R A D 
STALL (J)®0„0 
IF(CL-CLMAX)29>30»30 
CL=CLMAX 
STALL ( J ) =1 *0 
CONTINUE 

return 

E N D 
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SIBFTC Su'09 

SUBrqUT I NE ARC'S I w ( ARG>pH I s INDEX ) 


subroutine for computing arc sime 


arg is the sine of the angle 
dHI IS The ANGLE (principle value ,-pI/2 To pi/2) 
output from SUBROUTINE* 

INDEX IS 1 if ANGLE IS TO BE IN RADIANS. 

INUEX IS 0 IF AnGlE IS Tq BE IN DEGREES « 

HRl = lo 5707963 

A = 1 „ o*»arg*arg 
I F ( A } 5 ♦ 2 , I 
IF(ARG) A, 3 ,3 
PHI “HP I 
GO TO 6 
PHI --HP I 
GO TO 6 

PH I =ATAN ( ARG/SQR T ( A ) ) 

GO TO 6 
W R I T E { 6 » 1 0 0 ) 

00 FORMAT ( /5X*ERROR<. *oS'IN X GREATER THAN i„0*»//) 
GO TO 2 

IF( INDEX) 7 ,7, 8 
PH 1*90 o *PH I /HP I 
RETURN 
END 



-J O' 
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SIBFTC SU09 

SUBROUTINE ARC'S IM (-ARG » PH I s InDFX ) 
c subroutine for computing, arc sime 


c ARG is the sine OF the angle 

cdhi is the angle (principle value ,-pi/2 to pi/?) 

; OUTPUT FROM SUBROUTINE* 

C t m f! f y X t c S n ^ t c F ANGLE 15 To 5E IN RADIANS* 

c index is o ii- angle is Tq be in degrees* 


HP I = 1 o 5707963 
A = 1 „ 0«»ARG*ARG 
I F ( A ) 5 9 2 9 1 

2 I F ( ARG ) 4 , 3 , 3 

3 PH I “HP I 
GO TO 6 

4 PHI “—HP I 
GO TO 6 

1 PH I =ATAN { ARG/SGR T (A) ) 

GO TO 6 

5 WRITE{6»100) 

100 FORMAT(/5X*ERRORo.,oSIN X GREATER THAN 1*0*,//) 
GO TO 2 

I F ( INDEX) 7, 7 ? 8 
PH 1*90 o*PHI/Hpi 
8 . RETURN 
END 
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SIBFTC SU10 

SUBROUTINE matrix { a » NR ) 
DIMENSION A ( 6 9 7 ) > B ( 1 9 » i 9 ) ] ?.) 

NC = NR+1 ' ' 

DOl J== 1 s NR 

I N t J ) =0 

DO 2 I ® 1 9 NR 
CON = 0 « n 
I 1 = 1-1 
DO 3 J = 1 9 N R 
JJ=J -1 
I F ( I I )4,4,5 

5 DO 6 K = 1 9 I I 

I F ( ( < ) ) 6 ’ 9 6 

6 CONTINUE 

4 CONA=A ( J 9 1 ) 

IF ( CONA ) 8 9 9 > 9 

8 CONA«-CONA 

9 I F (CONA-CON ) 3 5 3 9 11 

II CON=CONA 
N ( I ) *J 

3 ' CONTINUE 

I F { CON ) 5 no , 60o s 500 
600 PRINT 601, f 

6 A 1 FORMAT ( IX 9*SINGULAR 9 I**,I 3 ) 
STOP 

500 NN=N(I) 

DO 12 J = 1 9 N R 
12 A I J 9NC ) =0 o 0 
A ( NN 9 NC ) = 1 „ 0 
D IV=A ( NN , 1 ) 

DO 13 L = 1 9 NC 

13 A ( NN 9 L ) =A ( MM 9 L ) /D I V 
DO 14 L= 1 , NR 

IF (L-NNJ 15 9 14 9 15 
15 CMUlT =A ( L 9 1 ) 

DO 16 J= 1 , NC 

16 A ( L 9 J ) = A ( L » U 5 — CMU|_T#A I NN * J ) 

14 CONTINUE 

DO 17 L = 1 , N R 
DO 18 J = 1 $ N R 



A ( L , J ) -A { L 9 J+l ) 
CONTINUE 
CQNT INUE 
DO 19 J= 1 , NR 
DO 20 L= 1 , NR 
NN=N(J) 
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$ I BFTC SUil 

SUBROUTINE FOURCS ( F »B s N».MOUT,(S!PR.INT) 

r > i ' 

C SUB ) FOURCS— —FoUR I er cosine series determination 
C FROM RALSTON AND WllF {-MATHc METHODS FOR DIGITAL COMPUTERS' 

C CHApTo 24 ZZ GoGoErTZE L 

DIMENSION F ( 5 o } , F F ( 5 0 ) , B ( 5 0 ) 

COMMON MZZZ 
700 FORMAT (lH/ 5 

7 a I FORMAT ( 4£15 „ 7 ) 

702 FORMAn5X^THETAr, 6 X^FOURIER F N * 6 X*HOR I G I NA L FN* 4 X*F 0 Ur I ER C 0 EF */ ) 

IF (NOUT-N-) 10 9 10,11 
11 NOUT =N 

10 ZR0=0 o 

ONE»loO 
T WO = 2 e 0 
IZRO»0 
I ONE =1 
I TWO =2 
HP I ® 1 o 5707963 
NH B N 

NS* I TWO* I N H"I ONE) 

FNS«NS 

FR=TWO/FNS 

PIN»PI*FR 

DO 201 1*1, NH 

CJ“ I"I0NE 

CS I =C0S ( p I n#C J ) 

CCSl=CSI+CSI 

CA*ZRO 

CB=ZRO 

DO 202 J * 2 , N S 
JK=NS--J+lTWO 
IF(JK-NH)220,220,221 
221 JK» I TWO*NH~JK 

220 CC=F ( JK)+CA*CCSI~CB 
CB»CA 

202 CA»CC 

FF(I )=FR*(F(lJ+CA*CSI-Cb) 

I F ( I -I ONE 5 290,291,290 
291 FF( I )=FF( I ) /TWO 
290 JN=NS—I«*I + IONE 

IF(JN) 203,204,204 

203 FF( I )»FF{ IJVTWO 
2a^ CONTINUE 

201 CONTINUE 

IF(NPRINT 0 1)2,6,7 
6 M*N 

GO TO 1 
7 M=NOUT 



f\J LO 4> 


• 90 


1 cn = n 

PRINT 700 
PRINT 702 
DTH=P I / ( Cl!— 1 o 0 ) 

T H“— DTH 
DO 3 I»1,N 
TH=TH+DTH 
SUM a 0 o 0 
DO 4* J = 1 , M 
CJ=J-1 

SUMbSUm+FF ( j ) *CqS ( C j#JH ) 
PRINT 701»TH,SUm»F( I ) »FF( I ) 

DO 8 K = 1 9 NOUT 
8 P(K)=FF(K) 

RETURN 

E N D 
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-“T r suiz 

SUtnnUT I Snrovfv C Cl) » XS'C ,np ) 


v ;| aC IT I ES IvOUChU BY Dl STo I • j T I nn nF ?nUPCE hImGS 
1 J ’ '• S I on T*- ( ) •»^D (?”5 ) »BT { -») $ A ( £ } »Of>< ) »UnD( •?*•- ' 

f r\i t. :/-% f . t * * “7 

' t = : - ... I 41^n?.5 
* f i \ = j r ? n o /•- q 
' >' 2 \ =— T C-“- 1 , ?6 

' 3} =-10*% ‘US 

■■ ' < s . i =TC*? „ P43 
'• ( I' ) =-TC*i 0 °1. " 
o r - o w e 

n T T = 7 o a 

1) pT=p I / B n 

xs=ycc-,5" 

• v.v + j 

, '= , ' , + 4 

15 J=1 

r«f.n=n B 

15 rs OfJ)=° 0 

Bo IP J = ? , v 

TH( J ) seTt-l ( J_t )+PT 

1?~ T rTU( j) 

5 T = 5 T M ( T ) 

5. \! 2 =5 TV ( T/?„ ) 
or. ( j i = a ( i /cm.? 

oo 18 T=2 ,5 
' . ’ T - T — 1 

r “ ‘ ~ T M 

T f * ! )< = T M + I M ^ 2 
iv y = j \? y 


{ J ) =00 ( J s 4-FV* M T ) * ( V.‘2*tf Cfv > 
r n m j j vy p 

or ( j ) = on ( j J *.$T 
rr tj vi in- 
i’ ' ( l '.'+ 1 ) =ot + 1 0 


? ? 


C> 1 


i-p 


nr{M + i )=nr> (M) 

( 1 ) = A ( 1 I 
T“t 1 )=0n 
r T O = O j /DT t 
T r> — a 

a 

'-'-OTP/S .. 

U T =o , 

C c T = -2 o ,#Y C 
C5T2=CST#CST 

nr £Q JP=1 ,'t.rrj. 

r=rn (J r , 

" = 1 
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ZWP*CD* ( CTP-CST) 

PP=(ZWP*ZWP) + ( ( R + 1 .. )#{r+i 0 ) j 
PM= ( ZWp*ZWp ) + ( (p-i„ ) 

AKZ = 4 0 -JfR/pp 

CALL ELLIPS ( AK2 »Zk jZE) 

PP-SOR T ( pp ) 

U=ZWP*Z E/PP/Pm 

52 IF(TP-TH{NT))55,5.4»53 

5 2 NT »NT + 1 

’30 TO 52 
5'3 0 = Q D ( N T ) 

00 TO 56 

jp 0.njt=oIu 1> + <TP “ TH(NT “ l),#<OD<HT, “ QD(NT - 1 ) ,/<TH(NT) - T HiNT«i)) 

60 Td=Tp+DTd 
<«2 

UI=UI+H*(0I(i)Ho*Dn2l+DI(3]| 

D I ( 1 ) = 0 1 ( 3 ) 

T F ( Tp**p I ) 5l»5l , 

61 UV=U I -&CD/P i 
CqD { JR ) =UV 
TP=0 o 

U I =0 o 

69 ' CONTINUE 

RETURN 
END 



n n 
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SIBFTC..SU13 

bobROuTlNE VTXRNG (CD»Xi>C»Ip»PP»CN»UG»XpRES»p) . 

C VELOCITY INDUCED BY DlSTnibUTlOM nF VopTfx’ RlwGS 

C 

DIMENSION TH( ?°5 ) »GD( 205 } »DI (^ ) srp ( 25) ,UG( ) 9 Cn ( 6 ) »XprES ( 2 5 ) 
DIMENSION FU(5^J »FV(5«)»CP(6) ,P(6*6) 

r 

COMMON MZZ Z 
C 

C 0 = C N < 1 ) 

C 1=CN { 2 ) 

C2«CN(3) 

C3=CN(4) 

C4=CM (5 ) 

C 5 = C N ( 6 ) ‘ 

SET Up A TASLE OF GAM: A 0 /GAMMA* S I M ( THETA ) VERSUS THETA 
C 

DG=100„ 

DTI-200* 

PI“3ol415926 
13 DT=pi/DG 

N = DG 
NM-N 
N-N+l 
K = N+4 

DO 12 J=l,< 

TH( J)=0 o 
' S T = 0 o 0 
12 GDI J)=0 o 

DO 19 J=.2*N 
TH( J)=TH( J-D+DT 
120 T=TH(J) 

C T = CqS ( T/ 2o ) 

CT = CT/SIn { T/ 2 o ) 

a=co*ct 

S T = S I N ( T ) 

3=C1*ST 

C=C2*S I N { 2 „*T ) 

D = C3*SIN ( 3o*T ) 

F=C4*SIN ( 4 0 *T) 

F=C5*SIN(5 0 #T) 

GD ( J ) = A+B+C+D+E+F 
GD( J)=GD( J)*ST 
19 CONTINUE 

TH(N+l)=PI+ls 
GD ( N+l ) = GD ( N ) 

27 GD ( 1 ) =2 o *C0 

TH( 1 )=0 U 

31 DTP=PI/DTI 

TP = 0 o 
H=DTP/3 o 
UI=0 o 

C 

mmmMm 
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c 

10 IF Cir > 70, n ,11 

11 XS=XSC“i 5 
CST*-2»*XS 

C ST2 = CST#CST 

c COMPUTE INFLOW VELOCITIES 

40 DO 6 9 J R » 1 9 I R 

R=PD( JR) 

50 0 = 1 

' T = 1 

^ - '■ ■ ■ . ' o 0 J = K 3 3 

5 7 p = S I M ( Tp ) 

CTP-COS (TP ) 

CTP2=CTP#CTP 

X 2 = C S T 2*-4 2 „*CST*cTp ) +CTP2 
X2=X2*CD*CD 
PP=X2+ { (R + l 0 )*(R+l e ) ) 
dm=X 2+ ( (p^l,)*( p«*l ,, ) ) 

AK2c4 0 *r/dd 

CALL ELLIPS(-AK2,ZK»ZE) 

PM=2»*(R-lo )/PM 
PP“1 ./SqrT (PP) 

U=Z«-ZE*( 1 o+PM) 

U=U#RP 

52 I F ( T P*-*T H { N T ) ) 5 5 » 5 4 9 5 3 

5 3 NT =NT + 1 * 

GO TO 52 
54 G=GD ( NT ) 

GO TO 56 

6 5 G=oD ( mT — i ) + ( Tp ^Th ( nT—i ) ) # ( Gu t m T )' A GD ( mT_i ) ) / ( Th ( m T ) _Th ( m T-j ) ) 

56 DI(J)=G*U 

60 TP = T P+DT p 
K = 2 

>JI=UI+H* { D I ( 1 ) +4 c #D I ( 2 ) +D I ( 3 ) ) 

01 ( 1 } = DI .(3 ) 

IF(Tp-pi ) 51 >51961 

61 UG( JR) =CD/2„./PI*UI 
TP“0* 

U I =0 „ . 

69 CONTINUE 
GO TO 999 

C C/vipUTE VELOCITIES AT THE DUCT pEFrEkiCE CYLINDER 

r . .. ' 

v_. 

70 I R = ~ I R 

DO 99 JR= 1 » I R 
X = XPRES ( JR ) 

u - 1 
IS ~ J- 

N T = 1 
XS c X-„ 5 
CST=-2^*XS 
CST2=CST*CST 
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77 DO 71 J=K»3 
CTp=COS(Tp) 

CTP2=CTP*CTP 

X2 = CST 2*-*4 2 o -*CST*cTp } +CTP2 
X2=X2*CD*CD 

PP=X2+4 0 

AK2=4„/PP 

CALL ELLIPS ( Ax2 »2 k »ZE) 
p d =i 0 / sort { pp ) 

U = ( Z K- Z E ) * P P 

72 if (Tp^Th(\tT) )75>74s>73 
73 NT-NT+1 

GO TO 72 
7 4 G - G D ( M T ) 

GO TO 76 

75 G.sGD(nT«! ) + (Tp-Th(nT-i J ) * ( Gu [ |mT ) -GO ( NT-j ) ) ' ( TH ( ( \,T )*■!,-: Z,,!-! ) ] 
76 DI(J)=G*U 

71 T P = TP+DT P 

K = 2 

0 I =U I + H* ( D I ( 1) A 4 e *Dl (2)+DI (3) ) 

D I <1 )-DI (3 ) 

IF(TP~PI }77,77>78 

78 UG< JR)=CD/2 0 /Pl*Ul 
TP = 0 . 

UI=0 o 

99 CONTINUE 

GO TO 999 
C 

c compute f <n) and f-star(n) EqUrier coefficients 

c 

200 MZZ=MZZ? 

MZZZ=1 

DO 210 J = 1 9 6 
210 CP(J)=0=o 

DUM=0o 0 
DO 211 L“1.6 • 

211 DUM=DUM+CN ( L ) *P ( 1 , L } 

CP ( 1 ) = ( CN ( 1 ) -DUM ) /2 o 0 
DO 212 K = 2 «6 
DO 213 L= 1 » 6 

213 CP ( K ) =Cp ( « ) +CN t L ) *P l K »L ) 

212 CP ( K ) = < CP ( k ) -CM < K ) ) / 2 « 0 

N=IR 

EN=N 

D T A a P I / ( E N~ 1 o ) 

TA=— DT A 

DO 220 JR “ 1 » N 

TA=TA+DTA 

CSTA=COS(TA) 

CSTA2*CSTA*CSTA 
DUM® 0 o 0 
DO 215 K=2 3 6 
J-K-l 
TJ = J 

DUM®DUM+CP ( K ) -*COS ( T J*T A 5 


215 
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20y 


202 

203 

204 

2 0 5 
206 ' 
201 


208 


220 


216 


999 


FV( JR) =CP ( 1 )+£>Um 
< = 1 ' 

NT = 1 

DO 20l j*sk» 3 
CTP=C05 ( T p ) 

CTP2=CTP*CTP 

■X2“CSTA2*-4 2o*cSTA*C'Tpj+CTp‘? 
X2=X2*CD*CD 

P P“X2+4 a 0 
AK2=4 0 /pp 


CALL E L L I P S ( A.k; 2 s Z << j 2 E ) 
PP^l o/SQRT (pp j 
U = ( ZK W ZE ) *pp 
I F ( T P-^TH (NT) ) 7 0 5' •> 2 04 , 2 03 

nt=nt+i 

GO TO 202 
G=GD(NT) 

GO TO 206 

GeGD ( M-T«j_ ) -r { T.p_Thj ( (v.j T _ - 
DI ( J)=G*U 

tp=tp+dtp 

K = 2 


* t G D ( N T ) — G D ( t\ T _ j ) ;• / ( T h ( kjT ) _T H ( n T-j j j 


U.I»Ul+H*-(DI(l)+4o*0I (2)+DT(3 i ) 
D I ( 1 ) = D I ( 3 ) 

IF CTP-PI )207,207»208 
f r U(JR)=CD/2o/pl#UI 
TP=0 cO 
UI=0 o 0 


CONTINUE 
NOUT =6 


NPR = 0 

CALL FoUrCS (FU»FU»n.MOUT» NPP) 
DO 216 J«l,6 
UG? J)=FU( j) 

UG ( J+6 ) ®CP (J) 

MZZ7 =MZZ 


return 

END 



r— I vO h CG 
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$ I B F T C SU14 

SUdroUT I nE AlFrnG ( CO 9 1 R » SC 9 UgA » XrrES ) 

C AXIAL VELOCITY INDUCED by ALPHA VORTEX RlNGl 

r 

DIMENSION TH(r 05 ) »GA(.205) »DI < 3 ) ,XppnS( 2 = ) sSI 
D I MEWS I ON TCD [ 7 } 5 TSC ( 6 5 7 ) 

DIMENSION GS ( 6 ) » F ( 5 0 ) sHt"' (6 ) s F'H( 5„ ) 9 DVT { 3 ) 

C 

common mzzz 

DATA TLD/0o»0o5'90o75»1o0»1o5,2o«*5o0/ 

DATA TSC< 1 ,j ) » TSC (2 »D ,TsC(3»t> »TSC< 4 9i) ,TSs 

l0 o ) 0 o 9 0 <, 5^ 5 )0 C / 

umTA T5D ( 1 , 2 ) ,TSC(2 »TSC( 3 , 2 ) *TSC( 4 , 2 ) ,T_ 

1 9 — 0 o 2382 , A 0 c. 0 1 7 5 »0c > 0 o , 0 „ / 

umTA TS<-{ x , 3) , TSC ( 2 9 3 ) ,Ts.C( 3 » 3 ) 9TSC( 4 , 3 ) , I .• 

1 9 — 0 c 345 y** "» 0 0 042 s 0 o >0 e 9 0 a / 

DATA TSC ( i 9 4 ) 9 ISC ( 2 » 4 ) 9 1 S C ( 3 9 4 ) 9 T SC ( 4 9 4 ) 9 T . 

1 9—0 o 41 5 1 9 H — 0 oOaqO 9 0 o jUo 9 ^ 3 / 

DATA ToC( ^ , 5 ) 9 T S C ( 2 95^ 9 T^C( 3 ,^) ,TSC( 4 9 ^) , 1 s ! 
1 9-0 o 4834,-0 o 1230,-0 0 OI 2 I 9 Oo , 0 t / 
dmTA TSL{ 1 , 6 .).,TSC(2»6J *TSC(3 #6 ) , TSC (4,4) ,Ts 
1 9-0 o 5045 , A 0o 1666 9-0o0 3 l7 ,o» AA 35 > 0 C / 

DATA TSC ( t 5 7 ) » T SC ( 2 » 7 ) 9 T SC ( 3 9 7 ) 9 T SC ( 4 , 7 ) ,TS’- 
1 9- A o4590» a 0 o 26^7 9 —0 ol2469 A 0o rt 49490 c / 

C 

IF (MZZZ) 2 , 1,2 

1 continue 

COMPUTE SC(N) COEGFICIEwTS 

DO 11 J= 1 ,7 
JSC-J 

I F ( CD“T CD ( J ) ) 4 9 3 9 1 1 
CONTINUE 
DO 5 N = 1 , 6 
SC ( N )=TSC(M9 JSC) 

GO TO 10 
I F ( JSC-1 J 6 ,6,7 
RETURN 

DELSCa (TCD(jSC)-TCD(jSC-i ) ) /(TCD(jSC)- CC) 
J- JSC 

DO 8 N= 1» 6 

SC(n)= < Tsc(N » J- l 5 -TSC ( N 9 j ) ) /DELSC+ TSc ( n » v 
SET Up m TABlE qF GAmMA AlpHh /V*SI n ( THE.T a ) 

c 

10 Pl= 3 o 141 5 926 
DG~ 200 o A 
DT= PI/ DG 
<■ 201 

DO 12 J = 1 9 K 
TH( J)= 0 o 0 
12 GA { J ) * 0 o 0 

DO 20 J= 2 9 K 


I t 6 } » USA ( 2 5) 9 ST (6 ) 


(5»i),TSC(69i) / lo 

< 5 » 2 ^ » TSC ( 5 , 2 ) / 1 c, 
- ( 5 9 3 ) , TSC { 6 , 3 ) /( ■ . 

' 6 r 4 ) 9 l SC l 6‘9 4 ) / - S 
( 5 * 5 ^ 9 TSC ( 5 , 5 ) / 0 o 
(5^6^ 9 T S'*- ( 6 9 6 ) / 0 & 
(597) 9 T SC ( 5 9 7 ) / 0 ... 


•5 

VERSUS THETA 
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1 3 


1 4 

2 0 


C 

C 

c 


47 


C 

C 


Th(J)= t h ( J— 1 ) A dt 
T= TH { J ) 

cT N* CoS(T/ 2 .0)/ S i isj ( T/? c O) 

CT« COS(T) 

ST ( 1 ) « SIN (T) 

ST(2)= 2 o 0*S T ( 1 ) #CT 
DO 13 N= 3j 5 
NM= N-.1 
M M M ® |\|— 2 

ST ( N ) ® 2 o 0*ST ( n,M ) -x-CT-ST ( mm»/ ) 

G A { J ).= SC ( 1 ) * CTN 
DO 14 N= 2 s 6 

■3A ( J ) “ GA ( J ) + S C ( N ) * ST(N-l) 

G A ( J ) = GA(J)+ ST { 1 J 

continue 

GA ( 1 ) = 2 0 Sc ( 1 ) 

TH ( K + l J * P 1+ i O o 
GA ( K+l ) = G A ( K ) 

compute g stap(n) and hin) foupier coefficient" 

M Z Z Z * 1 ' 

DTI« 200 o 0 
DTP« pi/ OTI 
H= DTP /3 0 p 
TP= OoO - 
UI= OoQ 
VlTa OoO 

R“ IcO 

R 32* R**lo5 , 

RP2= (F+l„0)**2 
N- I R 
CN= N 

DTA» PI/ (CN-loO) 

TA= —DT A 
DO 50 JR* 1 , n 
TA“ TA+ DTA 
CSTA= cost TA) 

K- 1 
NT= 1 

DO 41 J= K. 3 
CTP= COS(TP) 

X I = ( CT p*— cST a } * CD 
X2-XI*XI 
DEN* X2 + RP2 
AK2= 4 o 0#r/ DFN 

CALL ELLlPSt AK2» Zk» Zc) 


AKP* l,0 ft A K 2 
AKQ* 2 0 /*> AK2 
AK4= AK2* AK2 

AYEss ( t 8 o 0 *r*AKQ/Aka ) + { 2 oU*r } +2 o rt -Ao0/AK2* t " 
AYE* AYE+AKP''AK2*(8o*R+4o-i6o*p/AK2)*Z|< 
u “ AYE/Sqr T (DEN)/(X2 + t R“1 oO ) * ( p**l « } ) 


f- R ‘r 1 o } } ft L E 



AX- SORT ( A K 2 } 

AK.Z.® 4 o 0* r/ rp2 

SBET A= SqrT< ( icO- AkZ ) / ( 1 0 0-Ai<2 ) ) 

CALL ARCS In ( S BETA > BETA, 0) 

CALL A R L S I N ( AK » ASK* 0) 

CALL LAMBDA ( ASK, BETA, HLMB) 

C 

EPSR® 1,0 

A V T = 4 o 0 * ( 2K-ZE ) /AK+2 « 0*P I.*AK/SQRT ( AkZ. ) * ( 1 c ) * 

1 SqrT ( { i o 0-.akZ- ) / < AKA-AK2 ) ) 

AVT * AVT*Xl/{8oO*PI*R32)+ EpSp/ 4 o 0 

42 IF(TP— TH(NT) )45, 44, 43 

43 NT= NT+ 1 
GO TO 42 

44 G=GA (NT ) 

GO TO 46 

45 G= GA(nT*-i)+ ( Tp-Tn (nT_i ) )* ( GA < NT } -*GA ( N T- i j) / ( Th ' N >■••';',•!( NT-- ■ 

46 DI(J)«G*U 

DVT(j)=G*AVT 
41 TP aTp+DTp 

K=2 

U I ®U I +H* ( D I <l)«4o*DI(2)+DI (3) ) 

D I (1 )®DI (3) 

VIT=VIT+H*(DVT{1)+4«*DVT(2)+DVT<3) ) 

DVT ( 1 ) =DVT ( 3 ) 

IF (TP“PI ) 47,47,48 
48 F < JR ) =“CD*U I /2 o /P I 
FH( JR)«CD*VlT 
TP=0 o 0 

Ul~0o0 

ViTb 0 o 0 

50 CONTINUE 
NOUT =6 
NPR a O 

CALL F0UrCS(F,F»N*.N0UT,NPR) 

CALL FOURCS (FH»FH*N»N0UT*NPR> 

DO 51 J*1»N0UT 
GS ( J 1 =F ( J ) 

HF( J)*F.H( Jt 
UGA( J+6)*HF(J) 

51 UGA( J)®GS< J) 

: ZZZ=0 

return 

2 CONTINUE 

TP = 0 o 0 
Ul=0 o 0 
C 

c calculate the axia l velocity induced by the ap 0 ve vortex distribu 

C I T I ON 

c 

In DO 99JR=1,IR 
X®XPRES( JR) 

K*1 
NT = 1 
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XSsX- 0 5 
CSTa-2'..*XS 
R = 1 o 0 

77 DO 71 J=k»3 
C T P=COS ( TP ) 

X2=(CTP-CST)*CD 
X 2 = X 2 * X 2 

DEN=X2+ ( R+l „ )*(R+1„ ) 

A,<2»4 0 *R/DEN 

CALL ELLIPS ( AK2 sZk »ZE ) 

AKP=1oO-aK2 
AKQ=2 o 0-aK2 
A L,4=AK2*AK2 

AVE=( (8*0* R *Akq/Ak4)+(2«*R)+2o0-4o/Ak2*(2c*p+1o ) )*ZE 
A ^ t = AY£ + AKP/ A|< 2# ( 8o*R + 4.o — 16 0 *p/AK2 ) *7K 
U = AYE/Sq r T(DEm ) / ( X 2+ ( R— 1 a )*{R.i o 1 ) 

72 IF { TP~Th (niT) ) 75 s74>73 

73 NT“NT+1 


GO TO 72 

74 G=GA(NT) 

GO TO 76 


75 

76 
71 


G-GMnTbj ) + (Tp-TH( niT-2 ) )*(GA(nTJ-GA(nT-i) }/ (Th(nT ) ' T KwT.i ) ) 
D I (J) =G*U a N 1 ' J 

TP=TP +DTP 
K = 2 


Ul -U I + H* ( D I ( 1 ) +4o*DI ( 2 ) +D I ( 3 ) ) 

D I < 1 } =D I { 3 ) 

IF ( T P-P i ) 77 » 77 > 78 


78 

UGA( JR ) = 


19=000 


U I = 0 o o 

97 

CONTINUE 

999 

RETURN 


end 



n n nn 
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$ I 5FTC SU1 5 

c ^BrcUTimE GAmCYl ( CO ,Xc » R o , U G » I r »Xp R E5» ,U G p » RR p , R A } 

, C - SUbRoUTI NE p To CALCULATc The AXIA L VELOCITIES I nDUCEo £ { 
-^TrEnGTh Vq-rEX cylinders 
c " 

CO'Sn'Sz RPI25I,UG ‘ 25) ’ X ^ ES ‘25).UGP(25.25).RA(25) 

c 


PI=3o1A15926 

NR=N 

CR=CD*2o 

IF(IR)2»1»2 

1 Xg= ( xc-l o ) *CR 
C 

c compute the velocities induced at the ppopEllEr statiqnc 

20 DO 30 J=1»NR 
R=RP ( J ) 

IF(r)22»21»22 

21 U = SqrT ( ( XG*XG ) A L ) 

U=XG/U+1 0 o 
U = 0 o 5#U 
GO TO 30 

22 R?SiT 8 Si;?., R+1>) 

RM 1= ( R— 1 0 )*CR-l a ) 

SN5“SQRT ( XG2+RM1 ) 

SnB=XG/Smb 
SkG2=4 0 *r/ (XG 2oRPl ) 

CALL ELL I PS (SKG2»2K»Ze) 

SKG=S0RT (Skg2) 

CALL Arcs I N ( Si<G 9 askg »G ) 

CALL ARCSIN(SNB,3ETA,0) 

NDEX«0 

IF (BETA) 23 * 24^24 

23 BETA="BETA 
NDEX^l 

24 CALL LAMBDA ( ASKG sBETAjHLMB) 

IF(NDEX)25,26,25 

25 HLMB=**HLMB 
NDEX“0 

26 CONTINUE 
U“ 0 25*HLMB+ o 5 
TM=SKG*Xq*ZK/ ( 4 0 *PI ) 
f.M*TM/.SQRT (R) 

U=U+TM 

30 UG(J)=U 

GO TO 99 

2 NCYL-N 

Compute the velocity induced • at the rEFrence cylInderby 
the trailing vortex cylinders 


Constant 


DO 60 J=1,NCYL 



31 


42 

142 


242 


RGR=RA( J)/RRP 
RGR2*(RGr + i,.J* (RGR+lo ) 
RGRM = ( RGR — i o ) -:{• ( RGR-1 0 ) 
rTrGR=GqrT (rGr ) 

IF ( J-NCYL ) 31 , 32 ,32 

XP=l o 0 

P.GR=lo® 

R6R2=4oC 
PGR M = 0 o 0 
GO TO 33 
X P * X c 

DO 60 JR= 1 9 I R 
X=XpRES( jr) 

XG=X— XP 

xgr=xg*cr 
XGR2=XGR*XGP 
IF(XGR) 38*42*38 
I F ( RGR-1 6 ) 142 ,242 , 34 2 
D G ( J > = 0 ,0 
GO TO 60 
UG ( J } =0 o 2 5 


GO TO 6 A 

342 UG(J)»0o50 

GO TO 60 

33 CONTINUE 

G«Sq = 4 o *RGR / ( XGR2 + RGR2) 

GK=SGRT(GKSQ) 

CALL ELLIPS(G K Sq,Zk*2E) 

NDEX=0 

SNB»XGR/SGRT ( XGR2+RGRM ) 

CALL A R L S I N ( S m B » 3 E T A » „ ) 

CALL A R C S I N ( G K » A S K ? 0 ) 

IF(B ETA) 34,35,35 
3A BETA “-BETA 
NDEX=— 1 

3 - CALL LAMBDA! A$K, BETA, HLMB) 

I F(ND EX >36,37,37 

3 6 HL.MB®— HLMB 

BETA “-BETA 
NDEX=0 

37 .. CONTINUE 

I F C XGR > 41 , 42 ,43 

43 NDEX«-1 

HLMB“«*HLMB 

XGR=-XGR 

41 IF(PGR-1„ ) 141,241,241 

141 uG ( J ) = GK*XGR*ZK/ ( 4c *p I *RTr6r ) — *>.25*HLmB 

GO TO 45 

241 UG ( j )wGk*X'(5r*ZK/ ( 4o^pI ^RTrGr ) +0 ■© 5 O+ o 7 5*Ht MB 

45 I F ( NDEX ) 46 3 60 , 60 





46 fiSiS?8 GR 

47 IF (RGR-1. j 147,247,347 

1 47 UG(J)=-UG(J) 

GO TO 60 

247 UG ( J ) “-UG { J ) +0 o 5 

GO To 60 

347 UG ( J ) =~UG ( J ) + 1 0 0 

60 UGP ( J » JR ) =UG ( J ) 

99 - RETURN 
END 



--4 -J 
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$ I BFTC SU 1 6 

SUBROUTINE bncoef { c »bg I ) 

o Fourier analysis of radial velocity induced by an- 

C ACTUATOR DISK IN THE EXIT PLAN'S — >Q ( N ) 

DIMENSION 66(6*2) »BGI (6) *BG0<6> 

COMMON MZZZ 
PI =3 0 1A15927 
MZZZ=— 10 
N = 1 0 0 
KMAX=6 

DO 10 K-- 1 , KM A X 
In BGI (K)=U oa 

EPS »PI/I80.e 
PIM =P I —EPS 
C N = N 

DELTH = o 5*P I M/CN 
NP-N+ 1 
DO 6 1=1 > No 
DO 3 J = 1 * 2 
F I =1—2 
FJ = J 

Th=FI*PIM/CNaFJ*DELTH 
STH=S I N ( TH ) 

CTH*C0S(TH) 

XG=-0 o 5#C* (lo+CTH) 

XGSQ=XG*XG 

GK s Q= x GSq/ (Ao+XGSq) 

GK-SQRT (4o / (A Oi -,XGS0) ) 

CALL ELLIPS (GKSq^ZkG.ZCG) 

6GCao5/Pl*( ( Gk-2 o/Gk ) * ^kG+2 »/Gk*ZCG ) 

STHK“0o 
CTHK= 1 o 
DO 2 K = 1 » KMAX 

bg (k * j >=bgc*cthk 
con=sthk*cth+cthk*sth 

CTHK=CTHK*CTH-STHK*STH 
STHK=CON 
CONTINUE 
DO 5 K=1»KMAX 
I F { I<~1 ) 4 , 5 » 4 

A B G I (K) = 8GI (K}+DElTH*(BG( k »1 ) + ( BGO ( K ) +BG { K » 1 I +BG( K»2> )/3 « } 

5 BGO I K ) = BG ( k » 2 ) 

3 CONTINUE 

EXTRA =»»€pS/pI*AL0G{4 o /EPS*S0PT(2c/C) ) 

SI— 1 

DO 7 K* 1 > KMAX 
SI —SI 

7 BGI (K)=2o*(BGl (K >+SI*EXTpA) /ol 

BGI(l)«o5*BGI(l» 

MZZZ«1 

RETURN 

END 
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JIBFTC_SU17 

sudq 0 u T 1 me a n co£f (ngyl.n »ra:» x p » cd , RRP , MP u N , K PR I M T DS ,a . -• •; j 

DIMENSION F ( 50 ) s BS ( 6 ) » A ( 25 ? 6 ) sAS ( 25 9 5 ) ,ra ( 2 r ) 

COMMON M2ZZ ’ ' 

142 FORMAT ( 1H1 s4X ,*RUN NUMBER* , j 5 , 49X »*PAGE* ? I 3 9 2H-F > // ) 

242 Format ( iox»*foUrier cos i me series coefficients — output index 

113) 

NPAGE =0 
N0UT =6 
.2 0 CR=CD#2 o 0 

PI“3o1415R27 

21 CN®N 
DTH=PI/ ( CN-lo ) 

TH-H^TH 

C 

c compute b**star(m) fouriep coefficients 

22 DO 30 J*1,N 
TH»TH+DTH 

CSTh=CoS ( TH ) 

XG«*—o 5*CR# ( 1 0 n CSTH ) 

GKSq=XG*X6 

G K S Q = 4 © 0 4- G K S Q ~ 

GKSQ=4oO/GK.So 

CALL ELLI PS(GKSQ s ZK 3 ZE) 

GK*SqrT(gkSQ) 

30 F { J ) «0o25+ (GK*XG/4 0 *ZK/PI ) 

I F ( NPR I NT ) 171,171,170 

170 NPAGE =NPAGE n l 
PRINT 142 ,NRUN, NPAGE 
PR I NT 242 , NPR j NT 

171 CALL FqUrCS (F >F ?N*NO^T »NPRlNT) 

DO 31 J=l»NOUT 

31 BS < J ) = F ( J ) 

C COMPUTE A ( M ) FOURIER COEFFICIENTS 

r 

x- 

32 DO 40 J=1 , MCYL 
RRA*RRP/RA t J) 

RRAP**RRA + 1 0 0 
RRAp*RRAp«RpAp 
RRA4«RRA*4 o 0 
RASQT=SQRT(RRA) 

TH«-DTH 

33 DO 38 K*1,N 
TH»TH+DTH 
CSTH=COS { TH ) 

XR=CR*(lo— CSTH)*-o5 
XG=XR i »(XP*CR) 

XG»XG*RRA 
XGSQ =Xq*XG 
GKSQ“RRAP a XGSQ 
GKSQ b RRA4/gKsQ 

GK s SQR.T (GKSQ). 

CALL ELLIPS (GKSQ»ZK»ZE)- 
F(K)*{ (ZE-ZK) /GK+(GK*ZK/2o ) ) /PI/RASQT 



■06 

33 CONTINUE 

I F ( NPR I NT ) 173,173,172 

172 NPA6E=NPAGE+1 
PRINT 142 ,N RUN ».N PAGE 

PRINT 242»NPRjNT 

173 CALL FOURCS ( F , F , N , NOUT » NPP I NT ) 

DO 39 K a l , MOU t 

39 A ( J 9 K ) = F ( K ) 

40 CONTINUE 
C 

c compute a-starin) Fourier coefficients 

42 DO 50 J* 1 9 nCYl 
RRA*rrp/ra(j) 

RPRIM= 1,0/RPA 
RRAp = RRA+l 0 0 
RRAP=RRAP*RRAP 
. RRA4=4 0 0 *Rra 

RASqTbSqrT (rra ) 

th^-dth 

43 DO 48 K = 1 9 N 

TH=TH+DTH 

CSTh-COS(TH) 

XR = CR*( 1oCKSTH]*o5 

XG = XR“( Xp*£P ) 

XG=XG*RRA 
XGSQ=XG*XG 
GKSQ=RRAP+XGSo 
GKSQ«RRA 4/GKSo 
6 K<=SQRT (GKSO) 

CALL tLLl PS (GicSq »2k'»ZE) 

NDEX-0 
SNB=RR.A^>1 o 0 
SNB«SNB*SNB 
SNB*XGSq+SNB 
SNB= 1 oO/SORT ( SNB ) 

SNB»XG*SNP 

CALL ARCS IN (SnBsBETAsO) 

CALL ARCS IN (G< 9 ASK »0 ) 

XG“XG*RPR IM 
IF (BET. A) 62, 63, 63 

62 BETA»~BETA 
NDEX=“1 

63 CALL LAMBDA ( a.Sk » BETA *HLMB ) 

IF (NDEX) 64,65 ,65 

64 HLMB“-HLMB 
N D E X = 0 

BETA«*BETA 

65 CONTINUE . 

IF (XG) 7 1,72 9 73 

72 IF(RPRIM«1 o 0) 74,75,76 . 

74 F ( K ) =0 o 0 

GO TO 48 
F(K)»-0 o 25 
GO TO 48 


75 
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76 F ( K ) -“0 0 50 
GO TO 48 

73 NDEX=*-*1 

HLMB=“HLMB 

XGa-XG 

71 IF ( RPRINM-1 „0) 77 » 78 » 78 

77 F ( K ) 5 * j GK^Xg*Zk/4 c 07pi*rAS0T 

F ( K ) =F ( K ) ^ ( o 25*HLMB J 
GO TO 79 

78 F ( K ) =w GK*Xg*Zi</ 4 j 0/P I *RASqT 
7 ( •< ) =F (K )*“» c 5 — ( o 25*HLMB ) 

79 I F ( N D E X ) 8 0 , 4 8 , 48 

Or XG~*"XG 

NDEX»0 

I F ( R PR I 0 )3l»82983 

81 F ( K1 s**-F ( K ) 

GO TO 48 

82 F ( K ) = “F ( K ) “0 „ 5 
GO TO 48 

83 F(K)—F(K)-1oO 

48 f(K)=wf(k) 

IF(NPRINT)175, 175,174 

174 NPAGE =NRAGE+1 

PRINT 14 2 , MrUnj , NPAGE 
PR I NT 242 ,NPRTNT 

175 CALL FoUrCS ( F , F , N , NqUT , NpR 1 1\! T ) 
DO 49 |< = 1 sNOUT 

49 AS ( J , K ) =F ( K ) 

50 CONTINUE 
RETURN 
END 
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$ I B F T C SU 18 


subroutine cncoef.(gv»xcj 

OI.vEnSIon A ( 6 ) » P ( 6 ) 9 AS ( 6 ) , bS ( 5 ) sCni ( 6 ) 9 Cn 2 < 6 ) »o ( 
DIMENSION CVG(6) 

'-OMMON M 2 AZ » CD » RU »r 1 »R 2 * R 3 9 P I »B»dS»A»As»P 
VG = 1 oO/GV ' 


6^7) ».>(6»6l 9 XC ( 6 ) 


C = 2 0 *CD 


c CALCULATE COEFFICIENT matrix 

c 

TR„* 2 „*R* 

35 C 4 D«C/ 8 . A 

36 C LN«ALOG ( 32 o 0 /C ) — 1 o 0 

Q t 1 9 1 ) = 1 o 0 *-p ( 1 s 1 )-Tr j*C 4 D*CLN-R 1 *C 4 D 
Q ( 1 9 2 ) = R 2 / 4 ° 0 *C 4 D— TRU-K-CA-D-sf-U o 5 -K-CLN 
Q ( 1 s 3 ) =~P ( 193 } A Rl* 0 o 5 *C 4 D+R 3 *C 4 D/ 6 o 0 
Q( 1 5 4 )=-R 2 *C 4 D/ 4 oO 
q ( 1 , 5 ) =-p ( 1 9 5 ) — R 3 *C 40/ 6 o 0 
Q ( 1 » 6 ) = 0 0 0 

'Q( 2 » 1 )=P( 2 » 1 )*‘*Tr 0 *C 4 D ~2 « #R 1 *C 4 D*ClN-P 2 *C 4 D 
G( 2 » 2 )=P< 2 »?)«-lo 0 -R 1 *C‘ 4 D-k-Clm a (p 1 a R 3 ) *CiD /4 o 0 
Q( 2 » 3 )=C 4 D*(-P 0 -R 2 / 3 c ) 

Q ( 2 94 ) =P (2 94 ) — C 4 D*( R 3 a r 1 +R 1 5 /pcO 
0 ( 2 , 5 ) =“R 2 *CAd/ 6 <; 0 
Q( 2 , 6 )=P ( 2 , 6 )-R 3 *C 4 D/ 8 o 0 

Q ( 3 9 1 ) =P ( 3 »1 )-.C 4 D*( R 1 +R 3 )- 2 o*'P 2 *C 4 D*CLN 
0 ( 3 , 2 ) =R 0 *C 4 D/ 2 * 0 — R 2 *C 4 D#CLN 

Q ( 3 9 3 ) =P (3 » 3 )-l* 0 -C 4 D*(Rl/ 3 o+R 3 / 2 « ) 

0 < 3 94 ) =C 4 D^ (R 2 / 8 ^-R 0 / 2 o) 

Q ( 3 9 5 ) = P { 3 9 5 ) + C 4 D-:;-( R 3/ 10 o-Ri / 6 o ) 
Q( 3 , 6 )=-R 2 *C 4 n/ 8 o 0 
Q( 4 ,l)=P( 49 l )«C 4 D* ( R 2 + 2 °*R 3 *C|_N') 
Q( 4 , 2 )=P( 4 9 2 ) + C 4 D* ( R 1 /4 p-R 3 #CLN > 

Q ( 4 9 3 ) = C A D* ( P a / 3 o —R 2 / 2 0 ) 

Q ( 4 9 4 ) = P { 4 9 4 ) cj 1 „ 0 — pl*CAD/So A 
G( 4 , 5 )=C 4 D-k-(R 2 / 10 o“R 0 / 3 o ) 

0 ( 4 , 6 ) =P ( 4,6 )— R 1 *C 4 D/ 8 o 0 
Q( 59 l)=P( 5 s l )x»R 3 *C 4 D 
Q ( 5 , 2 ) = R 2 # C 4 D/ 4 .„ C 
Q( 593 )=P( 5.9 3 )+C 4 D*(R. 1 / 6 o“ , R 3 / 2 0 ) 

Q ( 5 9 4 ) = L4D* ( p A / 4 c -R 2/ 4 < ) 

Q(5»5)=p( 5»5) -1 ,0-C40*Rl/15 c 
Q( 596 )=-R 0 *C 4 d/ 4 o 0 
0 < 6 9 1 ) =P ( 6 9 1 ) 

Q( 6 , 2 )=P( 6 » 2 )+R 3 *C 4 D/A i 

Q ( 6 9 3 } = C 4 D *R 2 / 6 - u 0 
q( 694 )=P( 694 )+ C 4 D*(pl/ 5 o A R 3 / 4 o ) 

Q ( 6 , 5 ) = (_4D* ( RO/5 „**R2/6 <, ) 

Q ( 6 96 ) =P { 6 9 6 ) -1 o 0 -C 4 D#Rl /8 * 

C 

c calculate constants 

c 

DO 40 1 * 1,6 

CNl ( I )=T rO*6S( I )-(2 0 *6( I ) ) 

GO TO ( 41 , 42943944 , 45 , 46 ) ,i 
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41 

42 

43 

44 

45 

46 
40 


51 

52 

53 

54 


55 


56 

50 

C 


6 1 


GO N }o I 46 CNl ( ^ " ( P1* bM 2 ) ) + ( 02*85 ( 3 ) ) + ( j? 3 *BS( 4 ) ) 

1I5 ) ) 1 J=( ~ N1 ( 1 )+R1 * (dS( 3 ) + ( 2 -*t>S( 1 J J )+R 2 *.<tJS ( 2 > +6S { 4) ) +03*(BS (3 )+f 
GO TO 40 

l^n 1 ,=< " N1( 1 )+ r 1 *<BS ( 2 )+SS ( 4 ). )■ +R2* ( BS ( 5 ) + ( 2 O *BS { 1 }• ) ) + ?3 * ( SS { 2 ) + PS 
GO TO 4C 

C,\!l( I)=L.mi{ I J +p 1 *.( BS { 3 } + 3S ( 5 ) )+R2*(BS(p )+BS( 6 ) ) + 2 «, *p 2 *BS ( l ) 

GO 1 0 4 a 1 

Cm 1 ( I ) s=(_f>) l ( I ) n 0 1* ( 6S ( 4 ) +BS { 6 ) ) +R2*8S ( 3 ) +P3*^S ( 2 ) 

GO TO 40 ' " ‘ 

CN 1 ( I ) = C f.i 1 ( I ) + p I*p5. ( 5 ) +R 2 *3S (a.)+r 3 * 8 S( 3 ) 

CONTINUE 
A S { 1 )* s AS(l)+l o 0 
DO 50 1=1,6 

C,'\l2( I ) = — 2 o * ( A ( I ) ) + ( TOU* { Ac ( I ) j ) 

GO TO ( 51 * 52 * 53 9 54 9 55 9 56 ) * I 
C N 2 (I >=CN 2 < I )+Rl*( 46 ( 2 ) )+R 2 *(AS( 3 ) )+(?•**{ AS ( 4 ) ) 

GO TO 50 

C,M2 ( I J = C N 2 ( I ) + E1* (2-* i AS ( 1 ) ) ^ AS ( 3 ) ) +R2*< AS ( p ) +AS (4 ) ) + * 3 * t A 5 £ 3 > +AS ( 
15)) 

GO TO 50 

CN2 ( I ) = CN2 ( I ) /nR 1"- ( AS ( 2 ) + AS ( 4 ) ) +R2* ( 2 o * ( AS ( 1 ) ) +AS ( 5 ) ) +R3* ( AS ( 2 ) f-AS 
1 ( 6)) 

GO TO 50 

Cn2( I ) =*Cm 2 ( I ) +R 1* ( AS ( 3 )+aS (5 ) )+R2*(AS(2)+AS(6^ ^ +R3*2 ° * ( AS ( 1 ) ) 

GO TO 50 

C N 2 ( I ) = C M 2 ( I ) +R 1 *(AS( 4)+AS(6 ) } + R2* ( AS ( 3 ) ) +p 3 * ( AS ( 2 ) ) 

GO TO 50 


Cm 2 ( I ) - CM2 ( I ) +R 1* ( AS ( 5 ) ) +R 2 * ( AS (4 ) ) +93* { AS ( -a ) ) 

CONTINUE 

AS ( 1 )*AS(l)-l o 0 


CALL MAT R I X ( Q » 6 ) 

DO 60 1=1,6 
XC( I)=0 OA 

CVG(I)=CN1(I)+CN2(I)/GV 
DO 61 1=1,6 
DO 61 J = 1 > 6 

XC( I ) = XC ( I )+Q( I ,J)*CVG(J) 

RETURN 

EMU 
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SIBFTC 


SUBROUTINE PRESS 


dimension 

D I M E m S I q m 

dimension 

0 I m E n S I o N 
D I M E N S I o ,\i 
DIMENSION 


UGp ( 25 


TCRAT ( 7 ) >X v ( 19 ) 9 OVA { 8 » 7 ) s F ( 8 , 19 ) 9 UP ( 25 ) »UM { ? 5 ) »SUMU ( 25 ) 
o I t 5 ) sUEdp (25 ) 9 VCdp ( 25 } sCsT (5 ) ,SAVU (75 ) 9 FX ( _> 5 } 

B (6 ) 9 BS ( 6 ) 9 Sa ( 6 ) » SAS ( 6 ) ,p (695) 

C($) 9A ( 25jf ) *^SC 2 5»6) » D ( 5 ) 9 Ds ( 6 ) » SC ( 6 ) ,GS( 6 ) »H(M 
Ro 25)9 o R (25 ) ^ t A(25 ) S I c B L D( 25) » pA ( 25 ) » XppEo ( 2 5 ) 9 pHl ( 5 ) 
Uu ( 2 ? ) »UQ0 ( 25 ) 9 USD ( ?5 ) 9 UV ( 25 ) »(JCB ( 25 ) » VCb ( 25 ) » 


dimension alpha (25) 


5*25 ) 1 * * * * * * * 9 UGA ( 25 ) »GV< 2 5 ) 9 


^RV( 25) 9 ' 


- -- • .Cpp ( 5 9 25 ) 9 CPM ( 5 9 25 ) 

' ALL ( 25 ) ? J S T |_ ( 2 5 ) 9 TAlk(2'D) 


COMMON MZZ2 9 C 0 9 RU 9 RI »R2»R3»pI ,B »bS» SA» SAS vp 

SJ . A Pl/NR u Nf NPLD * N z » «Z , N pp £S * I N T I M E » jvjUpp , NP A i , N pH I »nppInT 
COMMON^ NtAp 2 / C 9 A 9 AS 9 D 9 DS '9 SC 9 GS 9 H 

c 0MM0N/NEar3/prp ,Xo ,Z ,BLO, ob,bp >e>TA »TCBLi>* TC , R Cb R p,M P A, A L F ,XppES » 

1 RAjX^BsXR&cLCBC 5 RMAX 9 pH I »COR J ? CQRCB 

Common /neara/cg 9 00 d 9 jgd s uv 9 ucb »\/cb»ugp »uga , gv» gam»gr v,cpp »cpm 
■COmMON/iMEARS/aRJ 9 ARJPsEES 9 Rad 9 cl 9 alpha 9 stall 9 jstlVt ’ LK 

DATA TCr AT / Oo 06 9 0 w 08 & 0 o 12 ? Oo 15 9 0 o 10 9 0<^ 2 1 9 ^ o 2 A/ 

UATA XV/O 0 O 9 D 0 OO 5 »0»O1 2 5-*u«.U 25 , »0«C5 »0«075»0^ I0*0o 15 »0*2 »0. ?5 f 0»3 , 

1 _ ^ u “ A 9 0 = 5 9 U o 6 9 0 o 7 9 0 - 8 9 0 . 9 9 C o 9 5 3 1 0 0 / 

UATa DV A ( 1 9 1 ) 9 OVA ( 1 9 2 > 9 DVA ( 193 ) 9 OVA (3 94 ) 9 OVA ( 595 ) 9 OVA ( 795 ) 9 OVA ( i ,7 

1 xT’ 9 ? 2 ,2 1,015 s ' 1 '' ,364,0o 9 8490 ‘'696 9°o462 9 0o478/ 

, / uVrt ( 2 9 1 ) 9 OVA 1 2 9 2 ) >^a( 2 n) 9 OVA ( 2 9 4 ) S DVA( 2 96) 9 0V7( 2 9A) 9 OVA ( 0,7 

1 ) / 2 . 90 , lo 7 95 , lo 3 l , 0 w 97 l* 0 o 694 * 0 o 561 .' 0 . 478 / ' " ' 

iT/? °^ ( 3 ? }A L)VA( 3’2> 9 DVA( 3 , 3 ) 90 VA( 3 94 ) 90 VA( 3 , 5 ) ,UVA( 3 , 6 ) , OVA ( 3 , 7 

1 ) /lo9b8»lo475»-lol99»0 c 93A*Oo685»0o558*u»A 7 8/ 

UA A UVm ( 497 ) ,UVrt( 49 2 ) *WVA( 4 * 3 ) ,OVm( 4 , 4 ) 90 VA( 4 >5 ) 90 VA( 4 » 6 ) , OVA < 4,7 

1 ) /I 06U 9 1 0 312 9 1 0 112 9O0 9O 9O0675 9O o 557 9O 0 478 / 
umTa 0 V A ( 5 9 x ) 9 OVA ( 5 , 2 .) 9 UVA( 5 9 3 ) »DVA( 5 , 4 ) 9 Dv'A( 5»5 ) ,DVA( 5 , 6 ) , OVA ( 5, 7 

1) /I o342 9 1 o 178 9 1 6 028 »0 0 861 »0 o 662 »0o555 90 o 478/ 

DV * ( 6 9 ] ) 9 OVA (6 9 2 ) 9 OVA ( 693 ) 9 OVA ( 694) 9 DvA ( 6 , 5 ) ,DvA ( 5 , 5 ) ,DVA ( 5 ,7 
; 167 9 I 0 O 65 »o 0 9469U 0 8I8 »0o648»0o550»d o 478/ 

OV A ( 7 9 3 ) 9 OVA ( 7 > 2 ) 9 OVA ( 7 9 3 ) 9 OVA { 7 , 4 ) ,L)VA ( 7 ., 5 ) ,DVA( 7 , 6 ) 9 OVA {7 , 7 
>050 9nc 96 4 90^870 9Q o 77l90o63290o542j0 o 478/ 

F( 191 ) 9 FJL 9 2) »F(1, 3 ) »F( 1 9 4 ) 9F(i» 5 j ,F ( 5 , 6 ) 9F(i » 7 ) ,F( 1 » 3 ) »F{ l9 Q 

1 ~ ^ 9 1 p 3 ^^ , ll^ ,r p ? l2)9F(X9l3)9F(i9]_4)»F(i9l5)9F(i9i5),r(i,i7)9 

2 r ( 1»18 9 ( 1 9 1 9 _) /0 » 9 O " 9 38 9 1 « O 57 9 1 o 089 9 1 " 1 O 3 9 i * 1 07 9 1 o 10 1 9 1 o 098 » 

1 ^ oU Vl9lo0869ioU789i«0669lou539l-0429lo0289loUi390«.990 90 u 97490o/ 

uaTA F (29l)9F( 2 92)9F(293)9F( 2 94}’f :: ^29559 F( ? 96)9F(297)9F(293 )>F(299 
l_)_5r(29iU)9F(29li)9F(29 l 2 )9F(29l3 ) 9F(29 1 4)9F(2»15)9F(29i6)9F(7, 17 ), 
2 F (29l8)»F(2»l9) /0o«0o890»lo O 5 U 9lolU5 5 i o i289l=l339iol 3 O9l0l289 
1 1« 122, lo H4, lo 106 ,l°0899lo072,lc0549 1=0399 1=017,0° 9 34,0= 969,0°/ 

DATA F( 3 ,i ) , h (3 , 2 ) , F ( 393 ) 9 F( 3,4) 9F( 3 »5) 9F( 3 , 6 ) ,F(^, 7 ) »F( 3, 3 ) ,F(3, 9 
l) 9 M 391 O) 9F(3,h) »F( 3 » 1 2) » f ( 3 » 1 3 ) *F< 3 , j 4 ) *^( 3 , 15 ) ,F( 3>16 ) ,F( 3>17 ) , 
2 F ( 3 9 1 8 ) 9 F ( 3 » 1 9 } /u,> 9 V 0 8 n 0 9 lo ^^5 9 lo 1 149 1 o 174> 1 o 184»1 o 1 38 9 1 o 188» 

1 lo 183 9 lo 174,1c 162 9 1 = 135 9 l' U l „8 9'lo080»l o053'.9lo0.2290 »978 9 rt « 952 »Oo/ 

OaTA F ( 4 9 ! ) 9 F ( 4 » 2 ) *F (4 9 3 ) »'F,t 4,4 ) ->F'( 4 V 5 } »F I 4 9 6 ) 9 F ( 49 7 ) »F ( 4 , g ) »F ( 4,9 
l)_ 9 F ( 49 l 0 ) » F ( 4 , 1 1 J 9 F ( 49 i 2 ) 9 F (4»13^ 9 F '( 49 i 4 > 9 F ( 4 9 1 5 ) 9 F ( 4 ,i 5 ) ,F ( 4 , 17 ) , 
2F(49l8) 9F(4 9 3,9) / 0o90o7399Uc'9669l o li29lo2049l u 2249 A = 3339lo2339 

1 1 o 229 9 1 o 218 9 1 « 204 9 i A 1 70 » 1 1> 13 T 9 1 c098 »1 j 064» 1 o 024 9 u o 9 , ? 2 9 0 o 934 9 O 0 / 

DATA F ( 5 , 1 ) ,F ( 5 9 2 ) 9 F ( 5 9 3 ) 9 F ( 5 , 4 ) S F ( 5 , 5 ) 9 F ( 5 9 g ) 9 F { 5 9 7 3 9 F ( 5 9 3 ) 9 F ( 5 9 9 


0 aT A 
1) /l, 
DAT A 

1) /I 
D A T A 
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]) >. r" ( 5 9 1 G ) , F ( 5 s ]. l ) , P C 5 1 2 ) > F ( 5 , i 3 J 9 F ( 3 , ] 4 } » F ( 5 » i_ 5 ) 9 F ( 5 > ]_ 5 ) 9 '“( 5 , 17 ) 9 
2 F ( 59 18) * F ( 5 » 19 ) /Ou »ut 682 »u*926 » lolu3» 1 o 228>1»'264» lo 276* l»278-» 
i lo2^59lo262 ! >.i , o2‘+7».i<.2G59i t l54»lolI6»loU74slt.025»Uo966*0*914*0o/ 
GaTA F ( 5 9 j ) 9 F ( 592 ) 9 F ( q s 3 ^ > F ( ft s 4 ) * F { £ , 5 ) 9 F ( g , g 5 9 F ( g 9 7 5 > F ( 5 , 3 5 9 F ( g 9 q 
1 ) 9 F { 5 9 1 w ) 9 F ( 6 » 1 1 ) 9 F (5s 12) f F ( 5 9 1 3 ) 9 F { 5 9 t 4 } s F ( 5 9 n 5 ) 9 F { g , 3 g ) 9 F ( § 9 i 7 5 9 
2 fl (69lS)9F(6,i9) /0v9JF,630,0o8879l<'0879 ]'o.242»l‘'29T9l«2179l®3259 

1 1 u 320 9 1 c 3Q6 9 1 “ 290 9 1 " 240 9 ~ 17S 9 X ^ 1 33 9 1 “ 085 9 1- 0 Q27 9 0 * 9 37 9 0 “ 3 95 9 O u / 

OaTa i~ ( 7 9 1 ) s F ( 7 9 2 ^ 9 F < 7 9 3 ) » F ( 7 9 4 ) 9 F ( 7 9 5 ) 9 f ( 7 9 6 ) 9 F ( 7 » 7 ) »F(798>9F(799 
l) >f ( 7 9 1 U ) 9 F ( 7 9 1 1 ) 9 F ( 7 9 X 2 ) 9 F ( 7 , 1 3 5 4 F ( 7 , } 4 ) 98(7,15) 9 F t 7 9 1 ) » F ( 7 9 1 7 ) » 

2 F ( 7 9 18 ) ,F ( 7 •> 19 ) / 0 o s u 5 7 9 9 u i. 8 48 9 1 o 0 6 3 9 1 o 244 9 1 * 322 » 1 o 354 » 1 » 374 » 

1 c- 3 6 b 9 i o o 5 p 9 i - b .2 3 9 i»277 9 i » 20F 9 i v 13 1 9 1 0 0 9 7 9 loU329be 9*t4 » 0 08 7 9 » 0 <• / 

C 

700 format (///lux 9 *.duct ppcssurc di 8 trU 6 ut ion • calcut ion assumes t/c 

1 ** >F5 =3//) 

r 

V. 

I F ( MZZ 2 ) 2 9 1 9 2 

1 CONTINUE 
C 

CSET Up A TAbLE qF CONTINUOUS VELOClTU CORRECTION FACTqrS 
C 

DO 31 J = l,7 
TCP«TCRAT( J) 

. J T C* J 

I F ( T CW C R A T ( J ) 5 32 933 931 

32 IF(J-l) 132,132934 
31 CONTINUE 

132 PRINT 7 0 n 5 TCP 

33 DO 38 N = 1 , 19 

IF (N“7 )37,37938 

37 D V A { 8 9 N 5 = D V A ( J T C 9 N 5 

38 F ( 8 , N 5 = F { J T C 9 N ) 

GO TO 39 

34 DElT C= ( TCrAT ( j TC )— TCR AT ( jTC— 1 ) )/ l TC rAT ( jTC ) — TC ) 

J = JTC 

DO 35 N = 1 9 1 9 

3 5 F ( 8 9 N ) = l F ( Jj»1 » M > ®*F ( J 9 N ) 5 /DEl_T C+F ( j 9 N ) 

C 

C SET Up A TAbLE qF THE i> I SCoN’T I wUrjUS VElqCITY Cqrp-CTIpn FAACTqrX 

C 

DO 36 N= 1 , 7 

35 DV A ( 3 > N ) = ( OVA ( j— ] 9 K) ) -OVA ( J ,!\j ) ) /DElTC+DVA ( J 9 M 5 

39 CONTINUE 
RETURN 

C 

2 CONTINUE 

c look up Continuous velocity correction factors 

c 

DO 49 N=1 9 IR 
DO 42 0*1,19 
JK-*J 

IF(XPRES(N)-XV( J) 541947,42 
42 CONTINUE 

41 DELX=(XV( jjO-xyt JK-15 )/ { XVt J|<)-XppESIn) ) 

FX ( N ) = ( F ( 8 9 JK- 1 ) -F ( 8 9 JK ) ) /DEtX+F ( 8 9 JK ) 
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GO TO 49 

47 F X ( N ) = F ( 8 ? J < ) 

49 CONTINUE' 

SmALF=SIm ( AI.F/RAD) 

csalf=cos ( ALF/RAD ) 

compute: velocity i n ouced by trai l i n g vqp.tex cylinders 

call GAMCY L ( CD »Xp »NZ » RB »«JG » I'R pXprES »UGd ,ddd,rA) 

DO 1 0 J = 1 3 I p 

* SUMU(j)=loO 

DO 12 JalsIR 
Do 12 N = 1 9 N 2 

2 9u!/|U ( J ) =SUmU ( j ) +UGP ( W » J ) *GV ( ?< ) 

compute velocity Induced by the duct BqUmC vorticity 
IPC — I R 

'-All VTXrnG (C'DsXp* IRC»Rb»C»UGD'»XpRES 9 p ) 

DO 14 J= 1 5 I R 

4 SUMU { J ) =SUMU ( J ) +UGD ( J ) *GV ( NZ ) 

COMPUTc VEloCITY INDUCED bY THE CEnTrE dqDY 
IRC — IR 

CALL HUB ( CDs XR 9 XCB s LLCuC sRMAX> IRC j XPRES ,RB ?RRP sUCBp » VCBP ). 
DO 16 J=1,IR 

6u;'4U ( J IcSUmU { j ) +UCBo ( j ) #CorCB 
SUMU { J ) =SUmU ( J )*CSALF 
6 SAVU(J)=SUMU( J) 

I F ( SiMALF) 20, 2 1 s 2 0 

Compute velocity induced by alpha vortex rings 

20 CALL A |_ F R N G ( CDs iRs SCs UGA» XPRES) 

GO TO 23 

21 DO 22 J= 1 , Ip 
2 2 UGA ( J ) = 0,. 0 

23 DO 99 M= Is NdHI 

CSPH I = COS (PHI (M) /RAD) 

SNPHI= SIN (PHI (M) /RAD) 

DO 90 J= 1 , Ip 
90 SUMU ( J ) = SAVE { J ) 

DO 25 J= 1, IP 

25 SUMU(J)= SUmU( j) + 'UGA l j)* SnAlF*CSpHI 
30 CONTINUE 

Correct continuous induced velocity 

DO 48 N= 1 9 IR 
48 SUMU ( N ) = SUmU(m)* FX(N) 

Compute discontinuous portion of surface x velocity 

c 
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SIBFTC SU20 

SUBROUTINE OUTPUT 

DIMENSION U(6)3 bS(6)» SA(6)? 6AS(6)9 p(6? 6) 

DIMENSION. C (6) » A(25» '6 ) » Ao( ? 5 9 6 j 5 D( 6 ), DSIfe)? SC ( 6 ) 9 GS( 6)» H( 
1 .- 6 ) 

DIMENSION rc ( 05 ) » Bp ( 25 ) » BTa<25 ) » TCd^DIp?)? rA( 25 ) » XppES(25)» 

1 PHI ( 5 ) 

DIMENSION UG( 2 5). UqD£ 2 5), UGO<25)» UV< 25 )» UCB(? 5 ) 9 /CB<25> » 

, v_ ' up 25 a 25)s> UbA<25) 5 G v (25)* GrV( 25)» Cpp(59 25 )» CF'j(59 25) 
l' I MENS I ON ALPHA (25) » STALL ( 25 ) » JSTL (25 ) » TALK<20) 

DIMENSION CTDpdp) » CNDP ( 5 ) * CMOp(5 ) 9 e< 6) > ES(6), F(6)» FS (6 ) 

c 

common m22z*cd» ru* ri 9 r29 r3 9 pi* b, bs 9 sa 9 sas» p 
common/ neari/ nrOn»n°l ^ * ■ nz» mz* npres » ir» nt ime» nErr ♦ n d ag» 

i N P H I s NPRImT 

common/ near?./ c 9 a 9 as » d 9 dss sc» gs 9 h 

COMMON/ NEAR 3/ RRP 9 XP » Z> BLD > RB? BR s BTA , TC3LD » TC * RCBRP » APA » 

1 AlF vXprEr >rA 9 XCB 9 Xr 9 LlCBC 9 pmAX 9 PHI » CorJsCgrCS 
common/ neara/ ug 9 uqd, ugd, uv s ucb, vcb* ugp ? uga,gv9 gam s 

1 .. GRV 9 CPP, CPM 

Common/ nEAp5/ ARJ 9AR jp»EpS»RAD»CL»ALPHA9 STALL» jSTl»T \lk 
c 

.101 Format <*RUN numbr*9I5»49X 9 9*page*»i3//-) 

701 Format ( *ru n number *- 9 i 5-9 - »*optiqnAl odtput-#9].7X9*page*9I3//) 

702 Format ( / 17X1HT 9F5 02 » ih ) 5 29X1H c 9 f 5 „ 2 » ih ) ) 

703 Format 1 ihi » 5 X» *pUn n u m°£r*» i«>»:ux»*'Duct surface pressure distribu 

■1TI0N*»10X» *PAGE*9l3//) 

102 FORMAT ( 5X »20A4/'/ ) 

113 FORMAT (9X 9 4HH ( N) ) 

114 format ( 9X5 hsc { n) ) 

115 format ( 9 x 9 4 HC ( n) ) 

116 FORMAT (9X»4HR(N) ) 

117 FORMAT (9X 9 5H6 {|\!)* ) 

118 FORMAT (9X 9 8HSUM A £ N > ) 

119 FORMAT ( 9X 9 9HSUM A ( N ) * ) 

i?u format c 9X9 4HD ( n) ) 

121 FORMAT (9X5HD(w)*) 

122 FORMAT £ 9X6HA ( M »N ) ) 

123 format ( 9X.7HA £ M'»n ) *) 

124 FORMAT ( 14X 9 *EFFECTI'VE CAMBER *94Fl0o&//) 

143 FORMAT £/iQX»*FOURlER Co^INF SERIES CqEF F IC I EnTS* 9 / 15X 5HN»1 » 6 5X4HM« 

11 9 12 // ) 

148 FORMAT ( 7X , *DuCT » •, „ C/D* , 6X4HXP/ C6X3HT /C6X6HRTE/RP4X6HRCB /RP5X4HA 
IP / A/ lCXoF 1 0 0 6/ / ) 

149 FORMAT ( 15X5HALPHA95X1HJS9X2 HJ* >6X8Hj CqS(A) 9 2X8HJ 1 CoS ( A) / 
110X,F10„394F10«5) 

150 FORMAT ( 5X1HNAX4HR/RP4X5HUqD/V7X5HUGD/V7X4HUG/V8X5HUCE/V6X8HGAmmA/r 
IV) 

151 format ( /// 15X 1 5 s2X9«-iterat ions epsilon =*9F9o6) 

152 Format (/ 17X6HINFL0W 9 19X5HbLrtDE/5XlHN4X4HR/RP5X3HU/V8X5HGAM/V7X5HAL 
1PHA6X9HDELTA p/Q) 

153 Format ( /qx^hctp ( d ) 5X6 hctd ( p > 4x.7HC.TD ( p ) * 5X4HCTDP7X5HCTDP1 6XahcnDp 97 
1X4HCMDP ) 

154 FORMAT (2X4H (A) »7llPElIo4)) 

155 format ( 2X4H ( b ) ?7(ipeiio4)) 
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156 


157 

244 

245 

256 

2 c i 

252 
2 5 5 

254 

255 
2 36 

257 
260 
261 
C 


931 

930 


C 


900 


61 

901 

902 


based c 

1_SPEED*) R " 1 A " (b) Cotf-FICIENTS BASED ON 

1 * /20X » * I N * A N N U L . I L n 6 ^ ^ 2 5 1 3 j* LlFT C ° EFFiCIE N T HAS * 

FOR.MAT( 5X»6(ipE13.6) ) 
format ( F 9«4»5 t 1 x 2 (ipeio °3 ) > 1 

rCR:MAT(// 2 X 7 HAZlMUTH/3X8HPHI - /*c(- H P * 

i-ORMAT ( /5X-6.HX/C .ifoVfKi-nn ~ v/ ' 5 3h ^ — ^X * F 7 * 2 s 4 X 4H' 

format (/axahx/c l N) •«THCp t0 un ,ix> ) 

format 1/ 5 xahx/c' ; ? Sxhcp " ’77)? ! 0U I 1 ’ 1 x • 1 

format 1 !; ? 5 x x :h 7 vc . 1 !a p ( ‘{;! ’ 7 ™?!°?”^’ 

Format (/ / /qx*opcgsje>f cqef c ic tents ^ ^ix) ) 

FORMAT! //9 X^PdEc:.oUoe“cocFFHTFmtI f CN PROpEllEr 

FORHAT(I6,F9 0 5,5fipEl2?5n 0N FrEE STpEAf " 

format ( 2X ih* 1 3 s F9 .. 5 s. 5 ! ipe 12 v 5 ) ) 


IF(NTIME «50) 930 , 93 G. 931 

NERR= 1 ’ 

NPRINT= NPRINT+1 
SNALF- S I N ( a L F / PAD ) 

CSALF= C0S( A[_F / RAD) 

ARJV= ARJ/ CSAlF 
ARJVP= AR-JP/ rSiLF 
NCYL= HZ 


IFtNPRlNTo Eq a 
GO TO 90i 


1 o OR * NPR I MT 0 GE 0 1 1 ) GO To 9„o 


NPAG= NPAG+ 1 

Wp I TE (- 6 , 7 01) NR UN 5 

write (69 102) talk 

WRITE( 6 * 148) CD.Xp, 
W R I T E ( 6 » 124) R 0 , R 1 
WRITE! 6? 149 ) ALFs. 

write (6? 702 ) corj? 

WRITE( 6* 150) • 

DO 61 J= l» NZ 
X= UG(J)# GAM 


NPAg 

TC? RRP » pCBRp s AP A 
9 R2 , R3 

ARjV, ARJVP» ARJ 9 AR.JP 

corcb 


Y = UGD(J)* GA.v, 


R = R B ( J ) * p p p 

x 0 = uqd(j)* cor j 

XY- UCB ( J ) * COR.CB 

write ( 6 * 26O) j 9 r,xq*y 9 x»xy,grV(j) 

W R I T E ( 6 9 151) NT I ME , EpS 

! F ( NPRI NToGEo 10) GO TO 902 

GO TO 903 

'.MPAG* NPAG+ 1 

W R I T E ( 6 » 701) NR UN , MPaG 

Write (6, 1 02 > talk 

W.r I T b ( 6 » 148 ) CD > Xp 9 TC 9 RRP 9 pCBRP* APA 
W;RITE (6,149) ALFs ARJV ,Ar j’Vp , AP J > AR JP 
write (6, i43)ncyl 

|k;I T E ( 6 9 115) 

Wfe I T E ( 6 * 244) ( C ( J ) > J= 1 , 6) 


J FREE STREAM 
PROPELLER Tip 

■CCEECc.0 ClmAX 

‘ X ) J 

Tip SpEED*) 

•1 VELOCITY#) 



n n 


90 5 


925 


C 


90- 


924 


62 
6 9 


write (6, 244) ( sc(j) 5 t= 1t A) 

CONTINUE J 15 6) 

R I T E ( 6 s 116) 

k - u 61 

wrue! 6 6 ; 1 541 K >- *- i. *■ 

" 9 I T E ( 6 s 244 ) ( 5 A ( J ) , j= is 6 ) 

WRITE! 6 9 119, J 1 ° 

WRITE! 6 » 244) ( 6^S(Jl i- i < , 

WRITE! 6 s 120) J “ 1 ’ 6) 

write! 6, 244, (d!n) » 1 , 6 . 

WRITEI6 , 121 ) ' 

Write ( 6 9244 ) (dsuj, w= 1 , 6 , 

WRITE {69 113 ) ' ' 

WR I TE ( 65 24-4- ) ( h ( ) - m- 2 . £ < 

^ R I T E { 6 9 122) ’ ' V u) 

DO 924 Has c Y [ 

write (6 -> 244) ’(a( ; v»n,» n= l9 ,, 

WRITE (69 123) 6 

DO 925 Ms l» ncyl 
W R I T E ( 6 9 244) ! AS { M »N ) ,N» l, 6 ) 

NPAG= NPAG+ 1 

WRITE! 6 9 101) NR UN 9 NPaG 

w r 1 t e ( 6 9 102 ) talk 

1485 CDs x Ps TC9 ppp» pCBRP 9 p pA 
ALF ’ ARJV ’ ARJVP - ARJP 

x« 0,0 

Y = 0 o 0 

DO 63 J= 1, NZ 

X= GRV(J)/pt/arjp*3lD/Z 

Y= Y+ X 
9 = RB(J)* ppp 

DELP= GFV<J}* BLD/PI/ ARJP 
IF( . t STALL(J)) 62, 64, 62 

vvpiTl( 6, 260) J , R ,UV ( j ) , G V ( j , 9 A L pHA(j) 9 uE L d 
o U TO 63 ' 

Continue 261) j * R,UVCj)j GV(J), ALPHA (J, 9 OELP 
X «= c o 0 


Compute ducted propeller thrust coefficients 

CTCF= AR J#AR J*PI /APA/S o 0/CSALF/CSAlF 
CTDp ( 1 ) = Y*APA*CSALF*CsAlF 
DO 20 J= l, 6 
24 E(J)= SA { j ) +D ( j ) #CORCB+ b ( j ) * qAm 

CQMa »*CD*p I*GAM#CSALF-* CSAlF 

rTA P l 2) c l (i) " (4 ^°^ E(1) + 2 ^ E( 2) ^2 o*C( 2 )*E(1) 

DO 21 = N- C l 1) * {4o * H(1 ) + 2 o'O*' Sc ( 2 ) * Hti) 
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21 


C 

C 

C 

30 

31 

r 

c 

c 


32 
3 3 


34 


35 


C 

r 

c 


^i LF ;, LTALF+ sctN+1) * h{ N j- sc in)* h ( m+ i) 

CTDp (2) = CTDP(2)+ C { N+l ) *E ( N ) 

CTOp { 2 ) = CTO.p ( 2 ) *C-)N+P I*CO* ( 2 
CTA(_F= -pI*CD/2,0*SNAlF*SmA-lF 
CTDp ( 2 ) = CTDP ( 2 ) + CTALF 
X= RRP* RRP 


-C ( N ) * t ( N + 1 ) 

o*Sc.(i)+ SC l 2 ) ) *SnA|_F*SmA|_F 
CTAlF 


x* ioO-a 0 o/x) 

x= X* DElP* C..SALF* CSAlF 
CTDp ( 3) = CTDP { 2 ) + X 
CTDp < 4 ) = CTDP ( 1 ) + CTDP (2) 
CTDp { 5 ) = C T D P ( 1 ) + CTDP (3) 
DO 22 J= l, 5 


CTDp ( j+ 5)« CTDP ( J ) * CTCF 

COMPUTE DUCT NORMAL FqrCE COEFFICIENTS 


DO 3 0 J = 1 ,5 
CMDP t J ) = ^ 0 0 
CNDP U)= na 0 
SACA=SNALF*CSalF 
I F (SACA ) 31 ?23 > 31 
MZZZ=-1 


compute f(.n;) and f-ptor(n) Fourier coefficients 
N F = 5 0 

CALL VTXrmg(CD*Xp».nF»RB*C*U6D*XppES»dJ 
MZZZ=1 
DO 32 J« 1 , 6 
FS(J)=UGD(J) 

Ft J)»UGD( J+6) 

DO 33 J = 1 9 6 

ES(j)=SAS(j) +DS(j)*C0RC3 +(6S< j)+FS{j) )*GA m 
CND p (2 ) =SC ( 1 ( 4o*ES ( 1 ) + 2 v.*ES( ? ) ) +2«*Sc ( 2 )*ES'( 1 ) 

DUM«0.0 
DO 34 J = 2 s 5 

DJM-DUm+SC { J+ l ) *ES ( j )-SC ( j )*ES( J„l) 

CNDP ( 2 ) = { CNDp ( 2 ) + DUm ) *p I *CD*S aCA/ 2. 

CNDp ( 3 ) aC ( 1 ) * ( 4u*GS't 1 ) + 2 o*GS ( 2 ) ) +2 o*C ( 2 ) *GS ( 1 ) 

DUM«0 o 0 
DO 35 J=2,5 

OUm = DUm+C ( J + l ) *GS (.J )~C ( J ) *GS ( j+l ) 

CnDp (3 ) - ( C m D p ( 3 } +DUm ) * p I *CD#S ACAX 2 a *GAm 
CNDp (4) =p I*CD*SACA* ( 2o*Sc ( 1 ) „SC (2) ) 
CNDp(l)=CNDpt2V + CNDp{3)+CNDp(4) 

COMPUTE DUCT MOMENT COEFFICIENTS 

CMDp (2 ) =C ( 1 )* ( 4o*GS ( 1 ) +4o*GS ( 2 ■> +2o*GS< 3 ) )+C ( 2 )*' ( GS ( 2 ) -GS ( 4 ) ) 

1 +C ( 2 ■)■■*( 2 o*GS ( i )-QC ( 5 ) } ( 4 1 ^ t GS t 2 )“GS(6 ) >+C(5) *GS(3 ) 

2 +C ( 6 ) *GS { 4 ) 

CMDp (2) =CMDP t 2 ) *P I V4 »*GAm*SACa*CD*CP 
CMDp (3 ) =P I/2e*CD*CD*SACA*(2o*SC ( 1 )+SC (3 ) ) 

Cm Dp ( 4 ) = 3C ( I ) * ( 4 o^-ES ( 1 )+ 4 a *E^( 2 ^ / , 2 j *ES( 3 ) ) +SC 1 2 ) * ( to ( 2 )*»ES ( 4 ) ) 
1 +SC ( 3 ) *( 2 o*ES 1 1 5-tS(5 ) ) +oC t 4.) *< ESI 2 )«ES ( 6 5 ) 
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36 


37 


38 

23 


70 


72 


71 

73 

74 

75 


C 

C 


2 +SC(5)*ES(3)+SCC6)*ES(4) 

CM Up (4)=U4DP(4)*PI/4c*SACA*CD*CD 


DO 36 J=l»6 
E ( J ) =E ( J ) A F ( J ) #GAM 

CM^P (-5 ) =SC ( 1 ) *'< 4»*E ( 1 }+2 **£{ 2 ) ) +'2'o*6>C ( ? ) *E ( 1 ) 

DUM-OoO 


DO 3 7 J = 2 , 5 

DUM=DUM+SC ( J+l ) *E E J ) -SC E J ) #E ( J+l i 
CMDp E 5 ) =>**( CMOp ( 5 ) +DUM ) *P 1/2 «, #SACA*CD 
DO 38 J = 2 5 5 

CMDp ( 1 ! = C.MDp ( l ) +CmDd ( J ) 

CONTINUE 
PRINT 153 

PRINT 1 54 > E CTDp ( n ) >N-l 55 ) ’CMDP ( 1 ) sCMDp ( i ) 
CNDpU)=CNDPEl)*CTCF 
CmL>P( 1)=CMDP(1 )#CtCF*RRP/2oO 
PRINT 155? ( CTDp { N! 5 ;-N=6*l<J ) »CmDP ( I ) *CmDP ( 1 ) 
PRINT 156 
N=0 

NTOT =0 

DO 71 J«1»NZ 
I F ( STALL ( J ) )7l,71<,72 
N«N+1 
J,STL(N)=J 
NTOT =N 
CONTINUE 

IF (NTOT) 75 5 75 ,74 

PRINT 157» (jSTlIN) *N« l*NTot) 

CONTINUE 

I F ( N P H I ) 9 9*99,80 


COMPUTD DUCT SURFACE prEPSUdE COEFFICIENTS 


30 CALL PRESS 

DO 81 I =1 9 NPHT 
DO 81 J = 1 9 IP 
IF (XPRES ( J )~XP)81 ,82 982 
82 CPP( I 5j)»Cpp( I »J)+.DELP*CSALF*CSALF 

31 CONTINUE 
IFlNPRESoLEol) GO To 85 
DO 3 3 I = 1 9 N P H I 

DO 83 J=1,IR 

CPP ( I , J ) =CPP ( T 9J)*ARJV*ARJV/2 o 0 
8 3* CRM ( I > J ) =CPM ( I , J ) *AR JV*Ar JV/2 oO 

8 u NPAG a NPAG+l 

PRINT 703 9 NRUN , NPAG 
PRINT 102.TALK 

PRINT 148 9 CD » XP a TC , pRP 9 RCSrP , APA 
PRINT 1 49 , ALF , ARJV , AR J VP , AR J., AF JP 
PRINT 250 9 ( PHI ( J ) 9 J=1 9 NPH I ) 

GO TO (1,2,39495] 9 NPH I 

1 PRINT251 
GO TO 6 

2 PRINT 252 
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GO . TO 6 

3 PRINT 253 
GO TO 6 

4 PRINT 254 

GO TO 6 

5 PRINT 255 

6 CONTINUE 

DO 86 J=l»ip 

85 PRINT 245»XpreS( j|»( CpP { I , j ) , CPM II » J ) , 
IF(NPRES 0 LEol) GO TO 6 ' " 

PRINT 256 
GO TO 99 

8 PRINT 257 

99 RETURN 

END 



I = 1 » NPH I ) 


\ 



$ I B JOB MAP 

$ I BFTC MAIN 

DIMENSIONS ALPHAr (48 ) sHqLVEl { 16 3 

DIMENSION XV(4 8 .) ,0ELTftX(48) »Xp(64) »YV(48) *S<48> »Yp(64) 

DIMENSION Nl( 10) »L1 (10) 

DIMENSION DWASH(48') »C(6) 

COMMON/C.OM1/HolVEL»NNN 

common /com 2 / n»nn»m»xv»deltax>x'p»yv,s,yp 

C 0 M MON / COM 3 / C B »NF *N1 »L1 
• D I = ATAN ( 1« 0 ) *4 * 0 

C , OF VqrTICESs M = N0° OF SpAN STATIpnSs NN= NQ° OF CONTROL POINTS 

C ir N = N N HOLE IN WInG ABSENT 9 WInG FORCES CALCULATED 

R t AD 501»N»NN,M 
5 a 1 FORMAT {3i 3) 

PRINT 502 »N*M*NN 

502 FORMAT ( IX »*N = * » 15 »*M ®* » I 5 » *NN=* , I 5 ) 

READ 1 0 0 9 ( XV ( I) ,DE'lTAX( I ) ,Xp( I) »YV{ I ) » S ( I ) *Yp(I ) »I*i»'|) 

100 Format ( 6Fi0c6) 

PR I f\|T 3^0 s ( I sX v’( I ) J DElTAX ( I ) sXP ( I ) 9 YV { I ) 90 ( I ) 9 Yp ( I ) 9 1 =1 9 N) 

C CSAR = ROOT CHORD 

CBAR«40 o 0 
CB“CBAR /5 o 0 
DO 203 1*1, N 
XV ( I )=XV( I )/cB 
DELTAX( I )=DELTAX( I )/CB 
XP ( I ) =XP ( I )/Cb 
YV ( I ) = YV ( I ) /CB 
S ( I ) =S ( I ) /CB*0 o 5 

203 YP( I)=YP( I )/cB 
PRINT 301 

PRINT 3'JO j ( I »XVt I ) ,DElTAX( I ) »XP ( I ) 9 Y V ( I ) »S ( I ) ?Yp{ I ) 9 I =1 9 N ) 

300 F0RMAT(lX»I5»6F10o5) 

30 1 format 1 1X9 1*96*9 #x v ( 1 j * 94X 9 -*delTax* 9 4X 9*xp # » gx 9-k-yv# 9 sx » *2^* »9X » 
1#YP*//) 

I F ( NN— N ) 20 9 2 1 9 20 

20 NF=N+1 

READl01»(Xp( I) »YP(I )*I«NF*NM ) 

1 a 1 FORMAT ( 2F10, 5) 

PRINT 313 

PplNT312 9 ( I 9 Xp ( I ) 9 YP ( I } 9 I=NF?NN) 

DO 204 I = N F 9 N N 
Xp( I )=XPU J/CB 
YP( I ) “YP ( I ) /CB 

204 CONTINUE 

PRINT 312 » (I »XP ( I ) 9 Yp ( I ) 9 I =nF »NN ) 

312 Format ( 1X» I3 9i0X 9 El5o8 9 10X,El5.,8> 

313 FORMAT! lX9*I**9i2X9*XPLl:)* 3.0X »*Yp ( I ) # 9 *VAlUES FqR HOLE*) 

21 CONTINUE 

c ifnnn is +vf only vElo ci ty in the hole calculated 

C ' MNMIS 0 ONLY FORCES ARE CALCULATED 
C MAIN IS M/E BOTH ARE CALCULATED 

C Nl=NO « OF VORTICES AT A GIVtN SpAm STATION PLUS ONE 
C LI® NUMBER oF THE FIRST VqrTEX AT THE SpAN STATION 

C K2- N3“ OF THE lAST VorTEX AT THE GIVEN S-pAp) STATION 

READ 503 9 (Ml ( I 5 9l=l rM) 



U U i\J 


503 FORMAT ( 1OI5) 

P R I M T 5 0 7 

PRINT 5^6 * (I »N1 ( I ) »L1 ( r ) * I = X 9 M ) 

5 ° 6 F OR MAT { 3 ( 4 .X 9 I 5 ) ) 

^7 format ( iX9*statiqn*»2X 9 *ni*>7Xs*lt-*i 

A 1 * 1 o u 

Do 11 1=1948 

11 _ ALPHAR ( I )=Ai*PI/i8Oo0 

L ,_ Ir 15 zero or - v e then subroutine wing 

c i,u ;;cXl 15 +VE THEN subroutine effect IS CA| I eu 
,vEXT = io 0 p 
I F ( N EXT ) 1 9 1 9 2 

DATA 

ALF=88 j 0*p I / isO „ C 
GAM=61o53817*CO s <ALF) 

READ 1^23 t C C I ) »I=i » 6 ) 

i °2 format ( 6 F 1 O 05 ) 
c. 

CALL EFFECT ( AlF 9 GAm 5 C 9 D^ASH ) 

DO 12 1=1 9M 

1 2 ALPHAS ( I ) — PW/\S,H ( I ) +ALdHAr ( I ) 

PRINT 5O9 

PRINT ?u 8 , (ALPHAP ( I ) 9 1 — 1 9 N } 

5 U 9 FORMAT ( 1 X 9 * THE Down aw MATRIX*//) 

5^8 format ( xo c ?x , Fio 06 ) ) 

NNN=— IP 

1 CONTINUE 

call WING (ALPHAr ) 

STOP 

end 


S CALL f 1 
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SIBFTC WING 

SUBROUTINE WImG(ALPHA) 

ulHtNblON x y(4 S) » d ^LTaX(48) ^Xp( 6 4) 9 YV( 48 ) ,YP(6A) 9 S(4S) »F(43 , 43 ) 

1 5 ALl (4b ) 9 AMI ( 48 ) ,B( 48 >1 ) 

DIMENSION Ai_PHA(48) ’HOL v El( 16) 9 HF( l6 »4S) 5 GAmWG(43) 

DIMENSION CLL(11)»CML(H) ,Nl(lO),n(10) 

DIMENS I ONCLX ( 12 ) 

COMMON/COM 1 /HqlVEl » nnn 

C 0 MON /Comp/ N » v ,M > m » X V > D E L T A X » X p » Y V , S , Y p 
COMMON/ COM3 /CB »MF ,N1 ,L1 
:' = ATAN(l n 0)*4 o 0 

■ 200 1=1 


2 00 J=1,N 
(I)-XV(JJ 
• ( I )“YV( J) 

YOAR =YP ( I )+YV ( J ) 

AS 1 = Y + S ( J ) 

AY2=Y— S ( J ) 

AB3=SQR T (X**2+ ( Y+S ( J ) ) **2 ) 
AB4=SqrT ( X** 2+ ( Y-S ( j ) ) ) 
A35=YBAR+S ( J) 

AB6=YBARS»S ( J) 

AB7 = Sq.R I (X#*2+A85*#2 ) 
AB8«SQR T ( X**2+AB6**2 ) 
IF(XsEQoOoO)Fll a OoO 
IF(XoEQo A o0)F13=0o 0 
IF (XoEQoOoO) GO TO 215 
FllaABl/ (X*AB35 
F13-AB2/ (X*a'B4) 

215 F 1 2=1 o / AB 1* ( lo+X/AB3) 
F1A=1 0 /AB2*( 1 oa X/ABA) 

F 1 = F 1 1 + F 1 2«-*F 1 3 a F 1 4 

IF {XoEQoOcO}F21=OoO 
IF (X<.EQo0e0)F23=0 o 0 
IF(XcEQo0 o 0)G0 To 216 
F21=AB5/ (X*AB7) 

F23=AB6/ { X*AB8 ) 

216 F22ilo/AB5*<lo«X/AB7* 

F24«1 0 /AB6*< 1 0 + X/AB3 ) 
F2-F21+F22— F23-F24 
I F < I „GToN.) GO To 202 
FU,J) =F1+F2 


GO TO 200 
202 11= I— N 

HF ( 1 1 9 J ) = F1+F2 
200 CONTINUE 

PRINT 5l rt , (F< 1*1) *I=1»N) 
510 FORMAT (lX*13Fl* 0 5} 

D02q 8 1 = 1, N 
2*8 B(I»1)=0 o , 

CALL M A T I n V ( F »48»NjB»1»1 ) 
K-NN-N 

DO AO 1 I = 1 , M 
A L 1 ( I ) =0 0 0 
GAMWG ( I ) =0 o 0 
I F { IoGT c k) GO TO AOl 




HOLVEL ( I )=0«0 
4 A l CONTINUE 

DO 402 1*1, N 
DO 400 J = 1,N 
F ( I » J) =8 o 0*PI*F { I ,J] 

IF ( NNN ) 1 s 2 » 1 

1 GAMWG ( I ) =GAMWG ( I ) +F ( I , J ) *A LPHA ( J ) 

5 IF (NNN) 4,4 , 3 

3 GO TO 400 

4 CONTINUE 

2 F( I » J ) =F { I »J)/DELTAX( I) 

AL 1 ( I ) = A L 1 ( I) +F ( I »J)*ALPHA'( J) 

4 n 0 CONTINUE 

IF(NNNoGToO) GO TO 402 
AX.LE=YV( i )/6„ n 
4 a 2 CONTINUE 

AMI ( I ) = ( XV ( I >-AXLE)*ALl ( I ) 

IF (NNNcEOoO) GO Tq 14 
DO 403 1=1, K 
DO 403 J=l ,N 

HOLVEL ( I ) =H0lVEl ( I )+1oU/ 8u u/dI*HF( I *j)*GAmW6{ J) 

403 CO/TINUE 
PRINT 509 

P R I N T 5 0 8 9 (HOLVEL (I ) »Iel,K> 

508 F0RMAT(1X»8(1X»E15.8)) 

5 a 9 Format ( ix,*velocities i,m the hoi e due to the wing*) 

I F ( NNNoGT o 0 ) GO Tq 999 
14 CONTINUE 

PRINT 307 

PRINT 306 » ( AlI ( I ) » 1*1 *N> 

307 FORMAT (IX** LMATRIX)*) 

306 FORMAT (8( 1X9F15o8) ) 

PRINT 308 

PRINT 306* (AMI ( I) »I=1»N) 

308 FORMAT (IX,* M MATRIX*) 

C To GET Cll (LOCAL LIFT COEFF) 

DO 888 I J = 1 , M 

CALL VARS I M ( N 1 ( I J J 9 XV » CLL ( T J ) » ALI *L1 ( I J")' » ^ » XP »N »N » NN ) 
CALL CURFIT(XV.YV,L1<IJ).AL1*AA1> 

13 CLL ( I J)=CLL( IJJ+AA1 

888 CONTINUE 

PRINT 310 

PRINT 3 O 99 (CLLU ) »I = 1*M) 

309 FORMAT ( 1X»8E15o7) 

310 F0RMAT(1X,*L0CAL LIFT COEFF *) 

C To GET CML (LOCAL MOMENT COEFF) 

DO 389 IJ=l,M 

CALL VArS I m (N l (I J) »XV9CML ( I J) >Am1»l1IIJ) 90 ,Xp »N »NN ) 
IL = Ll(IU) 

AXLE=YY(lL)/6.0 

5*(XV( III) -AXLE) *AM 1 ( EL) 

CML ( I J ) =CML ( I J ) +AA 1 



889 CONTINUE 
PRINT 311 

311 format ( ix ♦* local moment coeff*) 

PRINT 309».(CML<'n»r = l»M) 

C TO GET CL AND CM 
N A=M+1 

DO 220 l=l, N a 
IF(IoEQoNA) GO TO 221 
;<3=L1 ( I ) 

CLX( I ) =YV ( k3 } 

GO TO 220 

221 CLX(I)=60„0/(CB*5o0) 

CLL ( N.A ) = A e 0 

C-ML { NA ) =0 o 0 
220 CONTINUE 

Q = ( C L L ( 2 ) “C L L ( 1 ) ) / ( CLX ( 2 ) -CLX < 1 J ) 

Ol«tCML(2 )-Cml(1) )/(ClX(2)-C L X(1) ) 

CALL VArSIm<NA»ClX»Cl»ClL»1>1»Xp>12»11»NN) 

CALL VARSIM<^A»CLX*CM»CML* 1*1 »XP »12»11»NN) 
CL=CL+ U ° 5* (0*ClX ( 1 ) +Cll { 1 ) )*ClX ( 1 ) 
CM=CM+0o5*(Ql*CLXll )+CML( 1 ) )*CLX( 1) 

PRINT 222 »CL »CM 

222 format ( ix > * cl=*»ei5.>7»* cM**»ei5„7) 

999 RETURN 

END 



TIBFTC INVERS 

SUBROUT I NE MAT I NV ! A » NN »N * 8 ,MM ,M ) 

£> $KSN S'OM A(nN, HN'j , J Pi V07 6*8) IfC /DT'-4? )Nt!£y{4jE£2^gS£f5p 

C INITIALIZATION " • 

15 DO 20 J = 1,N 
20 I P I VOT ( J ) =0 
30 DO 550 1 = 1 , N 

c search for pivot element 

40 AMAX =0 o 0 
4 5 DO 105 J=1»N 
50 I F ( IPIVOT! J)-l)60*105,60 
SO DO 100 K= 1 s N 

70 IF ( IPIVOT {<)-l ) 80, 100, 74„ 

8r> IF ( ABS ( AMAX )»ABS ! A ( J»K ) ) >85 9100,100 

65 IROW=J 
TO I COL UM = K 
9 5 AMAX = A ( J ,K ) 
iCO CONTINUE 
105 CONTINUE 

110 IPIVOT ( ICOLUM) =IPIVoT( ICqlUm)+1 
c interchange rows to pot pivot element on d+agona 
130 I F { IR0W*»IC0LUM) 140,260,140 
140 CONTINUE 
150 DO 2 00 L = 1 , N 
160 SWAP = A( I ROW , L) 

170 A ( I ROW , L ) = A! ICOLUM , L ) 

200 A ( ICOLUM 9L ) =SWAP 
205 I F ( M ) 26„, 260, 210 
210 DO 2 50 L= 1 , M 
22 0 SWAp=B ( I ROW , L ) 

2 30 B ( I ROW 9 L ) =B ( I COLUM , L ) 

2 50 B< ICOLUM, U=SWAP 
260 INDEX! 1,1)= IpqW 
270 INDEX! I ,2) = ICOLUM 

310 PIVOT! I)=A( ICOLUM, ICOLUM) 

C DIVIDE PIVqT roW BY PIVOT ElEmEN.T 

330 A! ICOLUM, ICOLUM) =loO 
-40 DO 350 L = 1 , N 

% ku A ( IColUm » L ) =A ! I Col u m , L ) I v oT ( I) 

355 IF ( M ) 380, 380, 360 

3 60 DO 3 70 L = 1 , M 

3 7‘^ ti ! I ColUm » L 5 =6 ! I CoL u m , l ) J ? I VqT { I ) 

C . REDUCE NON^PIVOT ROWS 
380 DO 550 L 1 = 1 ,N 
39U I F ( Ll“I COLUM) 4 OO 9 55^9 400 
400 T=A( Ll , ICOLUM) 

420 A(L1» I COLUM) =0.0 

430 DO 450 L = 1 , N 

■ 450 A ( Ll , L ) =A ( L 1 » L )~A ( ICOLUM 9 L) *T 

455 I F ! M ) 550, 55q, 460 
460 DO 500 L = 1 9 M 

5 00 B (Ll ,l*)=B ( L1»L)-Bt ICOLUM »L)*T 

550 CONTINUE 

C INTERCHANGE ColUMNS - 

600 DO 710 1=1, N 
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610 L-N+1~I _ F4020074 

62^ IF ( InDLX ( l 9 1 ) — I N'OEX ( L?2 ) ) 63 A 9 710? 63'^ F402007f| 

630 JROW«I NDEX t L * 1 ) f~4Q200 

6 40 JCOLUM= I NDEX ( L.*2) !*4Q 20071 

6d0 DO 705 K = 1 9 N F402007S 

6 6 ; J SwAP = A ( K > JROW ) <:r 40 2007 f: ? 

6 / 0 A ( .< ? JROW ) = A ( K $ JCOLU M ) r 4.Q 200 3 0 

7 v ^ A t K ? JC 0 L Um ) = S W A P p 40 2 0 0 3 7 

LOiNjTINUt F 4020012 

0 CONTINUE 

C - -TURN 
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VARS I M ( N » X , SUM , F-U s> L 1 » J1 » Xp » Il»I2»I3) 

C IF J 1 »0 THE N THc L AST VAi_Un Uj THE LAST I^TErVAl =OA,njD X =Xp{'2 + l) 

c n should be odd so that the mo o of intervals are even 
DIMENSION X( II J ,FU( 12) ,Xp( 13) 

SUM = 0o 0 
N 1 = ( N - 1 ) / 2 
H=X(LI+1)-X(L1) 

H=ABS (H) 

DO 100 11=1, NX 
L = 2*h- 1 
L2=L l+L-1 

I F ( J 1 a GT o 0 ) GO TO 300 
IF ( I I oEO.Nl )AK1*XP(L2+1) 

IF (I I oEQoNl ) A k2 = OcO 
I F ( I I otQoNl ) GO TO 301 

300 CONTINUE 

AK 1*X ( l2+2 ) 

A K 2 = F U ( L 2 + 2 ) 

301 A LP 1= ( X ( L2+1 )— X ( L2 ) )/H 

A L P 2 = ( A K 1 W X ( L 2 + 1 ) J /H 
ALP 1 =ABS { ALP 1 ) 

ALP2=A.BS(ALP2) 

Sl = FU (L2 }*ALP2* ( ALP l**2 + 0 « 5* ALP l*ALP2-0 » 5* ALP 2**2 ) 

S2=0 o 5*FU(l2+i ) * ( ALP1+ALP2 ) **3 

S3=AK2* ( ALP1**3-*ALP1**2*ALP2, n aLP2**2*ALP1 + 1o5*ALP1*ALP2-1o 5*ALP1 
1**2) 

S4=H/ (3oO*ALP1*ALP2 )*(S1+S2+S3 ) 

SUM=SUM+SA 

100 continue 
return - 

END 
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'ETRfTC cur 

SUBROUTINE CurFIT (X »YV»L1 »AL »ArEA) * 
c USING CURVE A_*X**<~0<,5)+B*X** ao 5+C*X**1. 'WHERE X=0 IS THE i.E 

DIMENSION ALF ( 3s>3 ) »X (48> »AE6 (r» 1> »Al t AP ) 

XLt = YV{ LI ) /<2 0 *3, ) 

DO 100 1*1,3 
K=L 1+ I —1 
DO 101 J*1 , 3 

a & f c i » j ) = (x ( k)-xle ) ** (float { j)-io5 ) 
ioi continue 

AEB ( I s 1 ) = A L { K ) 

100 CONTINUE 

CALL MAT I MV (AEF ,3 » ? » AE3 » 1 » 1 } 

XLE=X(Ll) n XLE 

AREA=AE£( 1,1) * X L £ 0 - 5/0 *5+AEB(2»l ) *XlE**1<> 5/1 * 5+AEB ( 3 , i j *v t c 
1**2 o 5/2 „ 5 L 

RETURN 

end 



« I 6 F T C EFFECT 

SUBROUTINE EFFECT (AlpHAsGAMjC ,DWASH ) 

blM^NSION XXp ( 2 1 ) » Y ( 2 1 ) * pH I p ( ? ] ) 9 V 2 ( 2 1 ) sF(2l921)9C(6) 

DIMENSION DWASH ( 48 ) 

DIMENSION XV ( 48 ) s ^ElTax (4q ) , xp ( 64 ^ 9 YV(48 ) 9 S ( 48 } s YP ( 6^ ) 
Common /com 2 / n>nn»M 9 Xv»dei_tax»xp,yv,s,yp 
common mzzz 

DUm=OoO 

MZZZ =0 

CALL EllIpS(dUm 9 DU M 9 DUM) 
p I = 4 oO*ATAN d, 0 ) 
c 

C DATA 

CD = 0 o 1 999 
XF c 20 o 0 / 40^0 
YF= 26 o 5 / 40 o 0 
CROOT = 40 0 0 

R TEb 20 o 125 / 2^0 

C 

c 

PRINT 1 O 5 » ( C( I ) , 1 = 1 
i °5 format ( ix »*c ( 1 > * * 96 Ei 5 o 7 ) 

DO 999 JK = 1 ,N 

RR=SQRT ( (Xp ( JK )-XF)** 2 +<Yp< JK)-YF)** 2 ) 

RR=RR*CROOT/RTE 

PHIasABS ( ( X p ( J x ) — X F ) / ( Yp ( JK ) — Y F ) ) 

PHI =AT^n ( pHl ) 

I F (Xp ( JKJ-XF ) 31 94 O 933 
33 I F ( Yp ( JKJ-.YF) 37 94O939 
37 PHI = lo5*PI —pH I 

GO To 40 

39 PH I =P I *0 0 5 +pH I 

GO TO 4 O 

?1 IF (Yp( JKJ-YF) 34940,36 

M PHI = PI* 0 o 5 _pH 1 

Go To 40 

4 pH I = 1 o 5 #P I+PH I 

GO TO 40 

continue 

X = 0 o 0 
NK= 21 . 

C 1 = CD *2 ® 0 

c PRINT 30 l 9 JK 9 RR 9 PHI 

C 3 O 1 F 0 RMAT(iX,* JK=*»I 59 * RR = * 9 El 5 o 7 9 -sf PHI=* 9 Ei 5 o 7 ) 

STEp=Cl/ FLOAT (NK- 1 ) 

XXp ( 1 ) =Cl /2 oO 
STEp 2 = 2 ° 0 *pl /F lOAT ( NK- 1 > 
pH Ip ( 1 ) =0 * 0 
Do 1 °° 1 = 1 9 NK 

DO lOO J=l »NK 
I F ( J o EG » l ) GO To 200 
pH Ip C J)=pHlp( j-1 ) +STEP 2 
200 continue 

IF(joGTol) go To 20 l 
I A* I — 1 

IFlIoEQoD GO TO 101 
XXP (I )=XXP ( IA)— STEP 

101 continue 

Tl = 2 ® 0*XXp ( I ) /Cl 
THETA=A R CoS(TD 
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101 


201 


100 

C 

C50 9 
C 

C505 


C 

C500 


203 

202 


3^2 

999 


I 1 


XXR( I )=XXP( I A ) —STEP 
CONTINUE 


T 1 = 2 o 0*XXp ( I ) /Cl 
THETA=ARC 0 S( T 1 ) 


6A„D=C(l)»TAN(PI/ 2; 0-THt T A/ 2 „0 , + c (2)*SI N( THETAi + C( 3 ),SI N(2 .0*THETA 

^F(aLTti) n ooo? )+C< 5) ^ InUo0 ^ THE ' TA)+C<6, ^ SI ^ (5o ^ THETA > 

IF l ABS (■ T 1 ) 0 GT o 0 o 999 ) GAMD=0 o n 

51=GAM0/ (pi*2oa*SQRT { (X~XXP( I ) ) **2 + ( PR+ 1 » A ) **2 > ) 

1 ( 4 ° U*do/ { { X — XXp ( l ) ) #■&? + ( pp+l o 0 ) **2 ) ) 

CALL ELLIPS(A<P,TK,T'E) 

C0N!lNUE IT< ‘ <1 " i '' +2 ”°* <BR ' 1,O,/ ‘ lX “ XXP ' Ill ** 2+ ‘ l,R - 1 -" 1 »2II*TE I 


tl-C0S(PHIP( J) ) * ( P R*COS ( PH I -p h I P { J ) )_1 0) 

E?i!*7=GAMjlEl/E 2 :Cnr 1 '°* 2 '" <fiRXCOSIPHI - PHIP(J,) 


CONTINUE 

PRINT 5 0 9 s ( V { I ) 9 I = i s n K ) 

format { 1 x 9 -*?• y ( i ) =># $ s e i 5 7 > 

PRINT 505 , (XXP< n 9 1 = 1 , NX) ' 

FORMAT(lX.*XXp APE — * , l if 1 A 0 5 ) 

CALL VARS I m INK »,XXP»UGM0»Y » 1 » 1 ,Xp »NK jNM ) 

PRINT 50 A »UGMD 9 GAM 9 GAMP 
FORMAT ( xx 5 * UGmO=* 9 E 1 5 « 7 9 # GAm=* » E l 5 e 7 9 * 

ugmd=ugmd*gam 

DO 202 I=1,NK 
DO 2 03 J = 1 5 n K 


GAmDs* 9 Ei 3 07) 


Y( J)=E( I 9 J ) 

CONTINUE 

Call VAr s I ,m ( NX 9 PH I P 9 v 1 9 Y 9 1 9 1 9XP 9 nk 9 nk 9 NN ) 

V 2 { I ) = V 1 

continue 

CALL VArS I m t NK » XXp 9 UGALF 9 V 2 9 l 9 1 9 Xp 9 w j< 9 m k 9 1\! m 5 
UGALF=-UGALF/{ 4 oO*pl )* at An (ALPHA) 

PRINT 302 9 UGMo 9 UGALF 

format ( ix»* ugmd**»Ei 5 07 »* ugalf**»e* *><>7) 

DWAS'H ( JK ) =UGMD+UGALF 
DWASH ( JK ) a-DWASH ( JK ) 

CONTINUE 

RETURN 

END 
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JIBFTC 03 

SUBROUTINE ELL I PS ( AkSH!»Tk»TE) 

c 

c : uc. 0 li_lips --table look ud og elliptic integerals 

r 

v. 

D TK EL" SION CKK( 100),CK(100) ,CE'(10 A ) 

r r n«! mzzz 

!F miZ-l ) 73 » 10 5 3 
In CCS.TI NUE 

i 

<-kk= argument of elliptic integrals 

-kk !i )=o 0 oo 

U;-C (:' 1 = 0, „1 
c (3 1 = 0 0 02 
•IKK (4 )=0 ,03 
CKK (5 ]=O c 04 
CKK (6 )=0 ca 5 
CKK (7 )=0„ n 6 
CKK (8 |=0 en 7 
CKK (9 )=0„ n 8 
CKK <1 0) = no 09 
CKK <1 1) = n ,l0 
CKK <1 2) =0,11 
CKK U 3) =0.12 
CKK <1 4) =0, 13 
CKK (1 5) =0 0 14 
CKK <1 6) =0,15 
CKK II 7) =0 o 16 
CKK <1 8) =0 o 1 7 
CKK <1 9) = 00 18 
CKK <2 0) =0, 19 
CKK <2 1} : O0 20 
CKK (2 2) =0 o 2 1 
CKK 12 3) =0 0 22 
CKK <2 4) =Oo 23 
CKK <2 5! =0 o 24 
CKK <2 6) =0,25 
CKK (2 7) =0 o 26 
CK.K <2 8) =0, 27 
O IK <2 9) =G c 28 
CKK <3 0) =0 o 29 
CKK <3 1) =0,30 
CKK <3 2) =0,31 
CKW3 3)= no 32 
CKK 43 4) = 0 , 33 
CKK 43 5) =0,34 
CKK 43.6! = n ,35 
CKK 42 7) =0,36 
CKK 42 -8) =0,37 
CKK43-9) =0,38 
CKKl-40 ) = 0 o 39 
CKK 44 1) =0,40 
CKK 44 2) =0,41 
CKK 44 3) =0,42 



CKK(44)=0 o 43 
CKK ( 45 ) = 0„44 
CKK ( 46 ) =0 o 45 
CKK { 47 ) =0 o 46 
CKK(48 ) =0 o 47 
CKK ( 49 ) =0 o 48 
CKK( 50) = no 49 
CKK(51)= no 50 
CKK { 5 2 ) = 0 o 5 1 
CKK { 53 ) = A „ 52 
CKK ( 54 ) =0 o 53 
CKK(55)=0 o 54 
CKK< 56 )=0 o 55 
CKK ( 57 ) =0 o 56 
CKK { 58 ) =0 0 57 
CKK( 59 ) = 0*53 
CKK { 60 ) = A „ 5 9 
CKK(61)=0 o 60 
CKK ( 62 ) = 0 0 6 1 
CKK(63)= n ,62 
CKK(64)=„,63 
CKK(65)=0o64 
CKK ( 66 ) =0o 65 
CkK( 67)=0 o 66 
CKK ( 68 ) = A o 67 
CKK ( 69 ) =0 o 68 
CKK (70) =0 o 69 
CKK { 71 ) =0o 70 
CKK(72)=0 o 71 
CKK ( 73 ) =0® 72 
CKK ( 74 } = n a 73 
CKK ( 75 ) = 0 o 74 
CKK(76)=0 o 75 
CKK ( 77 ) =0 o 76 
CKK(78)=0 o77 
CKK(79)=0o 78 
CKK( 80) =0o79 
CKK(81)=0„80 
CKK(82)=0*81 
CKK(83)=0o82 
C K K ( 8 4 ) = 0 o 8 3 
CKK ( 85 ) =0 o 84 
C K K ( 8 6 ) = „ <, 8 5 
CKK ( 87 ) = 0 o 86 
CKK (88) =0,87 
CKK( 89 >=0c88 
CKK £ 90 ) =0 o 89 
CKK ( 91 ) =0 o 90 
CKK(92)= ao 91 
CKK(93)=0«92 
CKK ( 94 ) =0 c 93 
CKK( 95 ) =0o94 
CKK ( 96 ) =0 o 95 
CKK ( 97 ) =0 o 96 
CKK ( 98 )=0o97 
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CK(52)=1„862641 
CK 03 )=1 o 87140 a 
CK £ 54) =1 0 880361 ' 

CK ( 55 ) = 1 0 889533 
CK ( 56 ) =1 o 898925 
CK £57)=1, 908547 
C'< £ 5 3 ) = 1 o 9 1 84 1 0 
C. '<!, t 5 9 ) — 1 0 9 2 8 5 2 6 
CK (60) = 1 0 938908 
C !<: ( 6 1 ) = 1 0 9 4 9 5 6 8 
C|< ( 62 ) =1 0 9605 21 
CK ( 63 ) =1 0 971783 
C K. ( 6 4 ) = 1 0 9 8 3 3 7 1 
Cy ( 65 ) “1 0 995303 
C K ( 6 6 ) = 2 0 0 7 5 9 8 
CK ( 67 ) “2 0 020279 
CK £68) =2 O 033369 
CK £ 69 ) =2 o 046894 
CK170)=2 o 060S82 
CK(71)=2 o 075363 
CK (72) =2.090373 
CK (73) =2„ 105948 
CK.(74)=2 0 122132 
CK (75)-2 0 138970 
CK(76)=2 0 1565l6 
CK(77)=2 c 174827 
CK ( 78 ) =2 0 193971 
CK(79)«2 o 214022 ' 

CK<80) =2.235068 
CK<81)=2 o 257205 
CK£82) = <l s 280549 
CK(83)=2 o 305232 
CK(84}«2.331409 
CK(35)=2 0 359264 
CK l86)«2 o 3890i6 
CK ( 87 ) “2„ 420933 
CK(88)»2 0 455338 
CK<89)=2. 492635 
CK(90)*2 o 533335 
CK(91)=2 o 578092 
C!<(92)=2 0 627773 
CK(93)=2 0 683551 
CK ( 94 ) =2 0 747073 
CK<95)-2 o 820752 
CK ('96)=2. 908337 
CK £97)=3 0 ^16H2 
CK £ 98 ) «3 0 155875 
C K £ 9 9 ) = 3 c 3 5 4 1 4 1 
CK ( 100 ) =3 ,>695637 

CE=ComP[_ETE EL'L.L-1'pTlC I^TcGEtAl qF SECqnD k^MO 

CE(1 ) = lo570796 
CE(2 ) = lo566862 
CEO ) = lo5629l3 



CE ( 4 ) = 1 *558948 
CE< 5 1 = 1, 554969 
CE(6)=1o550973 
CF(7)=1o546963 
CE(8 )«1o542936 
CE(9 )“1u 538893 
CE ( 10 ) =1 , 5 348 33 
CE( 111=1,530758 
CE( 121=1,526665 
CE ( 13 ) =1 , 5225550 
E ( 141=1,518428 
C E ( 1 5 ) = 1 ,514234 
CE ( 16)«l o 5l0l22 
CF( 171=1, 505942 
CE ( 18 ) =1 0 5Q1743 
C E ( 1 9 ) “ 1 0 4 9 7 5 2 6 
CE (201=1,493290 
CE (21) =1,489035 
CE(22)=1„4847610 
CE(23)=1 o 480466 
CE(24)=1„476152 
CE ( 25 ) =1 ,471818 
CE(26)=1, 467462 
CE<27)=1, 463086 
CE(28)-1 0 458688 
CE(29)=1 0 454269 
CE(30)=1 o 449827 
CE(31)»1 0 445363 
CE(32)=1 0 44 0,376 
CE(33)=i, 4363-66 
CE(34)=1„431832 
CE(35)=1, 477274 
CE (36 ) =1 ,422691 
CE( 37) = 1 o 4180p, 3 
CE ( 33 ) =1 ,413450 
CE(39)=1, 408791 
CE(40)=l o 404l05 
CE(41)=1 0 399392 
C E ( 42 ) = 1 o 3 9465 2 
CE (43 1=1,389883 
CE(44)=1 o 385086 
CE (45)=1 o 380259 
CE (46 ) =1,375402 
CE (47)=l,3705l5 
CE(48)=l,365596 
C£(49)-l. 360645 
CE (505=1,355661 
CE(51 1 = 1,350644 
CE<52)=1,345592 
CE(53)=1,3405 o 5 
CE ( 54}=l ,335382 
CE(55)*1, 330223 
CE (561=1,325024 
CE(57)=1,319788 
CE ( 58 ) = 1 , 3 145 l 1 



CE ( 59 ) = 1 o 309192 
CE ( 60 ) = I o 30383 2 
CE(61)=1o298428 
CE(62)=I 0 292979 
CE(63)-1 0 2S7484 
CE(64)“1 0 281942 
CE (65 ) = 1 0 27635a 
CE (66)=1 ,270707 
C E ( 6 7 ) = 1 c 2 6 5 0 1 3 
CE(68)=1»259263 
CE(69)=1o 253458 
C E ( 7 0 ) = i s 2 A- 7 5 9 5 
CE(71)=1o241671 
CE(72)=1. 235684 
CE (7? 5 -I* 229632 
CE(74)=l 0 2235l2 
CE .( 75 >=1 0 217321 
0 (76) =1 o 2110 56 
CE ( 77 ) =1 0 2047 14 
CE ( 7S ) = 1 o 198290 
CE (79) =1. 191781 
CE(80)=1 o 185183 
CE( 81) =1.178490 
CE<82)=1 o171697 
CE ( 8 3 ) = 1 0 164798 
CE (84) =1 c 157787 
CE<85)=1 „ 150656 
CE ( 86 ) =1 0 143396 
CE(87)=i D l35998 
C E ( 8 8 ) = 1 o 128451 
CE ( 89 ) =1 o 15.5742 
CE(90)=l o 112356 
CE ( 9 1 ) = 1 ® 104775 
C E ( 9 2 ) = 1 o 0 9 6 4-7 8 
CE(93)=l on 37933 
CE(94)=l o 079l2l 
CE(95) = 1. 069.986 
CE(96)=1 o 060474 
CE(97)=1. 050502 
CE(98)=l or ,39947 
CE (99)=i an 28595 
CE ( 100 )-l 0 015994 
C 

GO TO 30 

3 IF (AKSQ-. 99)20,20 5 21 
21 P ARA=0 .25* (I. c-AKSQ ) 

700 TEST =L OOE-O 7 

IF (PARA-TEST 5 701,702,702 

701 PARA =TEST 

702 ZLP=AL0G ( 4o /PARA ) 

TK =ZLP*0 o 5* ( 1 o + PARA )— PA.RA 
TE=1 bO+(ZLP*PaRA)~PARA 
I GO TO 30 

|0 JA=100.0*AKSQ 

1 J A = 1 + J A 
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IFCCKKC JA)«AKSQ)22>23»22 
2’ 3 , ,,Tj<-CfCl JA) 

tE=CE(JA) 

GO TO 30 

22 C C N= (AK^Q-CKK ( JA ) ) / ( CKK ( JA+1 ) -CKK ( JA ) ) 
T<-CK ( JA ) +COi\!* { CK ( JA+1 ) — CK ( JA ) ) 

TE=CE( JAJ^C0N*(CE( JA+1 }-CE( JA) ) 

GO TO 30 

13 I F ( AKSQ— o a 1 ) ?2 1 j 72 1 3 72 q 

721 =ARA=o 25 *aSSQ 
GO TO 700 

720 AKSq»1 * — a K Z 7 

GO T020 

30 CONTINUE 
RETURN 
END 
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FAN 0 

u 

t 

NASA T M— D 8 8 - 

1 c 0 


0 o 8 

0 c 06 

0 


0 o o 

0 o 0 

v 0 A 


0.78 

.a 0 C 0 5 

-4 

10 

13 

1 1 

/, 

c H* 


0o45 

0 0 5 

0 b 


0 0 8 5 

0 • 9 

v. 3 3 9 


0 0 3 35 

0 0 33 

3 0 3 5 


-0 a 299 

0o 294 

I08 


38 0 5 

b3 

c 

CD 

O 

i O 

O 



2 7o 5 

w 0 9 


0 0 0 9 

Ov 09 

o09 


0 

0 

0 

^0 

0 c 0 9 

0 . 

1 

0 

0.54 

88.0 


1 o 006 

0 o 0 
0 o 8 

1 

a o 5 5 
0 o 95 
0 . 3.25 
0 o 289 
36 o 5 
26, .5 

0 ® 0 o, 


9 « 1 0 8 „.oi 


U 0 6 
loO 

Oo 321 

0 0 2 85 
35 a 0 
2 5 .- 5 , 

0 o 09 
0 o 0 9 


A -65 
Oc. 3 17 
33 c. 8 

0 o 09 


0 o 7 
0 o 3 12 
32 o 5 
O 0 O 9 


0 o 75 
0 « 307 
31 c 2 
O0O9 


Table 1: Input Data to Part I, Fan. • 
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64 64 10 



H4 

0*83 

0*94 

1 * 48 

l*0a 

2*02 

2 o 6 2 5 

1 o 2 5 

3*27 


1*43 

4*70 

5 ■■ 3 

1 = 43 

6*13 

7 :■ 0 1 

1*80 

7*93 


0*68 

u 0 98 

: ,46 

o o 9 7 

1 c 95 

5 -- 

1*07 

3 0 02 

= 'o 

1*36 

4*38 

:r 1 0 
- ■- i -1 

1*6° 

5*98 

-6 -78 

1*76 

7*74 

, 6 3 

0*6 

1*06 

1 ,53 

1*03 

2*09 

2 o 5 8 

1*04 

3 0 1 3 


1*25 

4*38 

5 s 1 

1 * 44 

5*82 

6 7 6 2 

1*72 

7*54 

= 7 5 

0*4 

1*04 

1 ,32 

0 o 63 

1*67 

1 o 97 

0 u 6 5 

2 0 32 

2*64 

0*70 

3 * 0'2 

3*4 

0 o oO 

3*92 

4 = 33 

0 * 9 0 

4*82 

5 ,32 

1*16 

5*98 

6 o 6 1 

1 o 38 

7*36 

, 94 

0 * 2 5 

1*06 

1 u 30 

0*48 

1*54 

1*74 

0*^6 

1 *90 

2*66 

0 * 84 

3*0? 

3*42 

0*86 

3*88 

4*31 

0 c 9’ 

4o 8 

5*36 

lol 

5 0 9 

6*53 

1*3 

7*2 

1HU4 

u ,“21 

1*17 

1*33 

0*39 

1*56 

1*8 

0*35 

1*90 

2 * 64 

0 * 83 

3*03 

3*44 

0 n 8 5 

3 = 88 

4*3° 

0 o 9 ? 

4o 0 

5 *28 

1*05 

5 0 8 5 

6 , .3 9 

1*18 

7 c 03 

1*24 

° * 2 7 

1 0 4 

1 o 6 1 

0 - 43 

1, 83 

2 c 09 

o*53 

2 u 36 

2*62 

o o 64 

3 - 0 

3*42 

u ^ 8 5 

3-85 

4*33 

0 o 05 

4o8 

5*27 

0 0 a 2 

5 0 7 2 

6*3 

I0I8 

60 9 


0 = 615 

1*23 

0 * 6 1 5 

0 = 615 

• 1 *23 

0*615 

0.615 

1*23 

0*615 

0*615 

1*23 

0*615 

u * 6 1 5 

1*23 

0*615 

0*615 

1*73 

0*615 

1*78 

1*10 

1*78 

1*78 

1,10 

1 = 78 

1,78 

1*10 

1*78 

1 *78 

1*10 

1*78 

O 

-4 

OO 

1*1 

1*78 

1*78 

1*1 

1*78 

2*82 

IcC 

2*82 

2=82 

1*0 

2*82 

2*82 

1*0 

2*82 

2*82 

1*0 

2*82 

2*82 

1*0 

2*82 

2*82 

1*0 

2*82 

3*8 

1*0 

3*8 

3*8 

1*0 

3 0 8 

3,8 

1.0 0 

3*8 

3*8 

1*0 

3*8 

3*8 

1*0 

3*8 

3*8 

1 s 0 

3*8 

3*8 

1*0 

3 0 8 

3*8 

1 = 0 

3*8 

4 = a 

1 0 0 

4*8 

4*8 

1*0 

4*8 

4*80 

1*0 

4* s0 

4*8 

1*0 

4*8 

4*8 

1*0 

4*8 

4*8 

1*0 

4*8 

4*8 

1 0 0 

4*8 

4*8 

1*0 

4*8 

5*8 

1*0 

5*8 

5 * 8 

1*0 

5*8 

5*8 

1*0 

5 = 8 

5 8 

1*0 

5*8 

3 <* B 

1*0 

5,8 

5 0 8 

1*0 

5*8 

5^8 

1*0 

5*8 

5 0 8 

1*0 

5 * 8 

608 

1*0 

6*8 

608 

1*0 

6*8 

6 & 8 

1*0 

6*8 

608 

1*0 

6*8 

608 

1*0 

6*8 

608 

1*0 

6,8 

6©8 

1*0 

6*8 

608 

1 *0 

6*8 


Data to Part II, Wing. 


Table 2 


Input 
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1 = 43 


0 o 33 

1 = 65 

7 o p6 

1*2 


7 = 86 

1 = ? 2 


0 o 6 7 

2 = 32 

7 = 86 

1 o 2 


7 o 86 

2 = 71 


0 = 88 

3 = 2 

7 = 86 

1 ° 2 


7 = 86 

3=66 


1 = 05 

4a 25 

7 ,,86 

1*2 


7 o 86 

4 c 78 


1 o 0 5 

5 o 3 

7.86 

1 o 2 


7 o 86 

5 a 78 


1 = 4 

6 = 7 

7 = 86 

1*2 


7 o 8 6 

1 o 66 


0 o 32 

1=85 

° o 17 

1 o 4 


-7 = 17 

2 = 28 


1 = 15 

3 o 

° = 17 

1 o 4 


9 = 17 

3 = 76 


1 = 55 

4o 5 5 

9 o 17 

1=4 


9 = 17 

5 o 44 


1 o 95 

6,5 

9.17 

1.4 


9.17 

1 o 98 


0.38 

2 o 2 

I0o9 

2 o 0 8 


1C o 9 

2 o 58 


0 = 95 

3 = 15 

1 0 o 9 

2 = ug 


1 0 o 9 

3=76 


1 ° 3 5 

4 = 5 

10 = 9 

2 = 08 


1 0 o q 

5 = 26 


1 = 68 

6=18 

1 0 o 9 

2 = 08 


10 = 9 

7 

7 

7 Q 

9 9 

? 7 

5 

5 


1 

7 

13 19 

27 35 

' 43 51 

57 

6 1 
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48 64 10 


n 4 

0.83 

0.94 


..615 

1.23 

00 

o 

r — 

loQa 

2 o 02 


0.615 

1.23 

2.62 5 

1.25 

3.27 . 


0,615 

1 e 23 

3 ,96 

1 o 43 

4o70 


0.615 

1.23 


1 .43 

6.1.3 


0.615 

1 ° 2 3 

7 , 0 1 

1.80 

7,93 


0.615 

1.23 

o 5 M- 

0 o 68 

0.98 


1,78 

1.10 

1.46 

0,9 7 

1.95 


1.78 

1.10 

2 o 52 

1 . C 7 

3.02 


14 78 

1.10 

3 o 6 9 

1.36 

4.38 


1.78 

1,10 

5,12 

1.6 a 

5.98 


1.78 

1.1 

6 c78 

1.76 

7.74 


1.78 

1.1 

.65 

0.6 

1.06 


2.82 

1.0 

1 o 5 3 

1 o 03 

2.09 


2 o 8 2 

l.o 

2 o 5 3 

1.04 

3.13 


2 o 8 2 

1.0 

3 , 7 

1.25 

4.38 


2 o 82 

l.o 

5.1 • 1.44 

1^30 0^48 

5.82 2,82 l.° 

7.54 2.82 1.0 

1 .04 5.8 1 . 

1 o 06 4.8 1.0 

1.54 4.8 1.0 


1 o 74 
6 o 5 3 
1 o 04 
1.33 
1.8 
6 „ 39 
1 o 24 

1 o 6 1 

2 0 09 

6„3 

1 0 43 
1,92 

2 o 7 1 
3.66 

4,78 
5.78 
1.66 
2 .,2 8 
3.76 
5 o 44 
1.98 

2 o 58 

3 J6 
5 o 26 


0,36 
5 o 02 

0 o 2 1 

0 o 39 
0,35 
4.83 
0,27 

0 o 43 
0,53 

4 o 54 
0,33 
So 67 

0 o 88 

1 o 05 
1.05 

1 o 4 
0,32 
1,15 
1,55 
1 o 9 5 
0 o 38 
0 o 9 5 
1,35 

1 o 68 


1 o 90 
7 o 2 
1 o 17 
1,56 
1,90 
7,03 
1 © 4 

1.8 3 
2,36 

6.9 
1,65 

2,31 
3 o 2 
4,25 
5.3 
6,7 
1 o 8 5 
3,0 
4„ 55 
6 o 5 
2.2 
3 o 15 
4 0 5 
6.18 


■4* a ■* 

4® 8 
5.8 
5 08 
5 o 8 
5,3 
6,8 
6.8 
6,8 
6 o 8. 

' o 8 6 . 
,86 
,86 
7 o 86 
7,86 
7 o 86 

9.17 
9 e 17 

9.17 
9 o 17 
l«o9 
l«c 9 
l rt< >9 

l.c9 


i 

7, 

7, 


lo0 
1.0 
IcO 
1.0 
1.0 
loO 
loO 
loO 
1.0 
loO 
1.2 
1.2 
1.2 
1.2 
1.2 
1.2 
1.4 
1.4 
1.4 
1.4 
2. .8 
2,„8 
2 0 08 
2.08 


Table 3: Input Data to Part 2 ,- Wing 


0 o 6 15 
0.615 

0. 615 
0.615 
0,615 
0.615 

1,78 

1,78 

1. -78 

1 o 78 
1 o 78 
1.78 

2 o 82 

2 o 82 
2 o 82 
2 o 82 
2 o 8 2 

2 o 82 

3 o 8 
3 o 8 
3.8 

3.8 
4.8 
4.8 

4.80 
4.8 
5.8 
5.8 
5 o 8 
5.8 

6 o 8 

6.8 
6.8 
6.8 

7.86 
7.86 
7.86 
7.36 
7.86 
7.86 
9 o 17 
9.17 
9.17 
9.17 
10.9 
10.9 
1C. 9 
10.9 


with a duct. 
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3 ,02 

3 o92 
4 C 5 
5 ,0 

3 o 02 
3,83 

4 s 8 
3 ,9 
3 c 03 
3,83 



3 o 0 

3 o 8 5 

4 c 5 
5,0 

7 7 

1 7 

7159675 


3 “ 8 
3 o 8 
3-8 

3 o 8 
. 4o 8 

4 ° 8 
4,8 

4 •> 8 

5 o 8 

5 o 3 
5 c. 8 

5 -j 8 

6 ° 8 
6“8 
6,3 
6 o 3 

75555755 
13 19 23 27 31 35 4i 45 

0 o 63 9485 ?“*0 <>292423.4 0o I 8 QO 15 O-O 0 I 35 O 915 


0o 1013527 




i mm i 

' .v ; , ✓ \ 

■ I pj$ip I 

«|iM j pi a^ ailgii 


HHSGI 









SECTION ON A A 
ELEVATION 

all Dimension are in inches 
unless otherwise noted 



CHARACHTERISTiCS OF THE MODEL 










S33U030 N! & 


ix? -<r cn o 

oo pq m ro O 

S3HDNI N! GdOHJ 


FIG. 7_ BLADE FORM CURVE 
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REFERENCE CYLINDER 




FIG. 10 _ CENTRE BODY GEOMETRY 




iu L _ NOTATION USED IN DERIVING THE VELOCITY 
COMPONENTS INDUCED BY A VORTEX RING 
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x 

x b±_.L_, 


"■ . '■ ■ 
l ■ 

U*' 

■' ( 

X 


r' - 




0 ~ ' 

n ; 


1 

1 


DUCT REFRENCE CYLINDER 


FIG.13_ NOTATION USED IN DERIVING THE INDUCED 
VELOCITY COMPONENTS DUE TO THE DUCT 
BOUND VORTICITY 



-'start) 


MAIN 

_J5UQO_" 


Eli ips 


f lambda! 

1 SU 04 { 


i *T PKL ■] , 

I SU 02 I HUB 

f *! SU 05 

i r c l a i“i 

SU 08 ! . _aitrng 

j SU 1 4 

INPUT - ] '“PRESS' 

1 SUOI j H" "1 SU 19 


PKL 
SU 02 


J , p AMFFP 1 
LULQ£l_J 




HUB I FOURCS 

su o s r r su n 

"prop } 


GAMCVL 
SU 15 


r BNC'dEF|_ J FLOPS' 
" su 16 ! lSUOS 



! e Clips - ! 
i _SU 03 J 

[Tourcs ! 
! sun 


. ELLIPS ! 

il i SUQ3 1 

“ARC'S IN - ] 

I SU 09 i 

1 [lambda! 

^ — *1 SUQ4 j 


SRCRNGL- ELlIF-S 
SU 12 ; | SUQ3 


syj05j 


t_ J ANCOEF: 

0 SU 17 p ~T~ 


I ELLIES 

1 SU0 3 j 

r FOURCS! 

1 sun ! 


i SU 09 j 

[“lambda 

1 SU 04 f 


(Q) INITIALIZATION SECTION 


FIG. !5_ RELATIONSHIP BETWEEN SUBPROGRAMS 
OF PART I 










A 



(b) ITERATION AND CONCLUSION SECTION 
FIG. I5_ CONCLUDED 












FIG. 16_ RELATION SHIP BETWEEN SUBPROGRAMS 
OF PART EL 









I 



FIG. 17_ VARIATION OF WING AND WING WITH A 
DUCT C L WITH WING ANGLE OF ATTACK 



THRUST ib. 






F1G.19 -VARIATION OF FAN - IN - WING TOTAL LIFT COFFFiCIENT WITH 
WING ANGLE OF ATTACK. 


SYMBOLS DENOTE 




o^wda-1 x=Wddn 

b- 0 


RE MENTAL LIE T COE EFICIENT . WITH PROPELLEF 



C f =0.423 

q = 4j0 P S 
00 

rpm. =30 

a ® © J = 0.35 

C' F = 1.64 8 

q 00 =4.IPSF 

rpm = 6000 

y v = 0.35.3 a © 

0.6 • ® 

0-98 X a 



U- FAN 4 

N INDUCED UP WASH DISTRIBUTION ON T 




a. 

u 



